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Modern laser sources nowadays deliver ultrashort light pulses reaching few cycles in duration, high 
energies beyond the Joule level and peak powers exceeding several terawatt (TW). When such 
pulses propagate through optically-transparent media, they first self-focus in space and grow in 
intensity, until they generate a tenuous plasma by photo-ionization. For free electron densities and 
beam intensities below their breakdown limits, these pulses evolve as self-guided objects, result- 
ing from successive equilibria between the Kerr focusing process, the chromatic dispersion of the 
medium, and the defocusing action of the electron plasma. Discovered one decade ago, this self- 
channeling mechanism reveals a new physics, widely extending the frontiers of nonlinear optics. 
Implications include long-distance propagation of TW beams in the atmosphere, supercontinuum 
emission, pulse shortening as well as high-order harmonic generation. This review presents the 
landmarks of the 10-odd-year progress in this field. Particular emphasis is laid to the theoretical 
modeling of the propagation equations, whose physical ingredients are discussed from numerical 
simulations. The dynamics of single filaments created over laboratory scales in various materials 
such as noble gases, liquids and dielectrics reveal new perspectives in pulse shortening techniques. 
Differences between femtosecond pulses propagating in gaseous or condensed materials are un- 
derlined. Attention is also paid to the multifilamcntation instability of broad, powerful beams, 
breaking up the energy distribution into small-scale cells along the optical path. The robustness 
of the resulting filaments in adverse weathers, their large conical emission exploited for multipol- 
lutant remote sensing, nonlinear spectroscopy, and the possibility to guide electric discharges in 
air are finally addressed on the basis of experimental results. 
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I. INTRODUCTION 

Over the past two decades, ultrafast laser sources pro- 
ducing ultrashort pulses have come of age. Technological 
advances in this field have permitted the generation of 
light wave packets comprising only a few oscillation cy- 
cles of the electric field. In space, the extent of the pulses 
becomes focusable to a spot size comparable to the laser 
wavelength. In time, mode-locking and chirped-pulse 
amplification (CPA) technologies allow to access smaller 
and smaller durations and optical intensities locally ex- 
ceeding hundreds of terawatt (TW) per cm 2 at moderate 
(sub-mJ) energies. Power levels nowadays approach the 
petawatt range with femtosecond pulse durations (I fs 
= 1CP 15 sec). Ultrashort light pulses enable researchers 
to probe ultrafast relaxation processes on never-before- 
accessed time scales and study light-matter interactions 
at unprecedented intensity levels. Due to extreme tem- 
poral and spatial confinements, the pulse strength ex- 
ceeds that of the Coulomb field which binds electrons 
at their nucleus. It becomes strong enough to overcome 
the Coulomb barrier and triggers optical-field ionization. 
Pushed to ultrahigh intensities > 10 18 cm 2 , the availabil- 
ity of CPA lasers has extended the horizon of laser physics 
from atomic and optics studies to relativistic plasm as, 
nuclear and high-energy physics ( Mourou et q/J . 120061 ) . 

Before reaching such extreme intensities, progress in 
ultrashort laser systems has led to the ability to ob- 
serve tunnel or multiphoton ionization in several spectral 
ranges, before the ionization process reaches saturation. 
In this regime, a n intriguing phenom enon was discovered 
by Braun et al. ([Braun et all , Il995l ) in the middle of the 
nineties. By launching infrared pulses with femtosecond 
durations and gigawatt (GW) powers in the atmosphere, 
the beam became confined to a long-living, self-confined 
tube of light capable of covering several tens of meters, 
i.e., many linear diffraction lengths, along the propaga- 
tion axis. The mechanism supporting this "light bullet" 



FIG. 1 Principle of producing femtosecond filaments in air at 
the laser wavelength of 800 nm. Photograph at the right-hand 
side shows a transverse cut of the filament profile. 



results from the balance between Kerr focusing, which 
increases the local optical index with the wave intensity, 
and self-induced ionization. When an ultrashort pulse 
self-focuses and couples with a self-induced plasma chan- 
nel, its spatial profile exhibits a narrow extent along the 
optical path. Spectra broaden due to self-phase modu- 
lation (SPM), which is sustained by the mechanism of 
high-intensity clamping. This picture classically refers 
to what is commonly called a "femtosecond filament". 
In the diffraction plane, this filament is characterized by 
a white- light spot, surrounded by concentric "rainbows" 
with colors ranging from red to green. The high nonlin- 
earities competing through the filamcntation process pro- 
duce an impressive supercontinuum leading to white-light 
emission. They also affect the beam divergence through 
an apparent conical emission, as illustrated in Fig. [T] 

Femtosecond filaments have, for the last decade, 
opened the route to a fascinating physics. Their highly 
nonlinear dynamics sparked broad interest, first because 
of the capability of femtosecond pulses to convey high 
intensities over spectacular distances, second because of 
the white light emitted by the filame nts transforming in- 
frared lasers into "whit e- light lasers" ( Chin et 'all , [l 999b: 
iKasparian et all . l2003lh For appropriate beam configu- 
rations, long fs filaments can be created not only in the 
atmosphere, but also in noble gases, liquids and dielec- 
tric solids, as long as the pulse intensity does not reach 
the limit of optical breakdown, which implies the electron 
plasma density to remain at subdense levels. These opti- 
cal structures are subject to strong modifications of their 
temporal profile, triggered by SPM together with the ion- 
ization front and the chromatic dispersion of the medium. 
The dynamical balance between nonlinear focusing and 
ionization can result in a drastic shortening of the pulse 
duration, down to the optical cycle limit. This prop- 
erty opens quite promising ways to generically deliver 
light with durations of a few fs only, which should fur- 
ther impact the fields of high-order harmonic generation 
and sub-fs pulse production. Besides, femtosecond pulses 
with broad spatial extents create several filaments, whose 
mutual interactions support the self-guiding of the beam 
envelope and can preserve a confined state upon several 
kilometers. This recently led to develop ultrashort Light 
Detection And Ranging (Lidar) facilities, that exploit the 
white light emitted by these filaments, in order to detect 
and identify several pollutants within a single laser shot. 

Many applications have been inspired by this ultra- 
fast "light bullet", which justifies the present review. 
Before commenting on these, an accurate understand- 
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ing of the filamentation phenomenon requires a rigorous 
derivation of the propagation equations together with ba- 
sic tools for capturing the underlying physics. For this 
purpose, Section [IT] addresses the model describing the 
long range propagation of ultrashort laser pulses in any 
optically-transparent medium. Section IIIII reviews the 
basic phenomenon of wave self-focusing and its limita- 
tion by potential players such as plasma generation, chro- 
matic dispersion and optical nonlinearities. Emphasis is 
given to semi-analytic (e.g., variational) methods provid- 
ing qualitative information about these effects. Section 
IIVI lists the major phenomena driving femtosecond fila- 
ments, whatever the propagation medium may be. Sec- 
tion|V]is devoted to pulse shortening that can be achieved 
by letting femtosecond pulses ionize atom gases at ap- 
propriate pressures and to high-order harmonic genera- 
tion. Section IVII addresses different propagation regimes 
in condensed materials (water, silica glasses), from the 
laser-induced breakdown limit to Wshapcd nonlinear 
waves self-guided with no plasma generation. Section 
IVIII concentrates on the atmospheric applications of the 
white-light supercontinuum emitted by ultrashort fila- 
ments. Current techniques for changing their onset dis- 
tance and self-channeling length are discussed, together 
with the diagnostics used for plasma and optical-field 
measurements. Attention is paid to novel ultrashort 
LIDAR-based setups and their use in remotely analyzing 
aerosols, biological agents, dense targets through remote 
filament-induced breakdown spectroscopy. Their ability 
to trigger and guide electric discharges over several me- 
ters is also discussed. Section IVIIII finally summarizes 
the principal aspects of this review and presents future 
prospects. 

II. PROPAGATION EQUATIONS 

To start with, we derive the model describing the prop- 
agation of ultrashort optical pulses in transparent me- 
dia. Using conventional description of nonlinear optics, 
straigh tforward combination of th e Maxwell's equati ons 
yields (|Agrawall 120011 : iHe and Liul . Il999t IShenl . Il984fl 

V 2 E - V(V • E) - cr 2 d 2 t E = p {d 2 P + 9t J), (la) 
V -E = (p-V -P)/e a , (lb) 

where eo, /io and c denote the electric permittivity, mag- 
netic permeability and the speed of light in vacuum, re- 
spectively. The optical electric field E, the polarization 
vector P, the carrier density p and the current density 
J are real valued. For further convenience, we intro- 
duce standard Fourier transforms applied to the fields 
(E, P, J)t as 

(S, P, Jj(f, u) = ^J (E, P, J)\r, t)e^dt. (2) 

The current density J describes the motions of the free 
electrons created by ionization of the ambient atoms, 



for which ion dynamics is discarded. The polarization 
vector P describes the bounded electron response driven 
by the laser radiation. It is usually decomposed into a 
linear part Pl = PW related to the first-order suscep- 
tibility tensor x and a nonlinear one Pnl satisfying 
\P^\ ^ \Pnl\- For isotropic, homogeneous, non magne- 
tizable media and spectral ranges far from any material 
resonance, P can be expressed as a power series in E: 

P = P (1) (r» + P (3) (f» + P (5) (r» + . . . (3) 
with scalar components 

= £ o - X^l 1 ... aj {- UJ <r^i,---,^j) 

x E ai (r, lui) . . . E a] (r, lu 3 )S(lu - cj a )dwi . . . dtOj, 

(4) 

where tu a = Wi + + All susceptibility tensors x 
with even index j vanish due to inversion symmetry. The 
subscript p indicates the field vector component in Carte- 
sian coordinates and the indices <x,- have to be summed 
up over x, y, and z. The tensor x is diagonal with 

xjS = x (1) <W so that 

£«(f )W ) = e X (1 ^)£(^), (5) 
and the scalar dielectric function, defined by 

cH = i + x w M, (6) 

enters the wave number of the electromagnetic field 
k{uS) = \J e(uj)uj/c. Since x^ is complex- valued, the 
dielectric function e(u>) contains every information not 
only about the material dispersion, but also about the 
linear losses given by the imaginary part of X 
When losses are negligible, e(ui) is real and reduces to 
e(u>) — 7j 2 (w), where n(u>) here denotes the linear re- 
fractive index of the medium, which can be described in 
certain frequency ranges (far from resonances) by, e.g., 
a Scllmcicr formula. By convention, hq = n(u>o) for the 
central frequency luq = 2irc/\o of a laser operating at the 
wavelength A and fc = k(u>o). Without any specifica- 
tion, lo — 2nc/ A, w± — 2it/k± is the waist of the optical 
wave packet in the plane (x, y) and k± is the correspond- 
ing extent in the transverse Fourier space. 

A. Helmoltz equation 

By taking its Fourier transform, Eqs. (jTaJ) expresses as 
the Helmoltz equation 

[d 2 z + fc 2 H + Vi] E = -aw 2 ^nl + V(V • S), (7a) 
V-E= (e e) _1 (p-V-P N L), (7b) 
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where 



d 2 + dt stands for transverse diffraction and 



F Nh = Pnl + iJ/u = F®E, 



(8) 



gathers all nonlinear contributions through the function 
F that depends on E and can be viewed as the effec- 
tive nonlinear refractive index change of the medium (<g> 
is the convolution operator). Because Eqs. Q are usu- 
ally difficult to integrate in the full space-time domain, 
assumptions are requested to simplify them into more 
tractable form. The most fundamental of those consists 
in supposing that the wavefield keeps a transverse exten- 
sion always fulfilling 



k 2 Jk 2 {u) < 1, 



(9) 



i.e., for k(u>) located around wq, the transverse waist of 
the beam has dimensions larger than the central wave- 
length. The second one assumes small nonlinearities, i.e., 



e e(w) 



< 1. 



(10) 



1. From vectorial to scalar description 

The previous conditions make vectorial effects negli- 
gible for, e.g, a transversally-polarized light field E± = 
(E x ,E y ). Indeed, Eqs. (J7J) can be combined, with the 

help of the continuity equation dtp + V ■ J = expressed 
in Fourier variable, into the form 



[dl + k 2 {uo) + \7\]E = -p u J 2 



NL 



V(V ■ -Fnl) 
fc 2 (w) 



(11) 

whose last term scrambles nonlinear vectorial compo- 
nents. When we project the vectors E = (E±,E z y, 

p NL = (P^P^V, J = (J±,J z y, v = (v^A)t 

onto the transverse and longitudinal axis with unit vec- 
tors ej_ and e z , respectively, E z is found to scale as 
0(k±/k). This follows from a direct Fourier trans- 
form of Eq. (|7b[) for weak nonlinearities [Eq. (TIT))) ]. 
Expressed in Fourier space, the nonlinear coupling of 
transverse/longitudinal components d escribed by the 
scram bli ng term behave as Q(k 2 / k 2 ) ( Fibich and Ilanl . 



2001a. a; iMilsted Jr. and Cantrell 119961) So, these ef- 
fects become important in the limit k± — » k(u>) only. 
Reversely, if we postpone to further demonstrations that 
nonlinear compression processes are stopped before k± 
becomes comparable with k (by, e.g., chromatic dis- 
persion or plasma generation), the last term in the 
right-hand side (RHS) of pip is close to zero, implying 
thereby E z ~ 0. The field remains transversally polar- 
ized along the propagation axis, making the influence of 
V(V-^nl)/^ 2 negligible. Hence, as long as the nonlinear 
polarization and current density preserve the conditions 
© and IjlOp. vectorial effects can be ignored for purely 



optical or weakly-ionized materials as well. This prop- 
erty justifies the use of a scalar description for linearly- 
polarized beams having, e.g., E y = 0. 



2. Weak backscattering 

The question of evaluating backscattering waves may 
be cruci al in several area s, such as remote sensin g exper- 
iments ( Kasparian et all 120031 : lYu et al. , 2001), which 
spectrally analyze the photons returning towards the 
laser source. The amount of backscattered photons, how- 
ever, constitutes a weak percentage of those traveling 
in the forward direction through a transparent medium. 
The reason can be seen from the scalar version of Eq. 
(fTTj) expressed as 



D + (u)D-(uj)E = -V\E - p LO 2 F(g> E, (12) 

where D±(uj) = d z =F ik(u>). By substituting the solution 
E = U+c ik ^ z + U-e~ ik ^ z , Eq. CCH) expands as 



JKk(u)*[gj2 + 2ik(Lu)d z + Vl + p u 2 F]U+ 

+ [d 2 - 2ik(u)d z + Vi + p u 2 F]U- = 0. 



(13) 



Here, U+ and U- represent the Fourier components of 
the forward and backward running fields, for which we 
a priori assume |d z L/±| <C \k(u>)U±\ and U+ ^> U-. For 
technical conv enience, we assume F = FS(u>). Following 
Fibich et al. ( Fibich et all [20021 ) . we can integrate Eq. 
(fT3"|) over the interval z — it /2k < z < z + ■n/2k (one fast 
oscillation) and Taylorize U±(z) to evaluate 



2ik(co)d z U- 



-d z [V\+p uj 2 F]U^ 



(14) 



— k\ in Fourier space, the backscattered 



Since 

component has a weak influence on the beam dynamics 
if k\ <C k 2 {w) and as long as the longitudinal variations 
of the nonlinearities remain small. 



3. Unidirectional pulse propagation 

The limit §§§ moreover implies that the wave compo- 
nents forming the angle 9 = arcsin (k±/k) between the 
transverse and longitudi nal directions mostly propagate 
forwards since 9 <C ir/2 ( Feit and Fleckl H988). Because 
the propagation physics is mainly brought by the for- 
ward component, one has U- — > and E ~ U + e lk ^ z . 
With the above approximations, the operator D_(w) for 
backscattering mainly applies to the most rapid varia- 
tions of the field, expected to be driven by the complex 
exponential, i.e., D_(gj)E ~ 2ik{uS)E. The so-called Uni- 
directional Pulse Propagation Equation (UPPE) 



d z E 



2k(cu 



■V\E + ik(uj)E - 



ip uj 
2k{uo) 



XL 



(15) 
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then naturally emerges from Eq. (|12|) . Validity 
of this model explicitly requires that the second- 
order derivative in z of the envelope function 
U+ must be small compared with \k{u))d z U+\ 1 

since D+MD _M£ = e i fc <"Hd 2 + 2ifcM&lL/4 . 

dGeissler et all Il999l ; iHusakou and Herrmann! . 
2001). This approximation, usually expressed as 
|9 2 C/+| <C \k(uj)U+\, refers to the " paraxiality" as- 
sumption. It holds if the field envelope U+ does not 
significantly change over propagation distances of the 
order of A, for all wavelengths under consideration. 
Paraxiality is again linked to the weakness of both the 
ratio k±/k and the nonlinearities. Let us indeed assume 
the nonlinear function F clamped at a maximal constant 
level, F max . The forward component of Eq. (fT3|) then 
goes like 



(16) 



It is seen right away that this solution fits that of the 
paraxial model discarding second derivatives in z, i.e., 

U+ ~ c -i(k 2 ± -vou> 2 F m ^)z/2k(u J ) j ag long ag the twQ CQn _ 

straints © and (fit))) apply. 

Recently, K olcsik et al. ( Kolesik and Molonevl 120043 
iKolesik et all l2002t) proposed an UPPE model based 
on the exact linear modes of the Helmoltz equation. 
The key idea is to preserve the linear solutions, En n = 
e ±i^/fc 2 (aj)+v^z^ ag j on g ag p 0SS i^i e along the manipula- 
tion of Eq. ([7a|l . Rewriting Eq. fH} as 



(17) 



with D±(oj) — (d z =F iyjk 2 (ui) + Vj_), we may retain 
the forward running component E = ^e'V^^'+^i 2 
constrained to the "paraxial" limit D^(oj)E « 

2iyjk 2 {uj) + V\E. Repeating the previous procedure 
immediately yields 



d z E = i x lk 2 {oj) + V 2 L E- 



2 v /fc2(^ 



(18) 



Eq. (fT5|) allows to formally describe DC-field components 
(u) = 0), whereas Eq. (TT5|) is limited to non-zero fre- 
quencies strictly. Despite these minor differences, UPPE 
models (fT5|) and (fl~8|) become quite analogous when the 
condition k\jk 2 <C 1 applies. The major advantage of 
the UPPE models is to elude the formal use of a cen- 
tral optical frequency and correctly describe the complete 
spectrum of pulses in nonlinear regimes, even when they 
develop very large band widths. 

For practical use, it is convenient to introduce the com- 
plex version of the electric field 



E = yfc{£ + E*), £ = 



1 



9(w)Ee-™du, (19) 



where c\ = o;o/^o/2fco and 0(x) denotes the Heaviside 
function. Because £ satisfies £*(uj) = £(—ui)* (* means 



complex conjugate), it is then sufficient to treat the 
UPPE model (TT5)) in the frequency domain uj > only. 
The field intensity can be defined by E 2 averaged over 
an optical period at least, for a given central frequency 
luq. This quantity usually follows from the modulus of 
the time averaged Poynting vector. It is expressed in 
W/cm 2 and with the above normalization factor c\ it is 
simply given by the classical relation / = \£ | 2 . 



4. Envelope description 

When a central frequency ujq is imposed, Eq. (TT5)) resti- 
tutes the Nonlinear Envelope Equation (NEE), earlier de - 
rived by Brabec and Krausz ( Brabec and Krauszlll997h . 
We can make use of the Taylor expansion 



k(u) = k + k'uj + V, V = V —w" 



(20) 



n>2 



where u> = lo — ujq, k 1 = <9fc/9cj| w=W[) 
d n k/dio n \ ul=WQ and develop Eq. flTSJ) as 



and fe(") = 



d z E 



iV 2 

—±-+i{k + k'oj + V) 
2k(uj) 



x E(uj)e' lult duj + i^ 



k{u) 



J r NL(w)e" 



t duj. 
(21) 



We apply the property following which E(u>) is the 
Fourier transform of i?(i)e 4 " * in u>, so that Q corresponds 
to idt by inverse Fourier transform. Terms with k(oj) in 
their denominator are expanded up to first order in ui 
only. Furthermore, we introduce the complex-field rep- 
resentation 



(22) 



involving the novel envelope function U . Next, the new 
time variable t — > t — z/v g can be utilized to replace 
the pulse into the frame moving with the group velocity 



= k r 



Eq. ([21]) then restores the NEE model 



{id z +V)U 



where 



2k Q 



(ViC/) 



2k Jcl 



TT^(U), (23) 



+00 



V = y2(k {n) /n!)(id t ) n , T=(l + —d t ), (24) 



whenever \ko — wnfc'l/^o <C 1. This condition is met if 
the difference between group and phase velocity relative 
to the latter is small, which is fulfilled in a wide range 
of propagation phenomena. The operator T _1 intro- 
duces space-time focusing in front of the diffraction term 
[T _1 (V 2 L C/)]. On the other hand, nonlinearities with the 
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envelope function ^^(U) are also affected by the op- 
erator T, which refers to self- steepening. For dispersion 
relations truncated at some finite orders n < +00, the 
NEE model applies to optical fields with sufficiently nar- 
row spectral bandwidths. With full chromatic dispersion, 
it holds for describing light pulses down to the single cy- 
cle. 



B. Nonlinear optical responses 

We henceforth assume a linearly polarized field (along, 
e.g., e x ) and treat nonlinear effects within a scalar de- 
scription. For centro-symmetric materials, only one rel- 
evant component of the tensor remains in the cubic con- 



tribution P^ 3 ), e.g., x (3) = XxLx (|Agrawall[200ll ). For 
simplicity, we may first consider \^ as keeping a con- 
stant value for a spectral domain centered around ujq. Eq. 
then simplifies with a single component, noted xL , 

and in time domain one finds P^(f, t) = eoxS^ 3 - This 
expression holds whenever we suppose an instantaneous 
response of the medium, which ignores the contribution 
of molecular vibrations and rotations to X ■ Strictly 
speaking, however, the phenomenon of Raman scatter- 
ing comes into play when the laser field interacts with 
anisotropic molecules. This interaction can be schema- 
tized by a three-level system built from the rotational 
states of a molecule. The molecular scatterer has two 
rotational eigenstates, the ground state (level 1) with en- 
ergy Mli and an excited one (level 2) with energy hVl 2 
where h — 1.06 x 1CP 34 Js and £l m denotes the frequency 
of the state m = 1,2,3. Far above lies an electronic 
(or translational) state with energy Ml 3 ^ Ml 2 — Ml\. 
This molecule interacts with the laser field whose pho- 
ton frequency ujq fulfills fli3,fl 2 3 ^ ujq ^> fl 2 i, [il n m = 
— ^m], such that state |3) cannot be populated. Be- 
cause of the definite parity of these molecular states, the 
dipole matrix element fi\ 2 associated with the transition 
|1) — * |2) via a single photon is null, so that the rotational 
state 1 2) can only be excited via transition through a vir- 
tual state |3) [/ii3 ~ ^23 = A* 7^ 0]. Following this path, a 
Stokes photon with energy hw s — Hujq — Ml 2 i is emitted 
and the corresponding polarization vector involves the 
density matrix element associated with the states |1) and 
|2) as P Raman = x (1) [Pi2e tujRt + c.c.]E. Here, uj r = 21 
is the fun damental rotat i onal f requency and p± 2 is found 
to satisfy ( Peiiano et all [2003h 



P12 

T2 



,,2 Z?2 

j » K c -it/n 



(25) 



where t\ — \/lor and t 2 is the dipole dephasing time. 
Eq. (j2"5|) provides the Raman response 



P Ra 



2 X ( V 



E j e-^fsm(- — -)E 2 (t')dt', 
which originates from nonresonant, nonlinear couplings. 



Expressed in terms of the rescaled complex field 
£ [Eq. (|19Dl and with appropriate normalizations 



( Sprangle et all l2002h . it completes the cubic polariza- 
tion as 



+ OC 



P0) = 2n Q n 2 e ^/~cT / R(t ~ t')\£(t')\ 2 dt'£ 



(27a) 

+ 2n n 2 £o\fc[{l - x K )£ 3 /3 + c.c, 
R(t) = (1 - x K )5{t) + x K Q(t)h(t), (27b) 



h(t) = - 



2 7* 



3 tit 2 2 



l e -t/r 2 



sin(t/ri), 



(27c) 



with the definition of the nonlinear refractive index 
n 2 = 3xL 3 '/(4^o C£ o)- Here, contributions in 0(£ 3 ) are 
retained to further describe third-harmonic generation. 
Expression (|27a[) possesses both retarded and instanta- 
neous components in the ratio xk- The instantaneous 
part ~ S(t) describes the response from the bound elec- 
trons upon a few femtoseconds or less. The retarded 
part ~ h(t) accounts for nuclear responses, namely, the 
Raman contribution, in which fast oscillations in E 2 give 
negligible contributions, as T\ and t 2 currently far exceed 
the optical period ~ . 

The fraction of delayed Kerr depends on the molec- 
ular speci es under consideratio n. For air at 800 nm, 
Sprangle ( Sprangle et all [2002;) suggests ti ~ 62 fs, 
T2 ~ 77 fs and xk = 1/2. This choice is consistent with 



that proposed in experime ntal papers (|Nibbering et all , 
119971 : iRipoche et aU . Il997t) . When n ~ r 2 , the func- 
tion h{t) ~ (l/T\)e~ t l Tl can a lso be used in the ratio 
xk = 1/2 (jChiron et all fl999l ) . For condensed materi- 
als, the parameter ranges t 2 /ti — 2 — 4, t 2 = 30 — 50 
fs with xk = 0.15 — 0.18 h ave been suggested ( Agrawall 
12001b IZozulva et aUll999[) . Values of the nonlinear Kerr 
index n 2 can be f ound i n the literature jGong et al. 1 199S ; 



Hcllwarth et all fl990: 



Lehmeier et all 1985: Luo et 



119951 : iNibbering et al 
cm 2 /W for gases and 10 



it) 



19971 ). Comprised between 10 
16 cm 2 /W in dense media, they 



may, however, vary by a factor of the order of unity, de- 
pending on the procedure used for their evaluation (po- 
larization spectroscopy, self- or cross-phase modulated 
spectra, time-resolved interferometry) , together with the 
laser wavelength and pulse durations at which measure- 
ments are performed. 

Besides, the susceptibility tensor ha s nonlinear compo- 
nents x^ >3 ^ that satisfy the ordering ( Bovdl . [l992l : ISheiil . 
I1984T) 



p(fc+2) x (k+2) E k+2 \ E Vi 



p(k) 



(fc) 



P* 



Pa 



(28) 



where P a t — 3 x 10 10 V/m is the characteristic atomic 
electric field strength with intensity I a t > 10 14 W/cm 2 . 
Typically, the evaluation x^/x^ ~ 10~ 12 holds for 
nonresonant interactions in, e.g., ga ses. Desp i te th e 
lack of knowledge on the sign 

ofx (5) dPan et oZ.I . fl990h . 
the quintic susceptibility is often expected to satu- 
rate Kerr focusing and has, therefore, a negative sign 
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( Nurhuda et "aLl. l2002bri . Since the ordering (|28|) suggests 
that x i s rapidly decreasing with the order j, the Tay- 
lor series with respect to the electric field is truncated at 
the 5th order. Quintic polarization can be derived fol- 
lowing the same procedure as above, with \^ assumed 
constant in the frequency domain. By developing E 5 in 
terms of (£,£*) [Eq. (|19pj. the quintic contribution of the 
polarization vector then expands as 

P< 5 > = -2n n 4 e ^r + \ \£? ? + ^) £ + cc. 

(29) 

where n± = b\x^o |/ (4rigC 2 e§) . On the whole, the total 
nonlinear polarization vector reads as 



NL 



p(3) i p(5)^ 



(30) 



C. Plasma generation for singly-charged ionization 



When free electrons are created, they induce a current 
density J ~ q e pv e . This quantity depends on the electron 
charge q e = —1.6 x 10~ 19 C, the electron density p and 
the electro n velocity v e . J is computed from the fluid 



equations ( Esarev et all Il997t ISprangle et all 1996) 



dtp + V ■ (pv e ) =S, 



(31a) 



d t v e + (v e ■ V)u e = — [ E + — X — ) - v e v e - Sve/p- 
m e \ c J 

(31b) 

Here, iS represents external plasma sources and v e is the 
effective electron collision frequency. These equations can 
be combined to yield 



2 

d t J + vJ = ^-S + n, 



(32) 



where 



rl = — JxB-— (V- J)-{J-V)v e (33) 
m e c pq e 

represents ponderomotive forces acting on slowly- varying 
time scales. For linearly-polarized electromagnetic fields 
(E,B) oscillating at the high frequency loq, the driv- 
ing term IT admits envelope components containing gra- 
dients of the field intensity, radiation pressure due to 
electron collisions and changes in the electron den- 
sity. Ponderomotive forces induce low plasma currents, 
which can in turn ge n erate electromagnetic pulses (EMP) 
(ICheng et all 120011 . 12002ft and provide sources of co- 
herent sub-THz radiation (ITzortzakis et ^.1.120021). Nu- 
meric al simulations ( Pehano et a?.l . l2004ISprangle et all 
|2004| ) have shown, however, that for 100-fs pulses reach- 
ing intensities of 10 14 W/cm 2 and free electron densities 
p oj 10 16 cm~ 3 in the atmosphere, the efficiency conver- 
sion to EMP is of the order of 10~ 9 with local intensities 



attaining only 10 kW/cm 2 . In dielectrics, the plasma 
generates EMP intensities remaining about ~ MW/cm 2 
for peak laser intensities of ~ 10 13 W/cm 2 . These pon- 
deromotive terms can thus be ignored, as long as peak 
intensities are below 10 15 W/cm 2 . In this range, plasma 
density perturbations due to Langmuir wave oscillations 
and relativistic increase of the electron mass have also a 
negligible influence. Therefore, the equation for the cur- 
rent density reduces to Eq. (|32|) in which II = 0. At the 
lowest order in v e , the growth of the electron density is 
only governed by the source term S, i.e., 



d tP = S = W(I)( Pnt -p) + —pi - f(p) 



(34) 



that involves photo-ionization processes with rate W(I), 
collisional ionization with cross-section a, and a function 
describing electron recombination or attachment with 
neighboring ions denoted by f(p). Here, p n t and Ui are 
the density of neutral species and the ionization poten- 
tial, respectively, while p <C p nt . Typically, the recombi- 
nation function in gases has a quadratic dependency on p, 
so that f(p) = /? r0 combP 2 with rr . crirnh \cm 3 /s] ~2x 10~ s 
at ele c tron temperatures T ? = 1 eV dMleinek et all 
Il998at ISprangle etall 120021 : ITzortzakis et all l2000bf) . 



Recombination times belong to the nanosecond scale. 
In dielectrics, much shorter recombination times are in- 
volved (r rccom b = 50 — 150 fs) a nd the density l i nearly 
decreases like f(p) = p/T recom h (lAudebert et cd. 
iPefiano et dl 120051: Pfzortzakis et aLl . l2001dl ). 

Besides, the electron collisional rate depends on 
the electron energy distribution function and temper- 
ature versus the ionization potential Ui. Assuming a 
Maxwellian distribution function for the electron velocity, 
this rate linearly varies like a \£\ 2 /Ui, as long as the elec- 
tron thermal energy is small compared wi th XJj. Here, a is 
the inverse bremsstrahlung cross-section (lLotzl . ll967al lbT). 
If we omit Ohmic heating, solving for the current density 
leads by Fourier transformation to 



J = 



-(v e + iuj)(pE). 



m e (v 2 + uj 2 ) 

The current density term in Eq. ((Ta|) transforms as 



(35) 



Pod t J 



+ 



C C 2 p c (l + V 2 /(jJ 2 ) 

after introducing the critical plasma density 
Wn 2 m e eo 1.11 x 10 21 

Pc = - 2 - 7— Cm 



](pE), (36) 



(37) 



at which the laser wave number vanishes. The cross- 
section 



1e 



m e e n G cv e (l + w 2 /f 2 ) 



(38) 



then provides the frequency-dependent collisional 
rate. This expression was earlier derived in the 
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limit of z er o ela stic collisions with ions by Kennedy 
( Kennedy!. 1995), Yablonovitch a nd Bloembergen 
( Yablonovitch and Bloembergenl . Il972f ) and by Feit and 
Fleck (jFeit and Fleckl . Il974l) Often linked to what is 
called the "Drude" model, it determines energy losses 
through plasma (cascade) ionization. 

In Eq. (f3"4"| . W(I) denotes the rate for photo- 
ionization. It is evaluated from perturbative theories 
valid as long as the electric field E is weaker than 
the atom field strength E at . This rate has been red- 
erived in Appendix [A] following Keldysh's and Perelo- 
mov, Popov and Terent'ev (PPT)'s theories apply- 
ing to atoms or diele c trics (c rystals') dKeldvshl . 
Perelomov and PopovL Il967t iPerelomov et all 



1965 



1966, 



1967h . To describe the ionization of complex atoms, PPT 
formula usually includes the so-called ADK coefficients 
[for Ammosov, Delone and Krainov ( Ammosov et all 
1986)], originally established in the limit of high inten- 
sities. Optical field ionization theories stress two major 
limits bounded by the " adiabaticity" Keldysh parameter, 



7 = w - 



^/2 m e U l 

\e7 



(39) 



namely, the limit for multiphoton ionization (MPI, 7 ^S> 
1) concerned with rather low intensities and the tunnel 
limit (7^1) concerned with high intensities, from which 
the Coulomb barrier becomes low enough to let the elec- 
tron tunnel out. Here, E p denotes the peak optical am- 
plitude (E p = \J2c\T). For simplicity, we ignore the phe- 
nomenon of above-thresho ld ionization ( Agostini et all 
I1979L ICorkum et all Il989h . through which the electron 
embarks more kinetic energy than KHujq — Ui. For laser 
intensities / = \£ | 2 < 10 13 W/cm 2 , MPI characterized 
by the limit 



7 » 1 => W{I) -> WWi = tr K I 



K 



(40) 



dominates, where K = mod(Ui/hu>o) + 1 is the num- 
ber of photons necessary to liberate one electron. For 
higher intensities, tunnel ionization starts to contribute, 
i.e., electrons tunnel out within one optical cycle. De- 
spite the complexity of ionization formulas, all of them 
exhibit common dependencies on the laser field strength 



exnibit com mon dependenci es on the laser neld strength 
([Perry et all ll988UlReissl . ll980l ). By producing charged 
ions of noble gases as a function of the laser inten- 
sity, quantitative measurements confirmed t he validity of 
Keldy sh theory for singly-charged ions only (jPerrv et all 
Il988al ). This theory basically applies to hydrogen atoms 
and discards the Coulomb field of the residual ion on the 
outgoing electr on. The PPT/ ADK or Krainov's models 
dKrainovl Il997h . instead, better reproduce experimen- 
tal ionization rates of atoms in, e.g., the t unnel regime 
and may even suit for higher- c harge d states (|Augst et ali , 
Il99lfc ICornaggia and Hermgl . l2000h . On this basis, ion- 
ization of diatomic molecules can be predicted through 
semi-empirical theories using the PPT rate in which the 
electron tunnels through a barrier with effective poten- 
tial Z g/r. Here, Z e g is determined by fitting the theo- 




Beam intensity [ W/cm ] 



Beam intensity [ W/cm ] 



FIG. 2 (a) Ionization rate of O2 molecules vs. laser in- 
tensity obtained from the PPT theory (solid curve), ADK 
molecular model (dashed curve), the fitting curve from PPT 
with Z cB = 0.53 (dotted curve) l|Talebpour et all I1999T ) 
and the MPI-like form ulation (dash-dotted curve) used in 
(|Couairon et ali, 120021 ) at 800 nm. (b) Keldysh ionization 
rates for crystals applied to fused silica (solid curve, Ui = 7.8 
eV) and water (dashed curve, Ui — 7 eV) at the same laser 
wavelength. 



rctical slope of the PPT formula onto measured ion sig- 
nals. For ion signals collected from the interaction of O2 
and N2 molecules with a 800 nm laser pulse, Z e g is in- 
ferred from PPT applied to a corresponding atom (Xe 
atom) with ionization potential close to that of dioxy- 
gen molecules (C/f 2 ~ 12.1 eV, Z cS = 0.53) and to Ar 
atoms with Ui close to the ene rgy gap of N 2 m o lecule s 



15.6 eV, Z qS = 0.9) (jTalebpour et all Il999h . 



The PPT/ ADK formula can alternatively be extended 
to O2 molec ules by selecting t he molecular coefficients of 
Tong etal. (|Tong et a/.l . l2002t ) within the so-called "AD K 
Molecular" ionization model (jNuter and Berget [2006) . 

Figure[2]illustrates ionization rates for some of the pre- 
vious theories applied to O2 molecules [Fig. 0a)] or to 
fused silica [Fig. 0b)]. PPT and ADK molecular ioniza- 
tion rates plotted in solid and dashed curves differ by only 
one decade. The dotted curve indicates the experimental 
fit from O2 signals with the PPT rate using Z e g = 0.53. 
A good agreement is achieved between this fit and the 
ADK molecular ionization curve. Note that as long as 
the beam saturates below 10 13 W/cm 2 , the MPI limit 
alone can be retained. Beyond this intensity, the photo- 
ionization rate departs from the law <jkI K and must in- 
clude tunnel contributions. To get qualitative behaviors 
only, the MPI rate may, however, be used at higher in- 
tensities, while keeping the interaction physics valid. 

Because free carriers are first generated by photo- 
ionization, we must take the corresponding losses into 
account. Via energy conservation law, self-consistent ex- 
pressions for these losses can be established. The tempo- 
ral evolution of the energy density w is determined by a 
local version of the Poynting theorem, i.e., 



jw{r,t) = J(r,t)-E{r,t), 



(41) 



from which we can compute the energy lost by the 
pulse w hen it extracts electrons through single ionizatio n 
process (|Kandidov et all l2003t |R ae and Burnett! . [19921. 
The amount of energy per time and volume units is then 
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given by J ■ E — U t d t pm where d t ppi = W(I)(p nt - p). 
Using complex- valued fields, the current associated with 
photo-ionization losses is merely found to be 



ko TT W(I) 

Ui — (pnt 



2uj p Q 



p)(£ + £*), (42) 



where fast oscillati ons are cance l ed in the denomina- 
tor. Geissler et al. ( Geissler et all fl999t ) derived similar 
losses by introducing the overall polarization vector for 
free electrons V — q e pj>\x, so that dtP — q e pp\x + J. 
Free electrons arise with zero velocity at the position 
x(t) ~ UiE / (2\q e \I) , yielding the same loss current. 

As a final result, the propagation equation within the 
UPPE description (fT5|) reads in Fourier space as 



^£ = 

dz 



2k(uj 



■V± +ik{uj) 



£ 



IPOLU 



ik 2 Q{uo) 



2k{uo) yfCy 
pi\ 0(w) 



e(w)Pi 



NL 



2c(w )fcH(l + ^) \Pc 



' e(^o) 
e( W ) 



where 



-H(te-P) + -^P 



'dt. (44) 



-Pnl(^*, [Eq. ([30])] and the expression containing the 
electron density p(r, t) [Eq. (|3"4"))] must be transformed 
to Fourier space. Treating the complex field £ only for 
positive frequencies is sufficient because of the symmetry 
£*(u>) = £(—uj)*, which imposes to select the parts of the 
nonlinear terms belonging to the frequency range u) > 0. 

For practical use, the collision cross-section a is stated 
for a central frequency ojq. We furthermore assume 
ul/wl < 1 and yj e(u} ) / e{u) w 1. The link to the NEE 
model is then straightforward, whenever the dispersion 
relation supports a Taylor expansion around ujq. By re- 
taining only waveforms beating at loq, the nonlinear en- 
velope equation for the forward component U is directly 
inferred from (|23[) as 



U I i n LUn 

■5-17 = 7n-T V± U + lVU + i—ri2T 
oz 2k Q c 



i^n 4 T\U\ A U- 
c 



{l~x K )\U\ 2 + x K / h{t-t')\U{t')Y dt' 



u 



l—^—T- L pU- ^pU 
2n 2 p c p 2 H 



Pm?k{\U\) 



u. 



j t P = W(I)( Pnt -p) + ^ P \U\ 2 - f(p) 



(45a) 
(45b) 



where t stands for the retarded time variable t — z/v g . half- maximum (FWHM) is At = v21n 2t p . The initial 



The function (3 M pa{\U\) = (p nt - p)[W(/)/|£/| 2 ac- 
counts for losses caused by photo-ionization. In the 
MPI limit (|4T)]) . this dissipative function takes the form 
Pmpa(\U\) -> I3^\U\ 2K ~ 2 where = Kfaj VKPnt 

is the coefficient for multiphoton absorption (MPA). The 
first term of the operator V corresponds to group- velocity 
dispersion (GVD) with coefficient k" — d 2 k/duj 2 \ LJ=UHr 
Equations (|4"5)) describe wave diffraction, Kerr focusing 
response, plasma generation, chromatic dispersion with 
a self-consistent action of deviations from the classical 
slowly-varying envelope approximation through space- 
time focusing and self-steepening operators [(T _1 Vj_£) 
and (T\£\ 2 £), respectively]. They are usually integrated 
numerically by using, e.g., initially singly- humped pulses 
taken with a super-Gaussian beam shape 



r 2N r 2 

U(x,y, z = 0,t) = U e w ° 



-id 



(46) 



which may be focused through a lens of focal length / and 
be temporally chirped if C ^ 0. Here, r = yj x 2 



y 



For Gaussian beams (N = 1), Uq = y/2P in /irwQ in- 
volves the input power Pj n , wq is the beam waist and 
t p the 1/e 2 pulse half- width, such that its full- width- at- 



level (t = — oo) of plasma density is zero. In linear prop- 
agation, such Gaussian pulses diffract over the distance 



z f = (f 2 /z Q )/(l + f 2 /z 2 ), 



(47) 



where zq = tzu^Wq/Xq is the diffraction range of the col- 
limated beam (/ = +oo). In nonlinear propagation, the 
interplay between all competitors in Eqs. ()45|) severely 
alters beam diffraction and can maintain the pulse in 
self-guided state over longer distances. As an exam- 
ple, Fig. [3] shows the dynamics of the maximum fluence 
(F = /_ °° \£ \ 2 dt) for an unchirped, collimated Gaussian 
pulse with 0.9 and 0.45 /iJ energies, input waist wq = 71 
jLtm, half- width duration t p = 42.5 fs propagating in a 
3-cm long silica sample. This beam is numerically sim- 
ulated from the UPPE model [Fig. E^a)] and the NEE 
including higher-order dispersion up to n — 5 [Fig.[3jb)]. 
The ionization rate is delivered by the Keldysh's formula 
for crystals [Fig. |2fb)]. Silica promotes strong chromatic 
dispersion, which is approximated by the NEE model, 
but not by UPPE that describes the full linear disper- 
sion relation. Despite this difference, we can observe that 
the beam fluence remains identical according to both de- 
scriptions (|43p and (|45|) . This property is true in several 
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FIG. 3 Maximum fluence for pulses propagating in silica glass 
with 0.45 \i2 (solid curve) and 0.9 /j.3 (dashed line) energies 
simulated from (a) the UPPE model, (b) the NEE model with 
dispersion limited to fifth order. 



media (dielectrics, liquids and air as well), whenever Eq. 
(I45al) includes dispersive co ntributions 



with n > 3 



Kolesik et al. I. l2003bl[2003) . 

Earlier studies on femtosecond pulse propagation 
in the atmosphere started with simpler models than 
Eqs. (l43l) or (l45l) by. e.g., ignoring time dispersion 



(Kandidov et al.. 1994). Their derivation was im- 



proved later dFeng et all 1 19951: iMleinek et all , Il998allbl : 
ISprangle et q/T " 19961 ). before being finalized into NEE 
and UPPE formulations. 



III. OPTICAL ULTRASHORT FILAMENTS: A FEW 
TOOLS FOR THEIR ANALYTICAL DESCRIPTION 

A. Kerr focusing and wave collapse 



For beams with no temporal dispersion, Eq. (|45aP with 
no other nonlinear response than the instantaneous cubic 
Kerr term (xk — 0), i.e., 

^-U = ^-VlU + *^n 2 \UfU = Q, (48) 

OZ ZKq c 

describ es the self-focusing of optical wav e-packets with 
r&2 > (|Askar'vanlll962tlMarburgerUl975l ). Self- focusing 
is a nonlinear phenomenon common to several branches 
of physics. In optics, it intervenes through the refractive 
optical index N opt — nQ+n 2 I that increases with the field 
intensity and forces light rays to refract onto the axis. 
This causes a compression of the beam in the diffraction 
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FIG. 4 (a) Principle of wave self-focusing. Insets (b) and 
(c) detail some intensity profiles of the solution to Eq. (|48fl . 
applied to silica (ri2 = 3.2 x 10~ 16 cm 2 /W). The initial con- 
dition is the Gaussian pulse (|46[1 where / = +oo, Wo = 130 
/an, Pin/Per = 3, t p = 85 fs and C = 0. 



plane, which leads to "wave collapse" when the Kerr non- 
linearity is not saturated. A necessary condition for the 
collapse is that the input power P in = f \£\ 2 dr exceeds 
some critical value, 



3.72X1 

87rno^2 



(49) 



computed on the Townes mode (see below). The beam 
waist decreases more and more, as the field amplitude \ U\ 
diverges. Figure 2] illustrates the self-focusing principle 
taking place in purely Kerr media. The two insets detail 
the collapsing evolution of a Gaussian wave-packet with 
initial condition (|46|) at z — and near the collapse dis- 
tance z c . This distance, also termed as "nonlinear focus" , 
locates the point at which the beam amplitude diverges 
along the optical path. It is give n by the well-kno wn 
semi-empirical Marburger formula ( Marburgerl Il975| ) 



0.367^ 



'fg-- 0.852 j -0.0219+^ 



(50) 



Below we briefly review properties of wave collapse and 
how self-focusing is affected by temporal dispersion and 
plasma generation. In order to argue on the relevant 
mechanisms participating in this process, we discard non- 
resonant terms causing harmonic generation and con- 
sider an instantaneous cubic response (xk = 0). We 
apply the approximation T -1 ~ 1 — (i/u>o)dt supposing 
uJotp 3> 1) neglect recombination and treat the plasma 
coupling term driven by MPI only and subject to the 
limit T^ 1 — > 1. Equations (|4"5| are then rescaled in di- 
mensionless form by using the substitutions r — > war, 
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c 2 ip and p — > (n 2 l p c /2z k )p, 



t — > t p t, z —* Azqz, U 
where C2 = A 2 , 
tended nonlinear Schrodinger (NLS) equation 



ir 2 non2WQ. They result into the ex- 



T = -5dfy -p^- e\ip\ A ip + iv\il>\ 



2K-2 



tpLOQ 



d t p = T\^\ 



2 A' 



(51a) 



(51b) 



where ^(ip) includes "perturbations" of a critical col- 
lapse, such as GVD, MPI, nonlinear saturation, MPA 
and pulse steepening. GVD and MPA have the normal- 
ized coefficients 8 = 2z k"/t 2 and v = 2z (3 ( - K ^ cf~ x , 
respectively. Quintic saturation is taken into account 
through e = 114,02/ n%. The rescaled MPI coefficient reads 
as T = (2zoko/riQp c )aKPnttpC2 ■ For Gaussian beams, 
the incident amplitude ([46]) reduces to 



|#s = 0)| 



167m «2-Pir 



(52) 



1. Principles of wave self-focusing 

With a purely cubic nonlinearity, solutions to the 
Cauchy probl em (151aD can blow-up (or collapse) a t fi- 
nite distance (|Kellevl . ll965HSulem and Suleml . Tl999h ■ To 
describe this singular process, let us consider the NLS 
equation ()51a|l with T = dft/j (5 = —1), so that the 
Laplacian — > V 2 = d 2 + d 2 + d 2 + ... formally ac- 
counts for the dispersion of a wave-packet along D or- 
thogonal spatial axes [r = (x, y, t, ...)]. The wavefunc- 
tion ip evolves from the spatially-localized initial datum 
ip(r, 0) = ipo(r), assumed to belong to the Hilbert space 

H 1 with finite norm = (HV'lll + II V?/>||!) 1/2 , where 

\\f\\p = (/ \ f\ p dr) 1 l p . Two invariants are associated with 
ip, namely, the L 2 norm (power) P and Hamiltonian H: 



P 



2 

2 • 



H 



|VV||| 



V'l 



(53) 



The following "virial" equality can be established 
dGlassevl fl977t IVlasov et fli.l . ll974h 



Pd 2 z {r 2 )=A{2H+{l-D/2)U\\t}, 



(54) 



where (r 2 ) = J r 2 \ip\ 2 dr/ P denotes the mean-squared ra- 
dius of the solution -0. By a double integration in z, 
Eq. ([54l shows that, whenever D > 2, there exist ini- 
tial conditions for which (r 2 ) vanishes at finite distance, 
which is the sufficient signature of a wave collapse. For 
finite norms P, the inequality P < (4/D) 2 (r 2 ) x ||VV>||| 
thus implies that the gradient norm diverges in collapse 
regimes. As H is finite, the collapse dynamics makes 
the L 4 norm \\ipWf blow up in turn and max r |-0| di- 
verges accordingly, by virtue of the mean-valu e theo- 



rem J {-ipl^dr < max r \tP\ 2 x P (|Kuznetsovl . ll996D . This 



leads to a finite-distance blow-up, at which the solution 



tp stops to exist in H 1 ( Rasmussen and Rypdall . Il986h . 
This mathematical singularity reflects the ultimate issue 
of the nonlinear self-focusing in the absence of saturation 
of the Kerr response and under the paraxiality assump- 
tion. 

While H < arises from Eq. (fB"4"]) as a sufficient con- 
dition for collapse, sharper requirements can be derived 
by means of the Sobolev inequality 



< C\\V1>\\% x ||V>| 



|4-£> 
\2 



(55) 



In the critical case D = 2, this inequality can be used to 
bound H from below, so that the gradient norm blows up 
only if P fulfills the constraint P > P c . The best constant 
in Eq. ([53)) is exactly Cb es t = 2/P c and it involves the 
quantity P c = J R 2 dr = 11.68, where R is the radially- 
symmetric soliton solution, called the "Townes mode" , of 
-R+ r^drrdrR+R 3 = (|Chiao et a/.l . ll964HWeinsteinl . 
Il983| ). P c justifies the existence of a critical power for 
the 2D self-focusing of optical beams in nonlinear Kerr 
media. In the supercritical case D = 3, a criterion for 
collapse, sharper than H < 0, can be established from 
a combination of Eqs. (O and (J55J) as H < P 2 /P for 
gradi ent norms initially above 3P 2 /P ( Kuznetsov et all , 
Il995f ). Here, P c again corresponds to the mass of the 3D 
soliton satisfying —R + r~ 2 d r r 2 d r R + R 3 = 0. 

Once collapse is triggered, the solutio n focuses self- 
simil arly near the singularity point (jRypdal et all , 

1985) 



(56) 



where £ = r/L(z), £(z) = J* du/L 2 (u) and the parame- 
ter A is positive for making the new wavefunction <fi local- 
ized. The function Liz) represents the scale length that 
vanishes as collapse develops, and <j) converges to an ex- 
actly self-similar form <j>{£) fulfilling 8q4> — > 0. For radial 
solutions Eq. (|51a|) transforms into 



»D-1 



+ |<^> + /?[£ 2 -£f# = 0, 



where £§. = /3 1 [A — iLL z {^ — 1)] is viewed as a complex 
turning point, with = —jL 3 L zz . As L{z) — ► 0, <j> 
can be treate d by means of quasi-self-similar techniques 
(|Bergel H998). The solution <f> is split into a nonlinear 
core, (j)c, extending in the range £ <ti £t, and a linear 
tail, 4>t ~ e~ A?T /' 3 /£ 1+^A / ,3 , defined in the complementary 
spatial domain £ 3> £t where the nonlinearity vanishes. 
The length L(z) is then identified from the continuity 
equation describing the mass exchanges between the core 
and tail parts of ip. The dynamics of self-similar collapses 
vary with the space dimension number as follows. 



• For D 



2, L(z) has a twice-logarithmic correction: 
Lqv^- z/ \ /ln In [1 / ( z c — z ) ([Fraimanl . 



zAAnlnJ l TTi^ 
119851: lLandman etail 1 1988: iMalkinl . Il99dh . As 

z — > z c , the exponential contribution of the tail 
decreases to zero, while the core converges to the 
Townes mode R. The power P relaxes to the 
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critical value P c and stays mostly located around 
the center. Self-similar relaxation of self-focusing 
beams to the Townes mode ha s been reported ex- 
perimentally (jMoll et al.l . l2003l ). 



• For D = 3, /3 attains a fixed point /3q ^ 0, leading 
to the scaling law L(z) ~ Lq^Jz c — z. The power 
is no longer preserved self-similarly in space, 
since P = L(z) J |</>| 2 <i£. This integral behaves 
as P ~ Pcorc(^) + PtB&(z), where F C orc(^) ~ L(z) 
vanishes, while Ptaii(z) contains almost all the 
initial mass as L(z) — » 0. A 3D collapse is thus 
accompanied by an expulsion of mass towards 
the large distances where it keeps a station- 



ary de nsity r 
199 It lleMesurier . 



co nst (|Kosmatov et al. 



2000; 



Zakharov and Kuznetsov 



■Vlasov et all 11988 : 
19861b 



2. Variational approaches 

Because the cubic NLS equation is not integrable at 
high dimension numbers D > 2, approximation methods 
describing the fate of singular solutions may be employed. 
Among those, "variational approaches" consist in build- 
ing a set of dynamical equations governing the size, am- 
plitude and phase of the beam. These beam parameters 
are evaluated from functional integrals (Action or La- 
grangian integrals) computed from a given trial function, 
which the pulse is supposed to keep along the optical 
path. Such methods give global, qualitative behaviors 
within a pedestrian way. Their principal drawback is 
that, as they need to fulfill the main conservation laws 
(e.g., preservation of the L 2 norm), they capture the en- 
tire initial mass of the solution and do not let the latter 
evacuate radiation to the boundaries or relax to an exact 
ground state. This discrepancy may partly be cured by 
accounting for corrective damping of the nonlinear core 
(|Arevalo and Beckeii l2005t iKath and Smvthl . Il995| ). but 
improvements become rapidly tricky and limited. Vari- 
ational principles applied to NLS solutions were exten- 

■ 1 1 1 1 

sivcly studied (Anderson, 1983; Anderson and Bonnedal, 

~9t lAnderson et all Il979t iBondeson et all I197S : 



1979 



Desaix et all 1991 ) and f urther exploited in differ- 



ent contex t s (ICerullo et all. 1996t iKarlsson et all 1992 ; 



Manassahl. Il992t iManassah et all. Il988l: iPietsch et al. 



199l tlSilberbereLll990tlVidal and Johnstonl . ll997h . Thev 



were also widely worked out by several authors in 
the present scope (ICouaironl. l2003bt Esarev et all Il997t 



iLehner and Aubvl . l2000t ISprangle et all Il996f ) . 



For example, considering the shape (pffj)) with N 
a Gaussian trial function in the form 



1, 



L(z)y/T(z) 



(58) 



with <j> = e -«i/ 2 -" 2 /2+*0M ; £ ± = r /L{z), r) = 
t/T(z) can be employed in the virial identities govern- 
ing the transverse and temporal mean squared extents 




FIG. 5 Peak intensities for a 1 mm-waisted, 100-fs collimated 
Gaussian pulse with 22 critical powers. Dashed line: Results 
from Eqs. J59) in the limits T,T _1 -> 1; Solid line: (3 + 
l)-dimensional numerical simulations (initial conditions are 
perturbed by a 10% amplitude random noise). 



(|Berge and Couaironl 120001: ISprangle et all l2002h . Sep- 
arate dynamical systems for the transverse width 
the temporal duration [T(z)] and the intensity factor 
[J'(z)] of the pulse are then analytically established as 



L 



ZZ 

T 



JT(0) 



ZZ 

Jz 



L 3 ATL 3 2(K + 1) 2 L 2K + 1 T K - 1 



S 2 . SJTjO) 
T 3 AL 2 T 2 
2v J K 

R 3/2 (J2 T ) K -^ 



(59a) 
(59b) 

(59c) 



where final dots formally omit contributions from MPA, 
Raman-delayed, quintic and steepening effects. With 
J(0) = 2v / 2~Pin/-Pcr imposed by the initial power, Eqs. 
(|59|) yield qualitative behaviors for the propagation of a 
single femtosecond pulse preserving its cylindrical sym- 
metry. It, of course, cannot depict any temporal or spa- 
tial splitting phenomenon, but offers a global evolution 
of the beam intensity. This evolution fairly supports the 
comp arison with (3 + l)-dimension al numerical simula- 
tions ( Champeaux and Berget [2005h ■ apart from an over- 
estimated nonlinear focus. Figure [5] shows an example of 
peak intensity reached by Gaussian pulses with param- 

100 fs at the 



22, 



w 



1 mm, t 



v 



eters P in /P a 

wavelength Ao = 800 nm in the limit T — > 1 and dis- 
persion limited to GVD. Collapse is arrested by plasma 
generation. The dashed curve corresponds to the solu- 
tion reconstructed from the two-scale variational method. 
The solid one shows the result obtained from numerical 
integration of the propagation equations. 

Besides, alternative procedures, such as perturbative 
methods, may provide qualitative information about 



13 



the changes in the scale L{z,t) induced by contribu- 
tions that can stop the collapse. Here, the radial size 
L(z,t) tends to zero in the diffraction plane as z — > 
z c , but it now depends on time. The principle is to 
perturb a collapsing state that naturally tends to the 
Townes mode i?, by means of L 2 orthogonal perturba- 
tive modes. These will provide an integral relation for 
the ^- dependent function [Eq. (1571)1. directly depending 
on T ( Fibich and Papanicolaoul Il999f ) . At the critical di- 
mension, the zeroth-order collapsing solution, whose core 
converges as 

ip -> i' s = L- x R(flL) exp « + iL z r 2 /4L), (60) 

has its mean radius L(z,t) modified by the perturbation 
\J-\ -C 1. Appropriate functions T are then capable of 
turning the sign of /3 = — jL 3 L zz from positive to nega- 
tive, which predicts the arrest of collapse. This method 
was applied to many "perturbations" of NLS, such as 
normal GVD, steepening and non-paraxial d eviations 
(|Fibichl . 119961 : iFibich and Papanicolaoul Il997l) as wel l 
as fo r NLS with attractive potentia ls ( leMesurier et all 
12004 ISchfodt-Eriksen et all l2001ah . It was later ex- 
tended to a variational system formally avoiding the con- 
straints e, 0, T <C 1 by introducing an amplitude fac- 
tor accounting for dissipativ e losses, ip — > <J 1 3{z, £)ip s 
(jBerge and Couaironl 1200 lbh . Applied to Eqs. (|5T]) for 
pulses having a Gaussian distribution in time, this pro- 
cedure involves the substitution 



J( z = 0,t) = F 
and leads to the dynamical equations 

J 2 



-2f 



/Per, (61) 



M 
4R 



L 6 L ZZ = 1 - J + 4e 



Jx 

J 



2vA 



L 2 2P C 

jK-l 



R 2 £d(:pdt (62a) 



L 2K-2 



B(L~ 2 ) t 

tpUlQ 



(62b) 



where M = j £ 2 # 2 (£)<&f, A= J R 2K d£/P c , B = 3 J + 1, 
and C = J R 2K+2 dt;/P c . Equation (|62b[) originates from 
the power variations, assuming J z /J -C 1. The first 
term refers to GVD, the second one to MPA and the last 
term corresponds to pulse steepening. 

Figure [5] depicts the normalized maximum intensity 
and temporal pulse distortions driven by MPI alone at 
z = 2.4 cm for the solution to Eqs. ([62]) in the absence 



of dissipation {J z = 0). The Townes mode is here ap- 
proached by a Gaussian profile for a beam containing 1.85 
critical powers. The temporal dependencies in L(z,t) al- 
ter the pulse profile. Note that discrepancies occur in 
the maximal amplitude, in the location of the first fo- 
cus and in the temporal shapes. Because they rely on 
trial functions that capture the same power and shape as 
initial, variational procedures cannot accurately describe 
fluctuations in the wave envelope. They should therefore 
be employed for what they are, i.e., as qualitative proce- 
dures yielding first indications on the beam dynamics. 
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FIG. 6 Normalized maximum intensity of Gaussian pulses 
in argon with power of 1.85 P cr integrated from Eqs. (|62p 
(dashed curves) and from a numerical integration of Eqs. (I45[) 
(solid curves). Inset shows temporal profiles at z = 2.4 cm. 




FIG. 7 Temporal profiles of Gaussian pulses (wo ~ 130 /im, 
t p = 85 fs) near the nonlinear focus with p OWGr Jin — 
for the Kerr response alone (dotted curve), for Kerr + normal 
GVD (dashed curve), for Kerr + GVD + T, T -1 (solid curve), 
for Kerr + GVD + T, T -1 + higher-order dispersion (dash- 
dotted curve). In all cases, plasma response is omitted. 



B. Saturation by temporal dispersion 

By mixing GVD and spatial diffraction when T ~ 1, 
the Kerr nonlinearity causes defocusing in time for k" > 
(normal dispersion) and temporal compression for k" < 
(anomalous dispersion), besides wave focusing in the 
transverse direction. The interplay of these processes 
results in the symmetric splittin g of the pulse alon; 
the tim e axis with normal GVD ( Chernev and Petro- 
Il992allbh and to a 3D spat iotemporal collapse witn 
anomalous GVD ( Bergel Il998[ ). In addition, when ultra- 
short pulses develop sharp temporal gradients, the op- 
erator T in front of the Kerr term (self-steepening) in- 
duces a shock dynamics: The field develops a singular 
profile with \i/j t \ — > +oo in the trail (t > 0) of the pulse 
( Anderson and Lisakl . 11983 ) . This dynamic s is reinforced 
by space-time focusing ( Rothenbera . Il992h . Figure [7] de- 
picts temporal profiles of pulses undergoing steepening 
effects in normally dispersive regime (silica glass at 790 
nm), together with a transverse collapse. 

To understand these dynamics, we can consider ultra- 
short pulses as being stacked along the temporal direction 
into different time slices having each their own power, 



e.g., P(t) = Pi n e for Gaussian profiles. Slices located 
at times t < correspond to the front (or leading) pulse; 
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those at t > constitute the back (trailing) pulse. Each 
time slice self- focuses at its respective singularity point, 
z c (t), according to Eq. ([50]) in which the ratio Pi^/P„ 
must be replaced by Pi n e _2t jP„. This scenario is known 



as the "moving-f ocus" model ( Lugovoi and Prokhorovt 
Il974l : IShenl . Il976l ) and yields simple comprehension el- 
ements to figure out the pulse distortions. Here, the cen- 
tral time slice focuses at the shortest distance z c (t = 0). 
Furthermore, z c (t) is positive for t > and negative for 
t < [i c (0) = 0], whereas z c (t) always remains pos- 
itive. Near the collapse distance, we can replace all 



time derivatives by d t 



-z c d z 



-zj) z 



+ zW 

(|Fibich and Papanicolaou! Il997t iLuther et all Il994bfh 
Power fluctuations l|62b[) are then given by 

jdisp 9 R 

^— = 25[-z c /L 2 + i 2 (l/L 2 ) z ] + —z c (L- 2 ) z , (63) 

where L z < in compression regime. With 5 > 0, 
Eq. (|6"3"|) describes defocusing around the central slice: 
Normal GVD transfers power towards non-zero instants, 
symmetrically located with respect to t = 0. Self- 
steepening and space-time focusing moreover produce a 
transfer of power from the leading (z c < 0) to the trail- 
ing portion of the pulse (i c > 0). This asymmetrizes the 
tempor al profile, which was ret rieved by direct exper- 
iments ( Ranka and Gaetal . |l998). Normal GVD alone 
"splits" a focusing pulse into two regular, symmetric 
spikes at powers < 2P C . For higher powers, the peak 
edges develo p shock profi l es and disintegrate into rip- 
plelik e cells ( Fibich et all 120031 : iGermaschewski et all 
1200 ll ). Because of the hyperbolicity of the opera- 
tor — <5<9 2 , one splitting event transforms the op- 
tical field in the (r t) p l ane into an A-shaped wave- 
form (iBerge ; et all 120021: IChristodqulides et all, 12004 : 



Conti et all 120031: iDi Trapani et ali 120031: iLitvak et al. 
2000aUbl: IZharova et a/.l . l2003f) . 



Normal GVD (k" > 0) and plasma formation compete 
at powers moderately above critical to halt the wave col- 
lapse. The stronger the GVD coefficient, the larger the 
power interval in which the collapse is arrested by pulse 
splitting. By solving the cubic NLS equation with nor- 
mal GVD, a boundary <5 C rit(j>), function of the ratio of 
input power over critical, p = P; n / P cr , can be calculated 
in such a way that initial conditions fulfilling 6 > 5 CT n (ft) 
will limit the Kerr self-focusing through GVD splitting 
i nstead of letting the s olution diverge into a singular state 
( Luther et all fl994ar i. Higher-order dispersion together 
with steepening effects (T, T^ 1 ) modify this curve in the 
sense of "delaying" the self - focus i ng threshold t o higher 
powers (|Skupin and Bergel 120061: ISkupin et al\ , l2006afh 
In particular, third-order dispersion tends to delocalize 
the pulse by pushing t he temporal centroid to the back 
(jFibich and Ilanl . I2004D . 

In contrast, for anomalous GVD (fc" < 0) power is 
transferred to center, as seen from Eq. ([63]) with S < 0. 
Ultrashort pulses thus collapse both in space and time. 
A mapping |<5| > <5 C rit (p) can ag ain be constructed on the 
basis of virial-type arguments ( Berge and Skupinl . 120051: 



to 




FIG. 8 Collapse regions in the plane (8, p). The solid curves 
extrapolate results from the propagation equations. The 
dashed curves represent theoretical limits estimated from Eqs. 
(|51[) . above which dispersion inhibits the self- focusing. Re- 
gions I, II and III are explained in the text. 



iKuznetsov et aLl . ll995h . It predicts pulse spreading when 
the dispersion length ~ t^j\k"\ is short enough to prevail 
over diffraction and Kerr nonlinearity. As for normal 
GVD, the theoretical boundaries are modified by higher- 
order dispersion and steepening terms to some extent. 

Figure [5] summarizes the zone of collapse and no- 
collapse in the plane (5, p) for both normally and anoma- 
lously dispersive media computed from Eqs. (|45|) . Open 
circles represent initial conditions that do not collapse. 
Close ones result in a strong divergence of the beam 
intensity. Three zones clearly occur: (I) A dispersion- 
dominated domain leading to pulse spreading; (II) A 
transition zone in which GVD and steepening operators 
inhibit the self-focusing; (III) A Kerr-dominated region, 
in which chromatic dispersion is unable to stop the wave 
blow-up. Here, the field intensity can increase by several 
decades before reaching the ionization threshold. Note 
that for anomalous GVD, the collapse process, which en- 
compasses spatial and temporal compressions, makes the 
intermediate zone II almost meaningless. 



C. Saturation by plasma defocusing 

The action of MPI is to deplete the pulse temporal pro- 
file through the sudden emergence of an ionization front 
near the focus point z c . In the absence of time dispersion 
and nonlinear losses, the plasma equation dtp = T\il)\ 2K 
can be integrated with the ansatz (|6ip as 



p(r,t,z) ~ J — 




L 2K (z,t) 



(64) 



where Erf(x) = (2/y / 7r) J* e~ u2 du denotes the error 
function. The last term of Eq. (|62ap then scales as 
L 2 ^~ K 'd^p where d^p ~ d^R 2K is negative. As L — > 0, 
this term efficiently competes with the Kerr contribution 

-at" 



J - P in e- 2t ~ /P CI . With the step function Erf(V2i«) + l, 
the plasma response arises like a defocusing plateau, as 
the field intensity reaches its maximum near the focus 
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FIG. 9 (a) Maximum intensity, (b) peak electron density and 
(c) mean radius (FWHM of the fluence) of a 150-fs, 0.5 mm- 
waisted Gaussian pulse with 4P cr in air at Ao = 800 nm. 
Dash-dotted curves refer to the case when T,T _1 = 1 and 
chromatic dispersion is limited to GVD (k" = 0.2 fs /cm). 
Solid curves refer to the full model (|45[) . (d) Temporal profiles 
at z = 0.8 m. The dashed curve refers to the distribution at 
2 = 0.4 m (identical for both models), (e) Radial profile vs 
time near z c . (f) Plasma responses associated with (d). 



point z c . All time slices belonging to the interval t > 
are defocused. At negative times, the pulse continues 
to self-focus and feeds plasma defocusing until forming a 
single time slice located nea r t = t* ~ — [In y/P^ / P^] 1 / 2 
at powers close to c r itical ( Berge and Couaironl l2001at 
iHenz and Herrmann! . Il999f ) . Thus, plasma generation is 
accompanied by a sharp duration shortening of the pulse 
(see, e.g., inset of Fig. [6]). 

When GVD and MPA come into play, the leading 
peak becomes unstable and refocusing of the trail oc- 
curs at high enough input powers. With MPA, the lead- 
ing component is partly damped in the intensity ratio 
J z jJ ~ — 2vA{J I L?) K ~ X , as the peak electron plasma 
forms a density plateau [see Eq. (|62b|) ]. Consequently, 
the electron density attains much lower levels and per- 
mits the emergence of a trail. This complex dynamics 
leads to two following signatures for plasma defocusing: 
(i) MPI shortens the pulse duration near the nonlinear 
focus; (ii) Because p scales as R 2K , the spatial zone 
of plasma defocusing takes place inside a narrow region 
(~ 1/ \/K times the beam width) of the intensity distri- 
bution, which creates spatial rings. This "self-guiding" 
mechanism is therefore not static. Certain pulse compo- 



nents are defocused to the benefit of the others. To clear 
it up, Figure [H depicts focusing/defocusing cycles affect- 
ing the peak intensity, electron density, beam radius and 
the distribution in space and time of a 150-fs pulse for 
atmospheric propagation. Note the rings formed just at 
the stage of plasma defocusing [Fig. M,e)] and the occur- 
rence of a sharp trail when steepening effects (T, T -1 ) 
are included [Fig.^d)]. 



D. Saturation by optical nonlinearities 

Additional defocusing quintic nonlinearities pro- 
duce staMe_s^htons_Jii conti nuous- wave (cw) me- 
dia (jTrillo and Torruellasl 1200 ll) . This property fol- 
from a balance between diffraction, Kerr fo- 

When a quintic 
-e\ip\ 4 ip), the beam 



lows 

cusing and nonlinear saturation, 
contribution is added {T 
starts to self-focus and evacuates part of power into 



of Eqs. (I5T1) ( 


Gatz and Herrmannl 1997; Kolokolov 


1976: Malkin. 1993: 


Vakhitov and Kolokolov, 1975 


Vidal and Johnston. 1996). This property is partly re- 



fleeted by Eqs. (f6"2"|) emphasizing an arrest of collapse like 
[3 = —jL 3 L zz ~ 1/L 2 , which induces regular oscillations 
of the size L(z) beyond the focus point. 

Reframed in the present context, several papers ad- 
dressed the question of limiting self-focusing by nonlin- 
ear optical sat urations, instead of relying on plasma gen- 
eration alo ne ( Gaeta and Wisd . l200lL iKoprinkov et all 
l2000l |2001( 1. In fact, it appears that the weight of x (5) 
susceptibility depends on the actual pulse length, on the 
ratio of instantaneous over delayed Kerr responses and 
on the ionization sources, that in fluence the ma x imum 
intensity attained by the beam (lAkozbek et all l200ll ; 
ICouaironL |2003al IVincotte and Bergel |2004T i. This prop- 
erty can be seen from the zeroes of the nonlinear refrac- 
tion index 



An = [(1 - x K ) + x K G(t)M 2 - e|^| 4 - p, 



(65) 



in which the Raman-delayed response has been intro- 
duced by changing the Kerr term \ifj\ 2 of Eqs. (fSTj) into 
[(1 - x K ) + x K G{t)}\ii)\ 2 . Here, the function G(t) fol- 
lows from computing the integral J h(t — t')\ip(t')\ 2 dt' ~ 
G(t)\ip\ 2 over the initial pulse profile (f52"|) . High-order 
saturation prevails over ionization and induces a soliton- 
like dynamics at high enough values of before the 
occurrence of an electron plasma. Relevance of higher- 
order optical nonlinearities above all depends on the ion- 
ization rate: For example in air, if the laser intensity sat- 
urates above 10 14 W/cm 2 by plasma generation alone, 
then even weak values of can soften this peak inten- 
sity. At lower intensities, quintic saturation has a more 
limited role. For instance, Fig. [TU] displays the peak in- 
tensities, beam radius (FWHM of the fluence distribu- 
tion), maximum electron density and energy losses of a 
70-fs, 0.5 mm-waisted unchirped Gaussian pulse prop- 
agating in parallel geometry in air. Limitation of the 
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FIG. 10 (a) Peak intensities, (b) pulse radii, (c) electron den- 
sity and (d) energy losses for Gaussian pulses with u>o = 0.5 
mm, t p = 70 fs, T = T _1 = 1 and different ionization models: 
PPT (solid curves), ADK molecular (dashed curves), ADK 
molecular + susceptibility (dotted curves) and the MPI- 
like approximation of Fig. [2] (dash-dotted curves). 



pulse growth depends on the ionization model employed. 
One of them involves the ADK molecular rate without 
(dashed curve) and with (dotted curve) a weak quintic 
saturation (714 = 2.5 x 10~ 33 cm 4 /W 2 ). The measured 
intensity peak of 5 x 10 13 W/cm 2 is refound in the latter 
configuration, susceptibility lowers the peak inten- 
sity and density, which stays close to 10 16 cm' 3 . Energy 
losses are weaker, which increases the self-guiding range. 

Other high-order effects may soften the Kerr focusing 
of nonlinear waves. For instance, the Raman response 
function (f2"T)) becomes all the weaker as the pulse is 
short, which affects the critica l power for self-focu sing like 
-Per -> P CI /{l-XK+x K G(t)} (|Ting et aUl2005aj) . Unlike 
long pulses, for which G(t) — ► 1, the non-instantaneous 
nonlinearity reduces the effective value of the Kerr index 
and results in delay ing the nonlinear focus from which 
a filament emerges (lAkozbek et all , l200lt IChiron et all , 
Il999t iLiu and Chml 12005 ). In connection, modulational 
instabilities have t heir growth ra te decreasing as the ratio 
xk is augmented (jWyllerl . 120011 ). 

Besides, the cubic polarization ~ E 3 generates third- 
order harmonics as 



E 



IT c^k(uj)z—iujt j_ f j ik(3u))z—3iujt , 



(66) 



at w = wo (lAerawai l200ll) . Two coupled equa- 
tions can then be derived under the basic assumptions 
of slowly-varying envelopes and narrow-spectral band- 
widths Autj/ujj < 1 (j = u, 3a;). The coupling is insured 
by cross-phase modulation (XPM) and four-wave mixing 
(FWM) induced by the cubic nonlinearity. By substitut- 
ing U 3uJ — > U^e lAkz , where Afc = 3k(uj) — k(3u) is the 
linear wavevector mismatch parameter, the propagation 
equations read like Eq. (|51|) where the function T for the 



- + 2|^| 2 ^ + ^ + (^-d t - 

3 Av (4zoj 



(67a) 
(67b) 



that includes the group- velocity mismatch responsible for 
temporal walk-off, Av — [w ff (3w) _1 — t> s (;j) _:l ] _1 . In 
self-focusing regimes, the third-harmonic intensity usu- 
ally con tributes by a little per centage to the overall beam 
fluence l Akozbek et all , [20021 ). Despite the smallness of 
the TH field, this component may act as a saturable non- 
linearity for the pump component when the term contain- 
ing |Afe| has an order of magnitude comparable with the 
FWM and XPM terms, tp 3uj ~ ^/3(4z \Ak\ + 2|</>„| 2 ), 
and all other contributions driving the TH component are 
neglected. Once inserted into t he pump wave equat ion, 
this so-called "cascading" limit ()Burvak et all , [l999h in- 
troduces a saturable quintic nonlinearity that lowers the 
peak inten sity and enhances the propagation over longer 
distances ( Berge et all I2005T ). 

So far, we have been dealing with linearly-polarized 
waves. However, by writing the electric field as £ = 
(£ x ,£ y ), each polarization component undergoes optical 
nonlinearities as 

pnl ~ \n^3 + + -tf^z. (68) 

with j, k = x,y (j k). A linearly-polarized state is 
described by the single-component NLS equation with, 
e.g., ip y — 0. A circularly-polarized state corresponds to 
the configuration ip y = ±iip x . In between, elliptically- 
polarized states facilitate energy transfers through the 
phase-dependent terms t 4)?jil>ZI'3 from one or t hogona l 



component to the other ( Schiodt-Eriksen et all . l2001bh . 
These energy transfers in c rease the power threshold for 
self-focusing ([Berge et all , I2003). For circular polariza- 
tion, the collapse threshold for the total power P = 
/(IV^I 2 + li'yl 2 )^ is increased by a factor 3/2 when 
passing from linearly- to circularly-polarized geometries. 
In addition, circular polarizations on femtosecond fila- 
ments self-channeling in the atmosphere produce sta- 
bler patterns than th ose issued from linear polarization 
(jKolesik et all , l200ll ). They moreover weaken the MPI 
efficiency and decrease ion signals by a couple of decades 
in the fluores cence of N 2 mole cules collected from in- 
frared pulses ([Petit et aLl . l2000h . 

Compared with plasma formation and chromatic dis- 
persion, quintic saturation and polarization effects have 
a limited impact on the nonlinear dynamics. Therefore, 
we shall henceforth focus on the first two players mainly. 



E. Self-phase modulation and supercontinuum generation 

The Kerr effect creates spectral broadening through 
self-phase modulation (SPM). The solving for id z 4> = 
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— \i^\ 2 ip indeed yields the exact solution tp = ip e l ^°\ z 
[t/j = tp(z = 0)], which describes a self-induced phase 
modulation experienced by the optical field during its 
propagation. This intensity-dependent phase shift is re- 
sponsible for spect ral broad e ning by virtue of t he relation 
Aui = -<9 t arg(V0 (jAgrawall 12001 Ishenl . 119841 ). Because 
the frequency spectrum is expanded by the nonlinearity, 
SPM leads to supercontinuum generation and white- light 
emission, as the wave intensity strongly increases through 
the self-focusing process. With plasma generation, this 
phenomenon can be described by employing the following 
approximations. First, we assume that the MPI response 
behaves like a static density plateau [see, e.g., Fig.[9jf)], 
so that the ratio p/\tp\ 2 is either zero or close to unity 
with d z (p/\ip\ 2 ) ~ 0. Second, omitting GVD and MPA 
for simplicity we retain self-steepening to the detriment of 
space-time focusing, i.e., we neglect diffraction in Kerr- 
dominated regimes. Phase and amplitude of solutions 
V> = Ae^ to Eq. (j5Ta|) satisfy <p e + j^dtip = A 2 - p 
and A z + -r^—A 2 d t A — 0, respectively. If we model the 

pulse initially located at t = to by A 2 — Aq/cos)li(t), 
r = (t — io)/ T 0: an d suppose that the amplitude does not 
change too much, then 



A 2 ~ A 2 /cosh{ T - 3QA 2 /A 2 ), 
ip ~ t p uj T J (l- j^j dr' 
— i p CLiosinh _1 [sinh(T) — Q], 



(69a) 
(69b) 



where Q = zA\jt v ix>oTo is linked to self-steepening. These 
expressions yield the spectral broadening 



AA 
Ao~ 



P_ 

A 2 



Q 2 - 2Qsinh(r) 
cosh 2 (r) 



-1/2 



(70) 



valid under the basic assumption Alo/ujo = — AX/Xq <C 
1 . As long as Q <C 1 and in the absence of MPI, varia- 
tions in wavelength AA/Ao ~ — Qsinh(T)/cosh 2 (r) start- 
ing with to = re present the early symm etric broaden- 
ing through SPM (|Yang and Shenl . 1 19841 ) . When MPI 
forms a defocusing plateau, AX/Xq then becomes larger 
in the region where p = (non-defocused leading pulse) , 
than when p — > A 2 (defocused trail). Consequently, as 
the beam reaches the first focus point z c , the dominant 
part of the pulse is the fr ont edge (p = 0) and MPI 
creates a primary redshift ( Lehner and Aubvl . I2000T ) . In 
the trail, full plasma coupling (p/A 2 — > 1) limits the 
spectral enlargement to the opposite side. However, self- 
steepening induced at increasing Q inhibits the MPI red- 
shifting and instead displaces more the spectrum to the 
blue side. As a result, this creates an asymmetric spectral 
br oadening with a prominent blueshift AA < for r > 
dBerkovsky et cd. , 20051: Champeaux and Berge , 2005 ; 



Gaetal l2000t IWard and Berget 120031: lYang and Shen 
19841) . 

Figure [11] shows spectral modifications around 790 
nm induced by SPM affected by space-time focusing 
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FIG. 11 Spectral broadening induced by SPM and steepening 
effects in fused silica: (a) with purely cubic nonlinearity (inset 
shows the spectrum in frequency for comparison), (b) with 
Kerr + GVD only (dashed curve) and when adding space- 
time focusing and self-steepening (solid curve) for Pj n = 3P C r, 
(c) same with plasma generation that remains inefficient at 
this power level (dashed curve, see Fig. |SJ , but strongly affects 
the spectrum when P ln = 6P cr (solid curve). For comparison, 
II If shows supercontinuum with (solid) and without (dash- 
dotted) T, T -1 in air for the initial data of Fig. [5] 



and self-steepening in the presence (or not) of ioniza- 
tion in fused silica. The properties stressed by Eq. ([70]) 
can be refound. Note the asymmetry in the SPM pro- 
file when the Fourier transform of the on-axis field in- 
tensity is expressed in wavelengths [Fig. UlT a)], since 
Au = — 2ircAX/ X 2 . For comparison, Fig. [TlT d) illus- 
trates supercontinuum reached in air with or without 
pulse steepening at 800 nm. In all cases, we can observe a 
clear amplification of UV/blue wavelengths due to steep- 
ening terms. Because the spectral component rapidly 
decreases outside the central wavelength, the principal 
source of ionization (scaling as I K in MPI regime), how- 
ever, remains given by that operating at Aq. 



F. Modulational instabilities: The route to multiple 
filamentation 

In the context of nonlinear optics, an intense wave 
propagating in a focusing Kerr medium can break 
up into several spots from the inhomogeneities affect- 
ing its initial distribution. This phenomenon was 
first evidenced by self - focus i ng experiments in l iquids 
( Bespalov and Talanovl . [l966l : ICampillo et a/] . Il973h and 
the resulting small-scale structures are usually named 
" filaments" . Because it mostly concerns spatial distor- 
tions, this process of " multifilamentation" can be un- 
derstood from Eqs. (|5ip . where temporal variations are 
discarded. (2 + l)-dimcnsional models freezing tempo- 
ral dependencies were used to simulate the early fila- 
mentation stage of high-power beams in various trans- 
parent media, yielding rather good qu alitative agree- 



ments w ith experimental measurements (iDubietis et all , 
2004a|lc3: Fibi ch et "oil. 12005112004 iKandidov et aLl . l2005l 
19991 : IPorras et all |20o3k In plasma regimes, 2D mod- 
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els can also simulate complex multifila mentation p a tterns 
created from terawa tt laser pulses ([Berge et all |2004 
ISkupin et all l2004bl ). For instance, by assuming that 
MPI mainly counterbalances Kerr focusing at a dominant 
time slice t — t c (z), the field envelope can be decomposed 
as if} = ip(x,y,z) x t c (z)] where the temporal distri- 
bution for the highest-intensity peak is modeled by the 
Gaussian i c (z)] = c _ [ t_tc ( 2 )] / T with constant extent 
r. Plugging the above expression into Eqs. (|5ip . comput- 
ing the expression of p and averaging the result over the 
entire time domain supplies the extended NLS equation 
for the spatial profile tp: 



id z i> = -viv - f{\W)i> 

2 



f(s) 



as 




tTs 



K 



„K-1 



K 



(71a) 
(71b) 



where the coefficient a averages the Raman delayed re- 
sponse. This model does not formally depend on the 
longitudinal location of the time slice t — t c (z). The 
only arbitrariness is the choice of the peak duration r. 
Because MPI shortens pulses to mean duration reaching 
1/10 of their initial values, the value r = 0.1 was found 
to provide the best approximations of fluence patterns 
developed by (3 + l)-dimensional fs pulses. 

Modeling filament formation then requires a perturba- 
tion theory involving a steady-state solution, expressed 
as ip s (r, z) — (p(r)e iXz . Defined in the limit of no dissipa- 
tion [y — » 0) , (p satisfies the differential equation 



- \<t> + vi^- 



(72) 



and A = const. Stability of 4> is investigated from pertur- 
bations v + iw with small real- valued components (v, w) 
acting against this stationary mode. Linearizing Eqs . 
(I7T1) yields the eigenvalue problem ( Kuznetsov et all 
1 19861) 

d z V = L w, i o = A-Vi-/(0 2 ), (73a) 

-d z w = L lV , L, = A - Vi - [/(0 2 ) + 2/'(0 2 )0 2 ], 

(73b) 

where Lo and L\ are self-adjoint operators with f'(4> 2 ) — 
df(u)/du\ u=r f,2. Combining Eqs. (|73|) . we then obtain 
d^v = —LqLiv, from which different kinds of instabilities 
may be investigated. 

Modulational insta bilities: Originally proposed b y Be- 
spalov and Talanov (jBespalov and Talanovi Il966t) . the 
modulational instability (MI) theory involves oscillatory 
perturbations with an exponential growth rate, that split 
the beam envelope approximated by a background uni- 
form solution. Pcrturbative modes are chosen as v, w ~ 
cos (k x x) cos (k y y)e yz and they apply to a plane wave <j> 

which satisfies V\4> — and A = f{4> 2 )- The growth rate 
7 is then given by 



Plane waves are unstable with 7 2 > in the range 
< k± < \ 2A and the maximum growth rate 7 max = A 
is attained for k± = /c max = \f~A. This instability pro- 
motes the beam breakup into arrays of small-scale fil- 
aments periodically distributed in the diffraction plane 
with the transversal spacing A mo d — 27r/fc max and longi- 
tudinal length ~ 7max- The number of filaments is close 
to the ratio Pm/-Pfii, where P$x is the power enclosed in 
one filament. Considering each filament with radial sym- 

A /2 

metry, the evaluation Pgi ~ 2ir J Q m r\4>\ 2 dr ~ 2.65P C 
holds for unsaturated Kerr media [/(s) = s\. 

Because they constitute rough representations of phys- 
ical beams, plane waves can be replaced by the soliton 
modes of the NLS equation (|72[) . The resulting instability 
appears when a soliton <fi is perturbed by oscillatory mod- 
ulations developing along one transverse axis. Perturba- 
tions are, e.g., local in x and they promote the forma- 
tion of bunches periodically distributed over the y axis. 
The operators Lq and L\ in Eqs. (|73p are transformed 
as L = A - d 2 x - f(<f> 2 ) + fc 2 and L x = L - 2/'(^> 2 )c/> 2 . 
Numerical computations are the n often required for solv- 
ing the spectral problem (1731 (lAkhmediev et all fl992 ; 
Rypdal and Rasmussenl[r989HZakharov and Rubenchik . 
19741) . 

Azimuthal instabilities: For broad beams, self-focusing 
takes place as a regular dis tribution of dots loc ated 
upon ring diffraction patterns ( Feit and FleckL[l988h . To 
model this instability, the Laplacian in Eq. (|71a[) must be 



rewritten as V 2 ; 



l d r rd r 



~ 2 dn, where 9 denotes 



the azimuthal angle. Within a first approximation, un- 
stable modes v,w ~ cos(M#)e 7MZ with azimuthal index 
number M break up a spatial ring, which is modeled by 
a uniform background solution <f> lying on a circular path 
with l ength s = f9, and mean radius f. Eqs. (173)) then 
yield (|Atai et all 1 1994 ISoto-Crespo et all 1 19921) 



7m 



2A 



M 2 

-2 



(75) 



and the maximum number of modulations on the ring is 
provided by the integer part of M max = f\f~A. 

In this regards, solitary waveforms having a ring- 
shaped distribution may also exhibit a phase singular- 
ity with an integer number of windings, m (topologi- 
cal charge). Such structures are termed "optical vor- 
tices" . They convey a cons t ant orbital angular momen - 
tum (jKruglov et all Il992t iKruglov and Vlasovl Il985l ) 
and are experimentally designed b y means of phase 
masks and hologr a phic t echniques dDesvatnikov et all 
120051: iPetrov et all Il998t iTikhonenko et all Il996l ). For 
cubic media, they undergo azimuthal MI, which make 
vortex solitons decay into M ma x ~ 2\m\ + 1 fila - 
ments (|Vincotte and Bergel 120051: IVuong et all 120061 ). 
For cubic-quintic nonlincarities, optical vortices may be 
linearly stable fo r \m\ = 1 dBerezhiani et all 12001 



Desyatnikov et al\, 2000|: Quigora-Teixeiro and Michinel 



19971 : ISkarka et all 12003 : rTowers et all l200lh . but 
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fej_(2A-fc 2 ), A ee uf\u)\ u=tj) 2. (74) they decay into M n 



2 1 ml filaments otherwise 
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FIG. 12 Multifilamentation patterns in Kerr focusing regime 
at powers above critical (Pi n < 30P cr ) for (a) a Gaussian beam 
(xo = 0.35, yo = 0, e = 0.1) and (b) an annular ring with flat 
phase (xo = 1.7, = 0, e = 0.01). These two beams undergo 
MI from a local defect with relative amplitude e located at the 
coordinates (xa,yo). (c) Vortex-shaped beams with charge 
\m\ = 2 initially perturbed by a 10% random noise. 



(jFirth and Skrvabinl . Il997t ISkrvabin and Firthl . Il998t ). 
The orbital motion confers robustness to vortex solitons, 
which can propagate p owerful beams beyon d several tens 
of Rayleigh distances (jMichinel et aZlfeOQll ). before they 
break up into filaments. 

Figure [12] shows the multifilamentation of (a) a Gaus- 
sian beam, (b) a flat-phase, ring-shaped waveform and (c) 
an m — 2 vortex breaking up into 2\m\ + 1 filaments in a 
cubic medium. MI is trigg ered either from an initial de- 
fect or from random noise ( Berge et a i. 2003). Note the 
robustness of the Gaussian profile, whose shape defini- 
tively differs from that of a plane wave. 

Modulational instabilities rather affect the early Kerr 
stage in the beam propagation. At later stages, the re- 
sulting filaments become fully nonlinear and they may 
interact mutually. From the interplay between diffrac- 
tion and nonlinearity, two filaments can fuse, when- 
ever their separation distance is bel ow a critical value , 
function of their individual powers (jBerge et all Il997t 
iMcKinstrie and Russelll . Il988ft . With no saturation, each 
filament whose power is above critical creates its own at- 
tractor, at which it freely collapses. By including satu- 
ration, filaments with powers above critical are able to 
coalesce into an intense central lobe an d develop vari- 
ous patterns (jKonno and Suzu 0, 119791) . Without any 
control, multifilamentation happens as an harmful insta- 



bility destroying the homogeneity of the energy distribu- 
tion. It has detrimental consequences on the quality of 
the beam shape. 



IV. UNIVERSAL FEATURES OF FEMTOSECOND 
Fl LAMENTATION 

For an efficient propagation, the beam intensity and 
electron density must both remain below their so-called 
" laser- induced breakdown" (LIB) limits. LIB means 
the breakdown of the beam that results in damage or 
tremend ous energ y abso rption in dielectrics as well as in 
liquids ( Kennedy! . Il995| ). It involves three major pro- 



cesses, namely, (i) ionization, (ii) heating of the "free" 
electrons by the radiation, and (iii) transfer of plasma 
energy to the medium. Once generated, the electron 
plasma absorbs and scatters the remaining pulse en- 
ergy. Assuming free electrons seeded by avalanche ion- 
ization, dtp — {a/Ui)pI, an incident pulse launched in 
tightly focused geometry will be significantly attenuated 
by energy absorption at density levels above a critical 
value, p CI . An estimate of the LIB intens i ty the n fol- 
lows as /lib ~ {Ui / <rt p )\n(p CI ) (iFeng et all KL995). For 
subpicosecond pulses, multiphoton processes mostly pro- 
vide the dominant breakdown mechanism, s uch that the 
estim ate /lib ~ {pa/ Pnt&Ktp) 1 ^ applies ( Feng et all . 
11997ft . Ranges of values for p cr are often close to 
10% of p nt (e.g., p cr — 10 21 c m~ 3 in silica or wa- 
ter) when using femtose cond pulses ( Lenzner et aLl . [l99"8t 
iNoack and VoeelL Il999ft . 

In what follows, characteristics of ultrashort pulse 
propagation are mostly depicted for peak electron densi- 
ties less than the LIB threshold. Density levels maintain 
the pulse in self-channeled state without significant con- 
version of beam energy into the plasma. 



A. Radial self-focusing and temporal splittings 

The self-channeling mechanism relies on the balance 
between Kerr focusing and plasma generation. The 
most illuminating example is given by "femtosecond fil- 
aments" in air, as presented in Fig. [1] Atmospheric fil- 
aments form microstructures exhibiting ~ 100 — 150 /im 
in FWHM diameter and ~ 1 mJ in energy at infrared 
wavelengths (800 nm). They are accompanied by a 
strong SPM-induced spectral broadening and cover sev- 
eral tens of met ers, whenever the input pulse power ex- 
ceeds ~ 10 GW (iBraun et all Il995t iBrodeur etOH Il997t 
iNibbering et a/.l . Il996ft ~ Propagation distances in excess 
of 200 m were even reported for 60 fs, 795 nm pulses 
with 4 mm diameter, traveling thus upon a few Rayleigh 
distances zq. These experiments revealed, after a possi- 
ble short stage of multiple filamentation, the emergence 
of one robust filament coupled to electron dens ities of 
a few times 10 16 cm -3 (jLaFontaine et all [l99 9b) . Fur- 
ther measurements of infrared filament s specified a peak 
intensity of about 5 x 10 13 W/cm 2 ( Kasparian et all 
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l2000a[ ). Similar behaviors involving slightly higher peak 
intensities were als o reported in N 2 m olecular gases at 
different pressures ([Becker et aZJ . l200lh . 

The "moving-focus" model was first revisited to 
yield qualitative explanations of th is phenomenon 
( Brodeur et 

Mi Il997t IChin et all Il999al) . Though pleas- 
ant and rather intuitive, this model, however, failed at de- 
scribing the self-guiding of light th at persists beyond th e 
linear focus of convergent beams ( Lange et all [l998b). 
The appropriate scenario for long distance propagation, 
called " dynamical spatial replenishment" , was specified 
by numerical simulations. During the focusing stage, the 
beam generates a narrow plasma that strongly defocuses 
the trailing part of the pulse and creates one leading 
peak. Once the intensity decreases enough (via, e.g., 
MPA), plasma generation turns off. The back pulse can 
the n refocus, which produ ces a two-spiked temporal pro- 
file (|Mleinek et all Il998allbl ) . In space, plasma defocuses 
only the inner part of the beam inten sity, so that the 
traili ng edge decays into spatial rings ( Kandidov et all 
Il994h . With increasing propagation distances and as the 
front pulse intensity decreases, the spatial rings coalesce 
under Kerr compression and allow refocusing of the trail. 
Although the beam radius may look " quasi-static" , the 
temporal and spatial distributions of the pulse strongly 
fluctuate (see, e.g., Figs. [5] and [TO]). At high enough 
powers, self-channeling is supported by several focus- 
ing/defocusing cycles triggered from a dynamical inter- 
play between multiple, stringent peaks alternating in the 
temporal pulse profile. "Pushed" by steepening effects, 
these peaks are responsible for the formation of optical 
shocks in the medium. This complex dynamics was re- 
ported by sev eral authors, dealing with pulse propagation 
eithe r in air (lAkozbek et all 20011: Berge and Couaironl 



I2OOOI : IChiron et all Il999b ICouairon et all l2002h. in ar- 



gon ce lls (jChampeaux and Bergel 120031: iNurhuda et ali , 
2002a). or even in fused silica jTzortzakis et all l2001dh . 

Because the fundamental scenario for self-channeling 
is that aforementioned, estimates for peak intensities 
(7 max ), electron densities (/? ma x) and filament radius 
(£min) can be deduced from equating diffraction, Kerr 
and ionization responses. This yields the simple relations 



Pn 



tpPntW(I ma ^), 



where 



2p c n n 2 

~ 7r(2fc Q n2/max/«o) j 



n 2 = 712(1 — x K ) + n 2 x_R-max t G(/j), 



(76a) 
(76b) 

(77) 



represents the effective Kerr index when the Raman- 
delayed response is maximum over the initial pulse pro- 
file and W^imax) reduces to W(I ma x) = ^kI^x in MPI 
regime. Some parameter values for different material 
used throughout this review have been indicated in Table 
I. From these values the above relations supply filament 
diameters of the order of 150/im in air at 800 nm and 



TABLE I Parameter values for dioxygen molecules (air) at 
800 nm ([/; = 12.1 eV, p n t = 5.4 x 10 18 cm" 3 ), fused silica 
at A = 790 nm and 1550 nm (U t = 7.8 eV, p n t = 2.1 x 
10 22 cm -3 ), and water at 527 nm (U r = 7 eV, p nt = 3.32 x 
10 22 cm~ 3 ). 
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10 — 20 /im in dense transparent media (dielectrics, wa- 
ter), which agrees with current numerical data. Note 
that for laser wavelengths about 1/im, the requirement 
([9]) is always fulfilled. As inferred from Table I, the criti- 
cal power for self-focusing exceeds the GW level in gases, 
but remains limited to a few MW in dense materials. 

Plasma generation is mostly expected to stop the di- 
vergence of the beam caused by self-focusing, except for 
specific values of the parameters P ln /P cr and 8 ~ k" 
that privilege collapse arrest by dispersion (see Fig. [5J . 
For An = 800 n m, nonlinear fluorescence techniques 
( Liu et all . l2005bl ) recently revealed that small-scale fil- 
aments created in water at powers lower than 8P cr did 
not need plasma saturation. With 38-fs, 110-^jm pulses 
and k" ~ 560 fs 2 /cm, the normalized dispersion length 
is 5 4.9, for which the plasma response may indeed 
be inhibited, following Fig. [5] Pulses are then subject 
to space-frequency coupling from the interplay among 
SPM, dispersion and phase mismatching, that produces 
conical emission (see below). The energy loss caused 
by the conical wave depends on the energy co ntained in 
the co lored part of the spectrum. Liu et al. (ILiu et 



l2005bl ) suggests that if the phase matching width of light 
waves whose frequencies experience self-focusing is small 
enough compared with the spectral width of the on- 
axis laser pulse, the spectral energy going into the self- 
focusing components is small and the energy loss caused 
by conical emission becomes large. This phase matching 
width is all the smaller as GVD is large and the spread 
of energy can stop the Kerr focusing without the help of 
plasma defocusing at powers moderately above critical. 

Nonetheless, there exist different scenarios having at- 
tracted attention in the past few years. For instance, 
research groups led by P. Di Trapani and A. Dubietis 
(|Dubietis et all l2004ai [20031 : iPorras et aLl . l2004h empha- 
sized the apparent absence of free electron emission in 
the self- guiding of 527 nm, 200 fs pulses with ~ 100 
/_tm waist and powers up to 10 P cr , propagating in wa- 
ter in loosely focused geometry. Numerical simulations 
outlined the filament reshaping into a nonlinear X (bi- 
conical) wave, driven by the combination between lin- 
ear diffrac tion, normal GVD, self-focusing and nonlin- 
ear losses ( Kolesik et all l2002t) . Whereas GVD cannot 
halt the beam collapse even in these configurations (5 
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FIG. 13 7*//max versus Ao for different pulse durations t p in 
air (70 fs, dotted curve), in fused silica (20 fs, dashed curve), 
and in water (170 fs, solid curve; 50 fs, dash-dotted curve). 



is too small), normal dispersion plays an important role 
by dispersing the most intense time slices of the pulse. 
Instead of an extensive production of free electrons, the 
nonlinear losses (MPA) inherent to this process more- 
over dominate at laser wavelengths as low as 527 nm. 
From Eqs. ()76a|) . it is always possible to compute nu- 
merically the ratio I*/I max that dictates the MPA ef- 
ficiency [~ Uip n t W(I)/I] compared with plasma defo- 
cusing [~ ikop/p c ] in Eq. (|4"5j) . This ratio involves the 
quantity 



r 



nlpcUi/(k t p ) 



(78) 



that scales like n /(\ t p ). The resulting curve is dis- 
played in Fig. Q2] for fused silica, water and air at differ- 
ent pulse durations. Keldysh ionization rate for crystals 
is used for dense media, while the PPT rate is employed 
for air (O2 molecules). This figure shows that the per- 
centage of MPA over the emission of free electrons aug- 
ments all the more as the laser wavelength is low and the 
pulse duration is short. In air, MPA becomes negligible 
compared with MPI. For dense materials, Fig. IT51 points 
out to a sudden change in the influence of MPA whenever 
Ao < 600 nm. This mainly follows from the dependences 
over Ao of both the critical plasma density ([57]) and the 
photo-ionization rate W(I). As a result, self- focusing is 
stopped at lower intensities a nd the free electron d ensity 
is maintained at lower levels (|Skupin et all l2006al ) . 



B. Robustness and multifilamentation 

Femtosecond filaments are remarkably robust, in the 
sense that they can propagate over several meters and 
even reform after hitting an obstacle. In the atmosphere, 
aerosol particles, like water droplets or dust, may have di- 
mensions comparable with the filament diameter, which 
could seriously harm the dynamical balance required to 
propagate filaments. However, even opaque droplets as 
large as 100 pm in diameter do not block filamentation. 
More intriguing, energy losses are limited to 15% of the 



filament energy ( Courvoisier et al ., 2003). The filament 
rapidly self-heals and seems unaffected by the droplet, 
so that an energy balance between the filament core and 
the unfocused part of the beam was conjectured to ex- 
plain the rebuilding of the pulse. Unfocused parts of the 
beam refer to the "photon bath", i.e., the low-amplitude 
background surrounding the filament core. 

Numerical simulations constrained to the radial sym- 
metry reproduced the smallness of the immediate en- 
ergy losses in air and the interaction dynamics was sim- 
ilar whatever the pulse self-channels in the first (lead- 
ing ed ge) or in the second (trail ing edge) focusing 
stage dKolesik and Moloney!. [200 4b). From water ex- 
periments (jDubietis et aLl . l2004bf) . analogous conclusions 
were drawn, which insisted more on the linear contribu- 
tion from the outer conical waves to the filament recon- 
struction. The "active zone" of the photon bath must, 
however, keep time slices with power above critical to 
re-localize the beam at center. The filament must also 
be somehow delocalized in space, in order to justify the 
smallness of energy losses during the collision. Answer- 
ing this point, simulations revealed that the breakup of 
rotational symmetry of pulsed beams can push the fo- 
cusing components, still untouched by the plasma, out 
of axis. The defocuscd time slices diffract to the low- 
intensit y background , but t he foc using ones exit out of 
center (jskupin et all 120031 . 120021 ). A femtosecond fila- 
ment thus consists of a "spatially-extended" structure, 
whose most intense components can move within a zone 
of several hundreds of /im around the origin. Through 
a collision, just a few components are involved, which 
explains the smallness of the losses. Three-dimensional 
simulations enclosing the filament inside a circular disk 
of 300 /im in diameter clearly showed that over the 2 cm 
of self-healing the beam components outside this zone 
play no signif i cant r ole in the filament reconstruction 
(jSkupin et a/1 l2004al) . Some of these results are sum- 
marized in Fig. 1141 Recent experiments using pinholes to 
block the photon bath confirmed these theoretical predic- 
tions. The energy reservoir actively feeding the filament 
core was measured within a diameter of 220 — 440 pm 
and was found to contain up to 50% of the pulse energy 
(iLiu et aLl . l2005d) . 

Whereas a single filament appears as a robust object, 
broad pulses having many critical powers become sub- 
ject to modulational instability. The background energy 
reservoir favors energy exchanges between the different 
optical cells, that focus under Kerr compression and de- 
focus by plasma generation. A sea of collapsing spots is 
thus nucleated, forming an " optically turbulent femtosec- 
ond light guide" along which non linear dissipation con - 
sumes small energy per filament ([Mleinek et all Il999[) . 
Inside the bundle, the filaments are able to merge and 
relaunch recurrent collapse events at further distances. 
They gather into long-range clusters composed of sec- 
ondary cells that occur around the primary filaments 
created from the initial beam defects (jBerge et aZj . 12004 
ISkupin et aZ.L ^od^rj ). For broad pulses (wo > 1 cm) en- 
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FIG. 14 (a,b) 3D surface plots of plasma strings in cases of (a) 
free propagation or (b) an impact with a 95 fim large droplet 
located at z = 0.4 m (see arrow), computed from (3 + 1)- 
dimensional numerical simulations. Note the convected com- 
ponents at z ~ 0.65 m, corresponding to the second focusing 
sequence when the trail peak takes over the front pulse and 
develops azimuthal instabilities in space, (c) Temporal pro- 
files occupying an active zone of ~ 300 /jm in diameter. 



closing enough power (e.g., more than 100 -Per), the self- 
focusing distance is close to the MI filament at ion distance 
0fii ~ l/r^/p, which varies like 1/Pjn when ever I -C 

*max 

(|Campillo et all Il974t iFibich et oLI . l2005h . 



For narrower beams, experiments exploring mJ , mm- 
waisted pulses (iHosseini et al. , 2004 at iLiu et al. |, 12004 
iTzortzakis et aLl . l2001a| ) reported the beam breakup into 
two filaments that continue to increase in intensity by 
Kerr focusing. Depending on their separation distance 
and individual powers, these filaments can merge. An 
example of this fusion mechanism is shown in Fig [TS] for 
two (plasma-free) optical spots, resulting from the mod- 




FIG. 15 Merging of two filaments created from mm-waisted, 
fs pulses with 5 mJ energy. Arrow indicates a scale length of 
2 mm. 



ulational inst ability of a 42-fs, 3 mm- waisted collimated 
pulse in air ( Tzortzakis et all , l2001a[ ). In the presence 
of plasma defocusing, "parent" filaments develop spatial 
rings whose radius increases with z. Two systems of rings 
can thus interfere to yield "child" filaments. Each focus- 
ing spot triggers a plasma sequence that defocuses the 
trail pulse and locates the cells at specific negative in- 
stants. The number of "filaments" is therefore higher in 
intensity th an in the beam fluence, which sums up all the 
time slices ( Champeaux and Berget 120051 ). 

Since multifilamentation emerges from the initial beam 
defects, several techniques to control this process in shot- 
to-shot, meter-range experiments have been proposed. 
Among th ose, we can mention the introduction of beam 
ellipticity dDubietis ef"afl.l2004ctlGrow and Gaetall2005| ) 
or of a tilted lens (jFibich et ali , l2004j ) along the opti- 
cal path. Strong field gradients or forcing aberrations 
through phase masks can also be used to cont rol high 
intensities over long scales ( Mechain et alll2004) . Other 
techniques consist in monitoring the filamentation dis- 
tance by changin g the beam diver gence angle with a de- 
formable mirror ( Jin et al\ , 120051 ). By this method, a 
single filament has been produced with powers as high 
as 420 GW, usually yielding multiple filaments in col- 
limated geometry. Alternative methods consist in tun- 
ing the focal distance of a telesco pe setup composed 
of divergent an d convergent lenses ( Fibich et all l2006t 
ILiu et ad 1200(1 ) . Interesting features have recently been 
published concerning the fusion, repulsion and spiral co- 
propagation of two filaments, depending on their rela- 
tive phase shift and crossing angle, which m ay further 
be u sed to optimize the propagation range ( Xi et all 
I2006T). In the same spirit, arrays of diffra ctive microlenses 



( Cook et all l2005h or periodic meshes ( Kandidov et all 
I2005T) into the propagation path helps in making the mul- 
tifilamentation " more deterministic" . As a matter of fact, 
these methods amount to altering the background sur- 
face of the beam from which small-scale filaments can 
randomly emerge from modulational instability. Alter- 
native w ays may be provided by femtosecond optical 
vortices ( Vincotte and Bergj . 120051) that keep the fila- 
ments rotating along a ring and are nowadays designed 
experimenta ll y on f s time scales ( Marivenko et all 120051 : 
IVuong et all [2006). Further explorations should con- 
cern optical smoothing tec hniques breaking the laser spa - 
tial or temporal coherence ( Marklund and Shuklal . 12006V 
These techniques are currently employed in laser facili- 
ties devoted to the inertial confinement fusion, in order to 



homogenize the beam distribution (Do nnat et al ., 1992; 
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lLabaune et al. , 1992). 



S(9 x ,9 y ,u,z) = \Q{6 x ,6 v ,u,z) | 2 , defining 



C. White light generation and Conical emission 

Spectral broadening has be en the topic of intense inves - 
tigations for several decade s (lAlfano and Shapiro . 19701 ). 
In 1995, Gil'denburg et al. (jGiFdenburg et a/.l . ll995f) per- 
formed 2D numerical simulations of ultrashort focused 
electromagnetic waves creating breakdown plasma due to 
tunnel gas ionization. High blueshifts of the pulse spec- 
trum up to Alu/luq = lo/ujq — 1 > 40% were reported. 
At that time, blueshift was attributed to plasma genera- 
tion, as the growth of free electrons increases the plasma 
frequency by w 2 c = g 2 p(/)/m e eo and implies a positive 
frequency shift A. i v + ~ fn dtodl (iPenetrante et a/.l . [l992t 
iRae and Burnettl Il992t lYablonovitchl . Il974fh Later, ex- 
perimental spectra revealed an " ultrafast white-light con- 
tinuum" . This occurs in a wide variety of condensed and 
gaseous media, whenever the input beam power exceeds 
the self-focusing threshold, which fix es a striking "band 
gap" for su percontinuum generation ( Brodeur and Chinl . 
|1998| . |1999D . Around the central wavelength of 800 nm, 
blueshifts (Aw + ) down to 0.5 /xm in wavelength enlarged 
the spectra at intensities lower than 10 -2 times the cen- 
tral one. Spectral enlargements appeared asymmetric, 
with a limited redshift, which could not be explained by 
SPM theory alone. Continua emitted from gases were 
found narrower ( Aw + ~ 0.5 ojq) than those for dense me- 
dia (Auj+ ~ 0.8 -1.6 w ). 

Because this process mixes all spectral components in 
the visible range, the core of the filament evolves like 
a white spot, whic h gives rise to a "white- light laser" 
(|Chin et all [l999b). Noting by <p(r,t) the phase of the 
field envelope, frequency variations are dictated by 



Q(6 x ,6 v ,gj,z)= / dxdydtU(x,y,z,t)e iiAjt - i9 * kx - ie y kv , 



Auj 



-d t ip ~ -k Q Azd t {n 2 I - p/2n p c ), 



(79) 



which varies with the superimposed actions of the Kerr 
and plasma responses. Near the focus point, only the 
front leading edge survives from this interplay and a red- 
shift is enhanced by plasma generation. At later dis- 
tances, second focusing/defocusing sequences relax the 
spectra to the blue side while the redshifted components 
decrease in intensity. A salient blueshift around the cen- 
tral frequency follows from self-steepening that creates 
a shock edge at the back of the pulse and amplifies a 
"blue shoulder" in the spectrum ( Akdzbek et all [2001V 
This spectral dynamics readily follows from Eq. (|70p . 
Asymmetries caused by the focusing events are impor- 
tant: They design the pulse spectrum which will be pre- 
served after the nonlinear stage of the beam. 

The pulse spectrum is usually computed from the 
Fourier transform of the field intensity, either at r = 0, 
or being averaged over the filament radius, Sn\(ui) — 
2tt fy hl \£(r, w)\ 2 rdr. More insights can be obtained 
by looking at the frequency-angular spectral function 



(80) 

where 9 X = k x /k, 9 y = k y /k are the angles at which dif- 
ferent frequency components uj propagate in the medium. 
The dependence of the angle 8 = 9 X = 9 y on wavelength 
defines the supercontinuum cone emission and is obtained 
by means of the expression A = 2ttc/uj. Surface plots of 
lnS^A,^) compared with the spatio-temporal intensity 
distributions figure out a conical emission (CE) driven 
by "wings" directed to non-zero angle s and connected 



with the p lasma-induced spatial rings (Kandido v et al. 



|2004L [2003). Nibbering, Kos areva et al. (jKosareva et al. 

], |l996) first revealed the existence 



119971 : iNibbering et al. 



of this conical emission accompanying femtosecond fila- 
mentation in air at 800 nm. Recalled in Fig. [TBT a). CE 
increases in the interval < 9 < 0.12° with decreas- 
ing wavelengths 800 > A > 500 nm and is independent 
of the position along the filament. This "conical emis- 
sion" corresponds to the angular divergence of spectral 
components. The (half-) angle 9 at which the radiation 
propagates is determined by the spatial gradient 

d 

6{r,t) = arctanf-fcn" 1 — ip(f,t)]. (81) 
or 

This angle depends on the spatial intensity distribution 
of the pulse, which strongly evolves in self-channeling 
regime, as illustrated in Figs. [I6jb,c,d,e). By compar- 
ing the radial dependences of the divergence angle 6(r), 
low-frequency components are located nearest the optical 
axis. High-frequency components lie near the periphery 
rings and their radiation is directed out of the axis. 

Besides Kerr and plasma effects, the material disper- 
sion can play a relevant role in condensed media, for 
which the GVD coefficien t is usually high. Kolesik et 
al. (jKolesik et all 2003a, b) compared the extension of 
supercontinuum wings generated in air (k" = 0.2 fs 2 /cm 
at 800 nm) and in water (k" = 500 fs 2 /cm at 400 nm). 
Numerical simulations accounting for the complete linear 
dispersion relation displayed evidence of that chromatic 
dispersion becomes a key player (for water above all) in 
determining the spectral extent of supercontinuum gen- 
eration: Full chromatic dispersion included in k(u>) en- 
compasses the T, T _1 operators when passing to an en- 
velope description (see Sec. [TTJ) . These operators induce 
steepening of the pulse profile and strongly "blueshift" 
the spectra. 

On the other hand, spectral broadening becomes en- 
hanced by harmonic generation. The question of cou- 
pling ultrafast IR components with self-induced third - 
harmonics was raised ten years ago (jBackus et all Il996h . 
While conversion efficiency as high as 0.1% was first re- 
ported for third-harmonic (TH) generation in air, higher 
efficiencies up to 0.2% were next achieved by using fo- 
cused pulses filamenting over only ~ 10 cm. At 800 
nm, the coupling of TH with the IR pump produces 
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FIG. 16 (a) Measured (symbols) and calculated (curve) CE 
angles in air. (b,c) Frequency angular spectra and (d,e) 
spatio-temporal distributions ol an atmospheric, fs filament 
starting from infrared pulses with 6.3 critical powers, 250 
fs FWHM duration and 0.18-mm radius at two different 
ropagation distances (b,d) z — 0.43zo; (c,e) z — 0.6zo 
Kandi dov et aZX l2003t ) . Courtesy of O. Kosareva. 



a "two-colored" femtosecond filament f rom the thresh- 
old intensity ^ 2 x 10 13 W/crrr* (lAkozbek et all 



l2002HAlexeev et all 120051 ; lYang et al. , 2003). Along this 
process, the pump wave injects part of its energy into 
the third harmonics. The amount of pump energy de- 
pends on the linear wave vector mismatch parameter 
Afc = [3k(ui) — fc(3w)] _1 fixing the coherence length 
L c = n/\Ak\. The smaller the coherence length, the 
weaker TH fields, since most of the pump energy is pe- 
riodically injected and depleted by the pump. In self- 
channeling regime, the balance between TH, pump wave 
nonlinearities and the linear mismatch parameter makes 
the phase difference Aip = 3ip(u>) — ip(3u!) be clam ped at a 
constant value equal to ir ( Akozbek et aLl . [2002h . Along 
meter-range distances, the two-colored filament propa- 
gates over longer scales than an IR pulse alone, as the TH 
component stabilizes the pump wave an d achieves 0.5% 
conversion efficiency (jBerge et all l2005f ) . Experimental 
identification of the TH component reveals some central 
spot surrounded by ring structures embarking most of 
the T H energy and haying a h alf-divergence angle of 0.5 
rnr ad (|Theberge et all l2005af) . This contributes to en- 
hance the total conical emission of the beam. 

Harmonic generation permits the occurrence of amaz- 
ing phenomena, such as, e.g., the emergence of new wave- 
lengths with non-trivial relative spectral intensities. Ma- 
nipulating intense, powerful filaments centered at 800 



nm results in supercontinuum generation extending down 
to 230 nm in air and yieldin g a continuous spectra l 
band of UV- visible wavelengths ( Theberge et all l2005bl) . 
This phenomenon was observed over laboratory scales 
as well as over LIDAR propagation dista nces > 200 
m (|Berge et all 120051 : iMeiean et all l2006bj ). Numerical 
simulations evidenced that UV-visible spectral broaden- 
ing is created by the overlap of redshift from the TH com- 
ponent and blueshift from the IR pump, Au>j — —dtfj 
(j = oj, 3u>), fulfilling the phase-locking constraint Aip = 
7T mentioned above. Also, TH-induced saturation rein- 
forces the defocusing action of the electron plasma and 
enlarges the conical emission to ~ 0.25° at 0.26 /j,m. 
UV-visible supercontinuum generation in air is plotted 
in Figs. I17f a,b) for a meter-range propagation of one 
filament. Numerical results reproduce the experimental 
data in the wavelength domain 200 < A < 500 nm at the 
same distances after the nonlinear focus (Az = z — z c ). 
Fig. ll7f c) illustrates the maximum TH and pump intensi- 
ties, calculated from a propagation model accounting for 
TH generation [see Eqs. (|57)) ] without steepening terms. 
Resembling spectral dynamics have also been reported 
from 45-fs, 1-mJ infrared pulses propagating in argon 
at atmospheric pre ssure, after subsequ ent compression 

. Simulations of 
were performed 
51) in the limit of 



by chirped mirrors dTrushin et 




these experiments (Akozbek et al. 
from the one-component NEE model 
weak THG. It was found that temporal gradients inher- 
ent to the T, T^ 1 operators are sufficient to amplify UV 
shifts and cover the THG bandwidth for very short pulses 
(< 10 fs) over short conversion lengths that involve a 
limited number of focusing/defocusing cycles. Huge UV 
shifts have already been displayed for air propagation in 
Fig. Illf d) . The issue of knowing under which conditions 
steepening effects or THG prevail in the supercontinuum 
generation remains open. 



D. Role of the laser wavelength 

Whereas emphasis was put so far onto optical pulses 
delivered by Ti:Sa oscillators operating around 800 nm, 
one could wonder whether the propagation dynamics 
changes at different laser wavelengths. Besides the 
Rayleigh length zq, several physical parameters become 
modified when varying Ao, namely, 

• The Kerr index 77,2, function of the cubic suscep- 
tibility , decreases inversely propor tional to Aq 



tibility \i , r; , decreases inversely propor 
( Agrawall l200ll ; iNibbering et all Il997f ) 



The photo-ionization rate and the number of pho- 
tons consumed by MPI transitions ~ Ui/fkuo vary 
with Ao- Even though MPI gain and MPA losses 
cannot be separated from each other, their respec- 
tive influence changes with Ao and stresses more 
losses at low wavelengths (see Fig. [T3l . 

The value and sign of the GVD parameter, k" = 
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FIG. 17 UV-visible supercontinua (a) measured in laboratory, 
(b) computed numerically, (c) TH (dotted curve) and pump 
(solid curve) intensities integrated numerically. The dashed 
curve represents the same quantity computed using the pump 
component alone. 



d 2 k/duu 2 \ w=iJ j , are modified by varying Ao- De- 
pending on the material considered, k" can turn 
from positive (normal GVD) to negative (anoma- 
lous GVD), which deeply modifies the spatio- 
temporal structure of the beam. 

On this topic, propagation studies mostly concerned 
the differences between infrared pulses at 800 nm in the 
atmosphere and ultraviolet ones at, e.g., 248 nm, cre- 
ated by an hybrid, frequency-tripled Ti:Sa/KrF excimer 
chain. Long distance p ropagation exp e rimen t s em ploy- 
ing ps, 248-nm pulses (jSchwarz et all l200ll 120001) re- 



ported that after a transient stage of multiple filamen- 
tation, the beam relaxed to a single filament upon prop- 
agation ranges of more than 10 m. The filament length, 
defined as the spatial region where the plasma couples 
with the light channel, only covered ~ 1 m, where the 
electron density, evaluated from the energy consump- 
tion per filament, attained 3 x 10 15 cm" 3 . Computed 
from the MPI rate, the peak intensity was deduced to 
be about 10 12 W/cm 2 . Also associated with this fila- 



ment, a narrow spectral bandwidth of AA ~ 10 nm was 
measured around 248 nm. Lat er, Mysyrowicz's group 
(|Tzortzakis et all l2000aL l2001b[ ) observed the filamen- 
tation of UV pulses with durations from 5 ps down to 
450 fs. Long, meter-range (~ 2 — 4 m) filaments coupled 
with 10 15 — 10 16 cm" 3 peak electron densities at inten- 
sities of ~ 10 11 W/cm 2 were identified. The "apparent" 
differences between UV and IR peak values follow from 
the pulse duration between UV and IR experiments, the 
plasma critical density p c ~ 1/A 2 , together with the num- 
ber of photons involved in multiphoton transitions. This 
varies from K — 8 at 800 nm to K — 3 at 248 nm for 
O2 molecules. Numerical simulations, however, showed 
that, for 50-fs FWHM input pulses, the filament dynam- 
ics and length (< 5 m), starting from the same nonlinear 
focus, were identical in both cases with peak intensities 
and electron densities attaining 4 — 8 x 10 13 W/cm 2 and 
10 17 cm -3 , respectively. Even though the ionization pa- 
rameters vary with Ao, the peak values are fixed by the 
balance p ma x/2p c ~ which remains similar at the 

two wavelengths. Averaging the peak quantities over a 
given filament width ca n justify the different measure- 
ments made in this field (Couairo n and Bergetl2002D . 

Besides, changes in supercontinuum generation were 
nume rically addressed by Akozbek et al. (jAkozbek et all 
120031) . who accounted for THG at other laser wavelengths 
such as 1064 nm and 1550 nm. The latter corresponds 
to eye-safety and is particularly important for LIDAR 
applications. For these wavelengths, the GVD parame- 
ters for both pump and TH frequencies are close to each 
other. However, the wavevector mismatch strongly de- 
creases from —5 cm" 1 (800 nm) to —0.65 cm -1 (1550 
nm), whereas the temporal walk-off parameter Aw = 



(cj)] 1 increases from 0.4 to 2.0 cm/s, 



respectively. Consequently, the coherence length sig- 
nificantly augments at large wavelengths, which allows 
higher TH intensities and energy conversion efficiency. 

More surprising features occur when the laser wave- 
length is selected in such a way that k" becomes neg- 
ative and leads to anomalous dispersion. In that case, 
the pulse ceases to be dispersed in time and undergoes 
a temporal compression, in addition to the spatial Kerr 
focusing. This happens, e.g., in fused silica at Ao = 1550 
nm, for which k" ~ —280 fs 2 /cm (see Table I). Re- 
cent experiments at this wavelength have shown that col- 
lapse events looked "extended" along the z axis, unlike 
the localized events promoted by normal GVD. Plasma 
halts the collapse at powers above critical, but anomalous 
GVD continues to transfer energy into the collapse re- 
gion, resulting in the formation o f longer filament s befor e 
the beam eventually defocuses ( Moll and Gaetal . 12004). 
The pulse can thus remain confined along several diffrac- 
tion lengths and develops very narrow, isolated temporal 
peaks. A strong temporal compression produces opti- 
cal spots whose duration is shrunk to the few-cycle limit 
( Berge and Skupinl . l2005h . Sharp shock profiles emerge 
from this 3D collapsing dynamics and pulse steepening 
along each focusing/defocusing event tremendously am- 
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FIG. 18 Temporal dynamics of 42.5-fs pulses with waist too = 
71 /zm propagating in silica at the laser wavelength (a) Ao = 
790 nm with power Pi n — 6P cr , (b) Ao = 1550 nm with Pi n = 
3P cr . Note that the physical value of P CI changes with Ao- 



plifies the blue part of the spectrum. Fig. [18] displays the 
temporal distributions of a femtosecond filament created 
in fused silica at 790 nm (k" — 370 fs 2 /cm > 0) and at 
1550 nm (k" < 0). At 790 nm, normal GVD stretches 
the pulse along the time direction. Reversely, at 1550 nm, 
anomalous GVD compresses it temporally and makes it 
shift to positive instants through third-order dispersion 
and self-steepening that push the pulse centroid towards 
the region t > 0. 



V. ULTRASHORT FILAMENTS IN GASES 

This section is devoted to novel trends opened by the 
self-guiding of femtosecond pulses over small distances 
(< 1 m), such as pulse shortening and high-harmonic 
generation. Most of these applications are conducted in 
pressurized gas cells, in which the raise of plasma excita- 
tions can be controlled by changing the medium pressure. 



A. Novel perspectives for pulse shortening 

In the last decade, an important challenge has con- 
sisted in testing powerful techniques of spectral broad- 
ening and ultrabroadband dispersion control, in order 
to compress pulses to durations of a few optical cycles 
only ( t c. = A p/c). Nisoli and co-workers ( Nisoli et ali , 
Il997al la. ll998f ) achieved the shortest pulse durations from 
~ 1 mJ, 20-fs pulses by manipulating the spectral broad- 



ening attained along a 60-cm long fused-silica hollow fiber 
filled with atomic (Ar, Kr) or molecular (N2) gases. Con- 
ditions for optimum pulse compression outlined an ap- 
propriate combination of SPM and gas dispersion man- 
aged by external chirped compensation systems. Pulses 
as short as 4.5 fs with output energies up to 70 fiJ could 
be delivered at 800 nm and 1-kHz repetition rate after 
using a chirped-mirror delay line. The potential scalabil- 
ity of this system to higher pulse energies was claimed 
to hold provided that two conditions are fulfilled: (i) 
The laser peak power must be smaller than P cr ; (ii) The 
peak intensity must not exceed the ionization threshold, 
to preserve a flat spectral phase for an efficient recom- 
pression. These two requirements could at that time be 
satisfied by a delicate tuning of the pressure parameter. 
For a gaseous medium with pressure p, the coefficients 
affecting Eq. (|45ap indeed vary as follows 



k" 



P, 



/3 ( 



K) 



P, 



l/n 2 ~ Pnt~P- 



(82) 

An optimal compression length can be estimated by 
L opt « (Oi-NLidisp) 1 / 2 where L disp ~ t 2 p /k" and 
Lnl — n2UJoPi n /cS e s are the dispersion and nonlinearity 
lengths, respective ly, with S e s being the effective area 
of the beam mode (jNisoli et aZ.I . Il997ah . This distance is 
tuned by the cell length and the local pressure p, in order 
to reach the smallest possible pulse duration. Stability of 
the beam was insured by selecting the EHn hybrid mode 
with intensity I(r) = Iq Jq (2.405r/a), where Iq is the 
peak intensity, Jo is the z ero-order Bessel func t ion an d 
a is the capillary radius ( Tempea and Brabed . [l998b). 
Keeping the beam stable inside the waveguide requires to 
limit the excitation of higher-order modes owing to the 
nonlinear spatio-temporal dynamics. The spectral broad- 
ening thus inherently depends on how close the initial 
peak power of the beam is to P cr ( Homoelle and Gaetal . 
l2000f ). Sub- 10 fs pulses are routinel y prod uced through 
these techniques ( Steinmever et ali Il999l ). They can 



even be applied to bulk (BK7 glass) media, from which 
pulse shortening by a facto r 3-5 and output e nergies of 
220 /zJ have been achieved (jMevel et aLl . l2003D . 

In this scope, plasma generation may not constitute 
a drawback for pulse compre ssion. Tempea and Brabec 
( Tempea and Brabed . Il998al) mentioned the possibility 
to produce few-cycle optical pulses by compensating the 
spectral chirp generated by the plasma nonlinearity with 
a dispersive line delay for intense pulses propagating in 
filled-gas hollow fibers. This scenario was experime ntally 
confirmed by a JILA team ( Wagner et ali , \2004 ). who 
succeeded in compressing 30 fs pulses to 13 fs in hol- 
low argon capillaries, even without any need of external 
dispersion compensation. It was recently improved to de- 
liver 5 fs, high energy (1 mJ) pulses following the same 
scheme supplemente d by amplifier s ystem s and chirped- 
mirror compressors (jVerhoef et aZ.I . I2006T ). In Sec. IIV1 



Fig. [9fd) recalled the basic property of pulse compres- 
sion through the combined effects of Kerr nonlinearity 
and self-induced ionization. Chromatic dispersion and 
nonlinear losses, however, saturate the defocusing ac- 
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tion of free electrons on the leading edge and enable 
the rear pulse to refocus at later distances. This pro- 
duces a two-peaked temporal distribution whose time 
extent is often of the order of the input pulse dura- 
tion. Achieving an efficient pulse compression, instead, 
requires to isolate one (at least domin ant) peak shrunk 
in time ( Champeaux and Berget |2~003( ). 

To realize this challenge, Hauri et al. used a configura- 
tion of two cascaded gas-filled cells with intermittent dis- 
persion compensation for producing sub-m J light pulses 
with durations down to 5.7 fs at ~ 800 nm (jHauri et all . 
I2004D . The cells were operated at different pressures. It 
was later suggested that this scheme could be improved 
when a single cell with a pressure gradient is used instead, 
in order to monitor the refocusing stage of the trail and 
create waveforms with a sin gle-cycle temporal signature 
( Couairon et all . [20061 [2005). De spite first results m ainly 
directed to higher pulse energies (jSuda et "aZl. l2005h . fur- 
ther efforts are needed to adjust the pressure gradient for 
optimum compression to extremely short pulse duration. 
Alternatively, Stibenz et al. recently demonstrated an ef- 
ficient pure self-compression obliterating the need for any 
kind of dispersion compensat ion, pressure gradie nts, or 



2006). The 



capillaries for beam guiding (jStibenz et all 
experimental setup involves beams in convergent geom- 
etry with energy up to 5 mJ, Wo = 11 mm, / = 1.5 
m, FWHM durations of 45 fs and input powers equal 
to 5 times critical. The beam is initially focused and 
propagates inside an 1-m long, 50 k-Pa argon cell whose 
center is positioned at the location of the geometrical fo- 
cus. Output pulses reach durations down to about 10 
fs and develop a strong blue-shift, as reported in Figs. 
[TW a.b) . FigurefTWc) details this self-compression mecha- 
nism, reproduced numerically from Eqs. ([45]) . First, near 
z c ~ 1.45 m, the front pulse focuses and plasma depletes 
its rear part. Second, at the linear focus (z ~ 1.5 m), the 
plasma density decreases, and the back pulse rises again. 
Finally, from z ~ 1.6 m, only the components close to 
center, keeping intensity levels above 10 13 W/cm 2 , are 
effectively trapped in the filament, whereas the temporal 
wings diffract rapidly. A robust, temporally-compressed 
structure of 11 fs forms in the core region of the filament 
[Fig. nW d)]. The pulse center relaxes to a narrow "waveg- 
uide" , preserving its energy at low plasma density levels 
(< 10 14 cm -3 ) after the linear focus. Fig. flWe) shows 
the on-axis intensity spectra, corresponding to this tem- 
poral compression. Not shown here, the spectral phase 
at the cell exit [z = 2 cm) becomes nearly flat in the blue 
wavelengths with tiny variations comprised in the same 
proportions as those measured by spectral phase interfer- 
ometry for direct electric-field reconstruction (SPIDER). 
The robustness for this self-compression mechanism was 
numerically found limited to rather long f /wq ratios and 
powers ar ound 5 P cr , in order to avoid another focusing 
sequence ([Skupin et all , l20"06bl ) . 

T o end with this as pect, let us mention that Chen et 
al. (|Chen et all . l2006h also succeeded in achieving pulse 
self-compression from 50 to 20 fs, by making 800 nm 
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FIG. 19 (a) Experimental pulse shape for the self- 
compression mechanism in a 50 k-Pa argon cell, (b) associated 
spectrum, (c) On-axis numerically-computed temporal evo- 
lution of 45-fs pulses in the (t, z) plane for a similar setup, 
(d) On-axis intensity and (e) spectral intensity at z — 1.6 m 
(dashed curve) and z = 2 m (solid curve, exit window of the 
argon cell). 



pulses with energy < 1 mJ pass through a thin BK7 
glass plate, with no subsequent dispersion compensation. 
SPIDER and spectral measurements revealed that self- 
compression takes place as long as the front pulse focuses 
while plasma cuts off the trail part. At this stage, the 
spectrum is marked by a strong redshift. By increasing 
the input peak intensity, the back pulse refocuses and 
blueshifts the spectrum. Accounting on the balance be- 
tween plasma-induced negative chirps and SPM-induced 
positive chirps, the spectral output phase was, again, 
observed to be flat. All these methods represent new 
promising techniques to produce high-energy few-cycle 
laser sources in the future. 



28 



B. High-order harmonic generation 

The possibility to reduce pulse durations to a sin- 
gle optical cycle provides ideal conditions for high- 
harmonic generation and laser - assisted X-ray photo- 
ionization ( Drescher et all [2001). Several fundamental 
atomic processes such as inner-shell electronic relaxation 
or ionization take place within a fraction of the oscilla- 
tion period of visible or near-infrared radiation and they 
require very short probes for being investigated. For this 
purpose, bursts of attosecond pulses (1 as = 10~ 18 sec.) 
need to be isolated. This is now feasible by combining 
single bursts emitted at extreme ultraviolet (XUV) wave- 
lengths and ultrashort lasers, in order to preserve the co- 
herence properties at so short time scales. In this micro- 
cosm, the ionization processes become crucial, and more 
particularly the transient stage of about half laser period 
along which an electron is liberated. Developments in at- 
tosecond X-ray science have been addressed in several ex- 
cellent reviews (IBrabec and Krauszl l2000t iPfeifer efal\ . 



120061 : IScrinzi et all 12006 ). Subfemtosecond light pulses 



can be experimentally obtained by superposing several 
high harmonics and making them emit simultaneously. 
By controlling their synchron ization, pulses of 130 as in 
duration have been achieved (jMairesse et al. I. I2003D . 

High-Harmonic Generation (HHG) describes the pro- 
cess by which laser light at central frequency ujq is con- 
verted into integer multiples of u>o during the highly non- 
linear interactio n with the medium. In 1993, Corkum 
( Corkuml . ll993f) proposed a quasiclassical theory for this 
process, which can be divided into three steps: 

(i) The electron is freed by ionization and driven away 
from the parent ion. 

(ii) Since the laser field changes its sign over times 
about r o c ./2, this electron slows down, stops at a certain 
position, then re-accelerates towards the ion. 

(iii) When the electron recombines with the nucleus, a 
photon with energy equal to Ui plus the electron kinetic 
energy is emitted. This gives rise to very high harmonic 
orders. 

Figure [20] illustrates a typical HHG spectrum. The 
harmonic intensity is rapidly decreasing after an almost 
flat plateau (A < 15 nm). It is terminated by a cut- 
off, signaling the highest harmonics that can be gener- 
ated. Among all possible electron trajectories, there in- 
deed exists a maximum kinetic energy corresponding to 
the maximum photon energy embarked by this process. 
The numerical solving for the electron trajectories with 
various phases leads th i s max i mum kinetic ener gy to be 
- 3.17 x U p (ICorkuml . Il993l : IZeng et all 120031) . Here, 
Up = q^Ep/Am e ujQ denotes the ponderomotive energy of 
the electron in the wave field (see Appendix [3} . Elec- 
trons thus produce harmonic photons up to the cut-off 
energy 



huj r = 3.17 x U„ 



Ui. 
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FIG. 20 Harmonic spectrum generated in neon by a 800 
nm Ti:Sa laser. 25th to 8 1st harmonics are displayed 
(Sali eres and Lewen stcin. 200l|). 



E(t) = E cos (cu t + (Pce), where the carrier-envelope 
phase fcE becomes an important parameter for the con- 
versio n efficiency into HHG when approaching the optical 
cycle (|Scrinzi et all l2006t l. 

In ionized media with symmetry inversion, harmonic 
peaks only exist at odd integers multiple of u . The 
above three-step process repeats every half-cycle of the 
laser field, so that the Fourier sum in the spectra 
makes even harmonics cancel each other along one cy- 
cle. The analytical description of this process requires a 
quantum-mechanical approach of the dipole moment as- 
soci ated with the harmonics . To this aim, Lewenstein et 
al. (jLewenstein et all . 1 1994 ) derived the time-dependent 
dipole moment d(t) = q e (ip(r, t)\r\ip(f, t)) that represents 
the expected position of the electron in the quantum state 
\ip(f,t)}. Details for computing d(t) have been given in 
Appendix |Bj The time-dependent dipole takes into ac- 
count ionization at time t', energy gain computed from a 
phase integrand and recombination at time t. Stationary 
phase conditions yield the information about the electron 
trajectories. By Fourier transforming the dipole moment, 
the harmonic spectrum can be calculated and it restores 
the cut-off law (j8"3"j) . 

The challenging issue is then to describe HHG by ac- 
counting for the propagation physics of ultrashort pulses 
havi ng initially a few fs durations. In 1998, Lange et 
al. ( Lange et all Il998ah employed self-guided fs pulses 
to produce HHG up to the 15th harmonics in noble gases 
and maintain the longitudin al coherence of the g enerated 
harmonics. Tamaki et al. (ITamaki et all Il999l) demon- 



The laser pulses employed in HHG experiments are 
generally short (< 40 fs), with the elementary profile 



strated experimentally efficient HHG through the phase- 
matched propagation of laser beams in pressurized Ne 
with peak intensities up to 10 15 W/cm 2 . By changing the 
propagation length, phase-matching magnifies the con- 
version efficiency around the 49th har monics by 40 times 
near the cut-off region. Tosa et al. ( Tosa et all . |2003[ ) 
reported the H13 harmonics triggered by the self-guided 
propagation of a focused pulse inside a 14-cm long cell 
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FIG. 21 Harmonic spectrum computed from the dipole mo- 
ment d(t) in pressurized argon for two different pump inten- 
sities. 



filled with Xe, in which the Kerr response was negligible. 
From these results, we easily guess that the next step in 
" attosecond" investigations should concern the optimiza- 
tion of HHG spectra by means of femtosecond filaments. 
The promising feature of pulse shortening induced by 
the balance Kerr/plasma could become an excellent tech- 
nique to produce an isolated XUV pulse. Moreover, the 
ability to keep the laser pump clamped at sufficient high 
intensity over several centimeters may magnify the har- 
monic signal. For instance, figure |2~T1 illustrates the am- 
plitude of the atomic dipole computed numerically for 
a plane-wave pulse traversing an argon cell pressurized 
at 60 kPa for two different pump intensities. The HHG 
spectrum becomes magnified at intensity levels close to 
100 TW/cm 2 . This preliminary result keeps us confident 
with obtaining harmonic signals with high enough energy 
level. 



VI. ULTRASHORT FILAMENTS IN DENSE MEDIA 

Below, we review major physical aspects of ultrashort 
filaments in solid dielectrics and in liquids. 



A. Pulse propagation in dielectrics and damages 

Transparent dielectrics such as SiC>2 samples were rou- 
tinely examined during the 90s, in the framework of the 
optical breakdown on the one hand (see Sec. IIV[) and of 
the carrier trapping dynamics of band-gap crystals on 
the other hand. One of the key parameters being the 
potential gap (which is the energy difference between va- 
lence and conduction band), several investigations were 
led to identify this gap in silicon dioxide and to mea- 
sure the ultrafast excitation and relaxation of an elec- 
tron gas pu mped into the conduction band by int e nse fs 
laser pulses (lAudebert et all 11994 iMartin et all . Il997t 
iTohmon et all Il989h . 

Meanwhile, other researches concerned laser-induced 
breakdown (LIB) in dielectrics. The threshold damage 



fluence varying as ^/t^, for long pulses in thermal con- 
duction regimes was seen to deviate from this scaling 
with fs pulses, for which electrons have no t ime to ef- 
ficien tly couple to the lattice. Du et al. ( Du et all , 
1994 ) reported a damage threshold fluence higher than 
the yJT^ prediction rule for pulse durations < 10 ps and 
they furthermore underlined that short-pulse damage ex- 
hibits a deterministic nature, unlike long pulses. A the- 
oretical model accounting for MPI, avalanche ionization 
and Joule heating stressed t he dominant role o f photo- 
ionization at fs time scales ( Stuart et all . Il996h . How- 
ever, both photo-ionization and avalanche come into play 
during the occurrence of damag e. Although the rol e of 
the former is initially dominant ( Ravner et all . |2~005[ ) . it 
can be masked by avala nche from an init ially-high elec- 
tron density (> p c /100) (jTien et aLl . ll999l ). To complete 
these two processes, the rapid electron decay over ~ 100 
fs scales through recom bination soften s the maximum 
electron density for LIB ( Li et all . Il999t ) . 

Investigations mixing propagation and damage in di- 
electrics started f r om the early 2000's. Tzortzakis et al. 
( Tzortzakis et all . l2001df ) observed the 1-cm long self- 
channeling of a 160-fs focused pulse conveying 3 critical 
powers at 800 nm. The pulse traveled across a fused sil- 
ica sample like a narrow waveguide with ~ 20 /im waist. 
Local heating caused by the high repetition rate (200 
kHz) and damage by local intensities > 10 TW/cm 2 were 
avoided by making the sample move in the (x, y) plane. 
Auto-correlation traces revealed a two-peaked structure 
and spectra were asymmetrically broadened, which sig- 
naled a pu lse splitting driven by plas ma. Earlier, Brodeur 
and Chin (|Brodeur and Ch in. 1999) noticed a similar su- 
percontinuum in glasses and explained it from the multi- 
photon excitation of electrons into the conduction band 
at the focus point. The low beam divergence was at- 
tributed to the Kerr-lens effect. As an example, Fig. 
2"27a) shows a transverse photograph of the self-guided 
filament measured at input energy of 2 /zJ in fused sil- 
ica. Fig. I2"27 b) details the filament waist in nonlinear 
propagation regime and Fig. I22f c) reproduces this waist 
numerically computed from Eqs. (|45[) . 

At higher local intensities, the breakdown limit is ap- 
proached. Because the plasma generated during break- 
down remains at the threshold, a precise control of the 
interaction region can be reached with ultrashort pulses. 
This property can be used for material processing, medi- 
cal laser applications as well as solid-state microelectron- 
ics. The major qualitative differences between damages 
caused by short (< 10 ps) and longer pulses (> 50 ps) 
appear in the damage morphology: Short pulses ablate 
the material, whereas long pulses produce conventional 
melting, boiling and fracture. With tightly focused fs 
pulses, permanent bircfringent structures embedded in 
bulk fused silica have b e en re a lized over microscopic di- 
mensions (jSudrie et all . 1 19991 l200lh . Arrays of paral- 
lel grooves formed transmission diffraction gratings with 
periodically-modified index changes. Two distinct types 
of damages can be created, namely, those consisting of lo- 
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FIG. 22 (a) Transverse photograph of a single filament propa- 
gating (from left to right) in fused silica, (b) Measured diam- 
eter of the beam comparing high (linear) and low (nonlinear) 
divergences, (c) Numerically-computed diameter of the non- 
linear filament (jTzortzakis et all l2001dh . 



cal isotropic increase of the refractive index at subcritical 
powers, and those inducing local birefringence at powers 
above critical. The damage track consists of diffusing 
zones of ~ 20 p:m transverse width. The total dam- 
age region can extend over a distance reaching 80 /im 
along the optical path. Numerical simulations accounting 
for photo-ionization, avalanche, electron recombination 
and Ke rr self-focusing reproduce d experimental damage 
tracks (jSudrie et all |2002t l200lh . along which electron 



densities up to 3 x 10 cm" J resulted from the com- 
bined actions of photo- and impact ionization. Account- 
ing moreover for Ohmic heating an d thermal cooling by 
collisional ionization, Pehano et al. (jPehano et al 
numerically solved the ID electromagnetic wave equation 
to quantify the transmission, reflection and absorption of 
100-fs laser pulses by a thin plasma layer formed at the 
surface of a dielectrics. With a fluence of 3.2 J/cm 2 , 
the interaction fully ionizes the surface and heats the 
plasma to ~ 10 eV. A significant transmission up to 30% 
of the laser energy is possible even when the plasma den- 
sity is above critical (w po > ujq). In the layer, the elec- 
tron density rapidly becomes supercritical and implies a 
sharp rise in the reflected pulse amplitude. Compared 
with MPI alone, collisional ionization increases the peak 
plasma density by 30%. 



B. Pulse propagation in liquids and applications 

Similarly to LIB in solids, the optical breakdown in flu- 
ids gained considerable interest, because it finds various 
therapeutic applications for, e.g., plasma-mediated laser 
surgery and prevention of ocular damages. First simula- 



tions in water ( Feng et al\ , Il997f) emphasized the domi- 
nant role of cascade ionization for long pulses, while mul- 
tiphoton ionization was expected to prevail for shorter 
pulses. Further studies (jNoack and Vogel Il999h evalu- 
ated the absorption coefficients during the LIB process 
for 100 ns down to 100 fs pulses. It was established that 
a critical density threshold for LIB was p CI — 10 20 cm' 3 
for long (ns) pulses, but p CI = 10 21 cm" 3 for short (ps) 
ones. In water and for 100-fs durations, LIB results in 
bubble formation supported by thermoelastic effects at 
pulse powers above the self-focusing threshold. This pro- 
cess also causes refractive index changes in the beam path 
upstream over several hundreds of microns. Within wa- 
ter droplets the tight focusing and the nonlinearities of 
the LIB process moreover generate a nanosized plasma 
hot enough to emit in the visible, and pr eferentially in 
the backward direction ( Favre et all I2002T) . 

The first experiments on fs pulse self-guiding versus 
the optical breakdown limit in water were performed 
by Liu and co-w orkers for different focusing geometries 
dLiu et aLl . l2002h . It was found that the shorter the fo- 
cal length, the larger the transverse size of the optical 
breakdown plasma is. Self-focusing drives the initial lo- 
calization of the plasma towards the beam axis. White 
light is generated along a short (< 1 mm) filament and 
deflected at a small constant angle only. The supercon- 
tinuum sources for the high-frequencies were identified in 
the rings formed by plasma defocusing and amplified by 
the back of th e pulse, where shocklike dyna mics blueshift 
the spectrum (jKandidov et all |2004 [20031 ) . 

Far below LIB limits, the self-guided propagation of 
femtosecond filaments in water was thoroughly exam- 
ined by Dubie tis, Di Trapan i and their collaborators a 
few y ears ago ( Dubietis et all l2004al . 120031 : iPorras et all 
|2004| ). By launching a ~ 3/iJ, 170-fs clean beam with 
~ 100 /im FWHM diameter onto a water-filled cuvette in 
loosely focused (/ > 5 cm) geometry, a single filament 
formed at the wavelength of 527 nm and was capable of 
covering up to 4 cm along the propagation axis, while 
keeping a mean FWHM diameter of a few tens of /jm. 
The experiment revealed that the filament dynamics was 
not sustained by a balance between Kerr-induced self- 
focusing and plasma-induced defocusing. Numerical sim- 
ulations outlined, instead, the spontaneous reshaping of 
the beam into a Bessel-type A- wave fulfilling the require- 
ments of minimum nonlinear losses, maximum stationar- 
ity and localization. This scenario has been discussed in 
Sec. HV.Al where the important role of the laser wave- 
length in selecting the key player able to saturate the 
wave collapse and support the self-guiding process was 
underlined. For this purpose, Fig. [23la) illustrates the 
pulse flucncc along the self-guiding range measured at 
the exit plane of a 31-mm long cuvette for Pj n = 6.6 P CT . 
Fig. I2"37b) shows the A-shaped spectral-angular distribu- 
tion of the radiation. Plots (c,d) represent peak intensi- 
ties and FWHM diameters reached with similar pulse pa- 
rameters with and without plasma gain, which evidences 
the weak influence of plasma defocusing in this dynamics 
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strated the possibility to arrange filaments in water into 
ID arrays of parallel, non-overlapping spots by clipping 
the impinging laser beam by a slit aperture built from ra- 
zorblades. Metallic wire meshes can also be used to gen- 
erate space-controlled 2D arrays of filaments. Recently, 
Heck et al. demonstrated the efficiency of an adaptive 
control over the position and extent of filaments in water 
tanks, through a closed feedback loop setup employing a 
spatial light modulator and a genetic algorithm that al- 
low to manipulate the amplitude and phase of the input 
pulse. 
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FIG. 23 (a) Normalized 3D beam profile at the output of the 
water cell (31 mm); (b) Spectra l-angular distribution of the 
radiation ijDubietis et al l l2003h . Courtesy of A. Dubietis. 
(c) Peak intensities computed numrically with (solid curve) 
and without (dashed curve) plasma gain from 170-fs, 100- 
/im-waisted pulses focused with / = 5 cm (Pi n = 10P cr ) in 
water; (d) Corresponding FWHM radius. 



(jSkupin et all l2006ah . 

Increasing more the pulse power leads to multiple fil- 
amentation, which has been addresse d in a few papers 
(ICook et al. I, 120031 : iHeck et all l200fA ISchroeder etal\ . 
120041 1. Bv using a cylindrical planoconvex lens, Cook 
et al. produced horizontal arrays of stable white-light 
filaments in water at 800 nm, allowing interference ef- 
fects between neighboring cells. Similarly to a pair of 
Young's slits, a filament pair creates interference pat- 
terns, which is the signature for a constant phase re- 
lationship between the supercontinua generated by the 
filaments (the same property applies to atmospheric fila- 
ments; see Sec. IVII.A.4[) . Those thus appear as coherent 
sources of white-light. Besides, Schroeder et al. demon- 



Apart from water, other fluids can support fem- 
tosecond filament at ion, such as alcohols like ethanol or 
methanol. Dyes may be introduced into these liquids, 
in order to visualize the filamentary evolution through 
one, two or three-photon fluorescence and modify the 
multifilame ntation patte r ns by varying t he dy e con- 
centration dGuvon et all [20061 : iLiu et all l2005al . 120031 : 
ISchroeder and Chinl . 20041 ). Reiving on the dye nonlin- 
ear absorption, the structural changes of the filamcnta- 
tion becomes visible to the eyes, due to the fluorescence 
from the dye molecules excited by multiphoton excita- 
tions at visible wavelengths. From a dilute solution of 
methanol and Rhodamine B, Schroeder and Chin figured 
out that femtosecond filaments propagate straight, may 
die off prematurely or fuse into new spots. Along the 
propagation axis, "mature" filaments ending after the 
bright fluorescence zone were identified. By means of 
the same techniques, Liu et al. observed multi-focusing 
events by increasing the beam energy in methanol doped 
with 0.13% of Coumarin 440 at 800 nm. The photon 
bath surrounding the 20-/im large filament core was nu- 
merically examined: It was found that the near-axis re- 
gion takes energy from a ring-shaped region limited to 
r < 60 /im around the filament core. Outside, the periph- 
eral domain (r > 60 /im) acquires energy. The resulting 
energetic balance preserves an almost constant energy in 
the near-axis region. Nonlinear fluorescence techniques 
were recently used to discriminate between multiphoton 
absorption and co nical emiss i on by z-scan analysis us- 
ing metal meshes ( Liu et a l.. 2005a). It appears that 
the energy loss caused by conversion into CE (40%) is 
much higher than that caused by MPA (~ 3%) at powers 
< 8P CT (see also Sec. IIV.A|) . Finally, Guyon et al. em- 
ployed Coumarin 153 in dye-doped cells of ethanol. At 
high dye concentration (4 g/1), the filamentation pattern 
was shown to self-organize into a latticelike (hexagonal) 
figure, by letting the two-photon absorption of Coumarin 
switch out filamentary sites in the (a;, y) plane. Pump- 
dump experiments furthermore revealed that, at all fila- 
mentary sites, excited Coumarin molecules could coher- 
ently relax by fluorescence and emit in-phase light stim- 
ulated by the dump pulse, i.e., a collection of filaments 
can be used as microscopic laser sources in dense media. 
Figure [2"4l summarizes some of these observations. 
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FIG. 24 Filamentation pattern obtained from a 1-cm long cell 
of (a) pure ethanol, (b) dilute solution of ethanol/Coumarin 
153 at 4 g/1, (c) far-field fluence distribution of excited states 
of Coumarin emitting simultaneously in phase. 



VII. FILAMENTS IN THE ATMOSPHERE: CONVEYING 
INTENSE STRUCTURES OVER KILOMETERS 

The final section is devoted to the medium which orig- 
inally served as the "birthplace" for the science of fem- 
tosecond light filaments, namely, the atmosphere. Spe- 
cial emphasis is laid on experimental diagnostics. From 
the theoretical point of view, filamentation in air is mod- 
eled frojnl^cjs^jQS^^smgclassical dispersion relation in 
air ( Peck and Reedeii Il972t ). Although of weaker per- 
centage (20% vs 80% in air), dioxygen molecules have 
a lower ionization potential than nitrogen and they pro- 
vide the dominant spe cies prevailing through ionization 
(|Couairon et aLl . l2002h . 



A. Long-distance propagation and white-light 
supercontinuum 

1. Temporal chirping and spatial lensing 

The propagation of high-power (TW) femtosecond 
laser pulses in air has attracted considerable attention 
from the pioneering obser vation of the white- light super- 
continuum beyond 10 km (jWoste et aLl . ll997l ). A number 
of important practical ap plications have been s uggested, 
including remote sensing (R airoux et al. , 2000), directed 
energy delivery and artificial lightning (jLaFontaine et all 
1999a) among others. Because the filament onset and 
length are key parameters for spectroscopic measure- 
ments and for depositing high intensities on remote tar- 
gets, monitoring these parameters are of utmost im- 
portance. To achieve high intensities at remote dis- 
tances, a negati ve frequency chirp can be introduced i n 
the laser pulse (|Alexeev et all 12004 IWille et all l2002fh 
In addition to transverse self-focusing, the pulse under- 
goes a temporal compression as it compensates the nor- 
mal group- velocity dispersion along the propagation axis. 
Chirping effects can be measured by evaluating differ- 
ences in the conical emission w ith and without pulse 
chirping ( Rodriguez et all [2004) . For instance, the ver- 



tical propagation of the Tcramobile beam, presented in 
Appendix [C] was examined from the ground using the 2- 
m astronomical telescope of the Thiiringer Landesstern- 





FIG. 25 Typical fs beam image of the Teramobile laser beam 
from Tautenburg observatory, (a) Fundamental wavelength; 
(b) Blue-green band of the continuum. The horizontal stripes 
across the pictures come from stars passing through the tele- 
scope field of view. Note the strongly nonlinear altitude scale 
due to triangulation. 



warte (Thuringia State Observatory, Germany). Col- 
lected images were a combination of both (i) cross-section 
images of the beam impinged on the bottom of clouds or 
haze layers acting as screens, and (ii) side imaging of the 
Rayleigh-scattered light from the beam over large alti- 
tude ranges. Some of these images are shown in Fig. [25] 
at both the fundamental wavelength and the white-light 
continuum. In this case, the altitude was retrieved using 
triangulation. Figure 125T b) demonstrates the efficiency 
of the supercontinuum generation, since its blue spectral 
signal has been detected from altitudes beyond 18 km for 
appropriate chirping. These experiments demonstrated 
for the first time the possibility to deliver high-intensities 
and generate white light upon variable km distances, by 
means of GVD/chirp precompensation techniques. The 
same observation assessed that the conical emission bears 
two thirds of the overall white-light energy, while the on- 
axis, forward-directed central component carries the re- 
maining one third. 

To modify the filamentation distance, a chirp is usually 
introduced by changing the distance between the gratings 
of CPA laser compressors. Applied to the Gaussian pulse 
(|46|) . this technique modifies the second-order phase con- 
tribution ip" = d 2 ip/duj' 2 \ UJ=U j . It enhances the input 
pulse duration and decreases the beam power as 



t c p =t c p =°VTVc\ 



pC = pG=*/y/Y+C* (84) 



at equal spectral content. As it linearly propagates, the 
pulse competes with normal GVD and reaches the mini- 
mal duration t m i„ 



tp ~° at the distance 



\C\ 



1 + C 2 2k" ' 



(85) 



Pulse chirping mixes two essential modifications. The 
first one is dictated by the phase ~ e~ lCt '*p. For input 
powers above P cr , a positive chirp C > delays the oc- 
currence of the self- focus point, whereas this occurs ear- 
lier along the optical path for a negative chirp (C < 0). 
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The second one results in diminishing the effective pulse 
power, which pushes the first focus to later propagation 
distances and reduces the number of filaments in case 
of multifilamentation, whatever the sign of C may be. 
In usual experimental conditions, this second effect pre- 
vails over the first one. Despite the strong nonlinearities 
driving the filament dynamics, the chirped-induced lin- 
ear compression stage still persists with C < 0. It gives 
rise to additional focusing events and keeps the beam lo- 
calized at di s tance s close to z min ( Golubtsov et a/.l . l2003t 
iNuter et all , 120051 ). Moreover, the generation efficiency 
of the supercontinuum may vary by several orders of mag- 
nitude compared with that of a transform-limited pulse. 

The onset distance and longitudinal extent of fs fila- 
ments are also conditioned by the initial spatial focusing 
geometry. No simple analytical rule exists on this point, 
because of the complex spatial distortions destroying the 
initial homogeneity of the beam. The filament length 
moreover depends on the energy consumed by plasma 
excitation and on the accessible peak intensity. How- 
ever, most of the experiments emphasize the use of rather 
large ratios / /wq, in order to avoid an immediate plasma 
defocusing in tightly focused configuration. This prop- 
erty can be refound by integrating the dynamical equa- 
tions ([55)1 derived from the two-scale variational method. 
The variational principle indicates that at given Pi n /P cr 
the ratio f /zo must be large enough to insure the self- 
trapping condition z c < f. Otherwise, for f / Zq -C 1, the 
beam diffracts just after the focal point. Figure I^STa) 
shows the beam diameter of 100-fs, 1-mm-waisted pulses 
in air at P ln — 10 P cr and various focal lengths /, which 
confirms the previous belief. Figure l26l b) represents the 
beam extent along the z axis of temporally-chirped pulses 
with waist wq = 3 mm, input energy of 1.9 mJ and same 
power ratio, for which = 100 fs while t^ = 70 fs. 



2. Plasma and optical field measurements 

Plasma detection relics on the existence of a difference 
of potential produced by the current generated by ion- 
ization of air molecules. Working in free atmospheric 
medium makes it easy to measure the conductivity of 
the plasma channel, that causes a drastic reduction of 
air resistivity after the passage of a self-guided filament. 
This conducting column can directly be evidenced by 
letting the filament pass between two copper electrodes 
drilled in their center and between which a DC voltage 
of typically 1000 V is applied (see Fig. [27}. The cur- 
rent circulating through the plasma column is then mea- 
sured by record ing the voltage induced across an external 
load resistance (jTzortzakis et aLl . ll999ll2001d) . Knowing 
the current density per ion (i = 3 X 10~ 14 A/cm 2 ) and 
taking into account the volume occupied by a plasma 
filament, peak electron densities of 10 16 — 10 17 cm~ 3 
have been reported in air at 800 nm. Another tech- 
nique consists in resolving in time small local changes 
in the atmospheric refractive index by diffractometry 
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FIG. 26 (a) Beam diameter of pulses spatially focused with 
Wo — 1 mm, / = 35 cm (dashed curve), / = 37 cm (solid 
curve), / = 50 cm (dotted curve), / = 5 m (dash-dotted 
curve); (b) Beam diameter vs. z > 5 m for temporally- 
chirped, 3-mm-waisted pulses with C = (dash-dotted 
curve), C = 1.02 (dotted pulse) and C = —1.02 (solid curve). 



( Tzortzakis et ~aH l2000bl) . The principle is here to use 
a probe beam that crosses the filament path under a 
small angle. The far field image of the probe forms fringe 
patterns, that yield a direct measurement of the accu- 
mulated phase containing the plasma-induced defocusing 
effect. Let us also mention sonographic methods, that 
take advantage of the sound signals along the p lasma col- 



take advantage o t tne sound signals along tne p ig 
umn ( Hao et al\ , [2005b; Ho sseini et al. L l2004ri n 



_ Yu et all 

|2003( ). This " acoustic" diagnostics employs a microphone 
placed perpendicularly to the channel and recording the 
sound signals by a digital oscilloscope. The sound emit- 
ted from a plasma string is a portion of the pressure 
modulation of an acoustic wave from which the absorbed 
optical energy and plasma density are deduced from the 
electric peak voltage at a given distance z. By doing so, 
the variations in the electron density can then be plotted 
along the propagation axis. 

Plasma lengths attained from single ~ mJ femtosec- 
ond pulses currently remain of the order of the meter. 
The combination of twin pulses, launched collinearly in 
convergent geometry and separated from each other by 
a suitable time d elay, c an, however, double thi s length 
(|Couairon et all . 120031 ; iTzortzakis et all , |2003|) . This 
process, called "concatenation" of plasma filaments, re- 
lies on locking together the ionized channels generated 
by each individual pulse. By tuning their focal lenses 
(/2 — /i ~ 1 m) and their time separation (~ t p ), the less 
powerful time slice ending the self-guiding of the first 
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FIG. 27 Schemes of experimental setups for (a) conductivity 
measurements and (b) CCD imaging for optical field measure- 
ments. 



pulse coincides with the most powerf ul one of the second 
pulse starting a new plasma column ( Bergel . |2004[ ) . 

Experiments over a horizontal path provided more in- 
formation about the length over which ion ized plasma 
channels are formed (Mec hain et 0*1 l2005al) . The influ- 
ence of an initial negative chirp enlarging the pulse dura- 
tion from 0.2 to 9.6 ps, corresponding to 190 and 4 P CI at 
constant energy (190 mJ), respectively, was especially ex- 
amined upon distances up to 2350 m. Whereas ionization 
is clearly observed over propagation scales < 50 m for 
short pulse durations (0.2 ps), it becomes more and more 
sporadic for durations above 3 ps and ceases for 9.6-ps 
pulses, i.e., at powers becoming close to critical. In this 



case, a few low-intensity (~ 10 W/cm ), mm-waisted 
spots survive inside the photon bath, still capable of cov- 
ering several hundreds of meters. The bundle does not 
convey enough power to trigger full ionization. Instead, 
the beam evolves in a regime along which GVD is able to 
take over MPI and sustain the beam in a confined state 
over distances of the order of the input Rayleigh length 
( Champeaux and Bergel . [2005h . When the initial power 
is too weak, broad beams may produce a few bright spots 
by modulational instability, but none of these is capable 
of developing extensive plasma sequences (see Fig. H]). 

Several diagnostics exist for optical field measure- 
ments. The first consists in recording intensity profiles 
by a thick glass plate placed on the propagation axis at 
~ 45° angle (Fig. |2"T]) . The weak reflection from the glass 
is then imaged with a high-aperture l ens onto a linear 
charge coupled device (CCD) camera ( Tzortzakis et all 
l2001al) . Detection is performed out of the highest- 
intensity region close to the nonlinear focus, in order to 
leave the entrance window of the glass plate undamaged 
and keep up the reflected beam undistorted. Also, the 
high repetition rate of the laser source may not avoid 
multi-sp ot measurements. Fo r this reason, Bernstein 
et al. ( Bernstein et all l2003f ) reported on single-shot 



measurements of self-focusing pulses that do not have 
the intensity required to produce ionization. Besides 
a CCD camera, a spectrometer and second-harmonic 
frequency-resolved optical gating (FROG) device allowed 
to measure the spatial, spectral and temporal distribu- 
tions of the pulse, respectively. The data, collected from 
an initially collimated Gaussian beam, showed spatial 
and temporal narrowing and spectral broadening at dis- 
crete energy levels preceding the ionization stage. They 
confirmed a critical power value for nonlinear compres- 
sion effects of 11.5 GW for 800 nm pulses, i.e., about 
~ 10—15 GW. Nowadays, several diagnostics such as 
SPIDER and crossed (X)FROG traces complete stan- 
dard auto-correlation pictures and spectra to catch the 
spatio-temporal structure of a pulse in the (ui,t) plane. 
To explore the filamentation stage along which inten- 
sities as high as 5 x 1 13 W/cm 2 are attained, Ting , 
Gordo n and co-workers ( Gordon et all [2006t iTing et all 
I2005allbh elaborated on a new method following which the 
filament is propagated into a helium chamber through a 
nozzle that creates a sharp air-helium interface. Because 
helium has lower Kerr index and higher ionization po- 
tential than air, nonlinear focusing and plasma defocus- 
ing are arrested at the transition in the chamber. The 
filament then expands due to diffraction to larger sizes 
and lower intensities. A calibrated portion of energy can 
safely be collected and imaged either directly by a CCD 
array or through an imaging spectrometer. 



3. Multifi la mentation 

Open-air terawatt laser facilities make it possible to 
observe optical focal spots formed by a myriad of fila- 
mentary cells for input powers containing several thou- 
sands of critical powers in air. The 2D reduced model 
(|71[) using the experimental fluence as initial condition 
actually reproduces the evolution of the multifilamcnta- 
tion pattern over long scales, saving computational re- 
sources when one neglects the temporal dimension. As 
an example, Fig. [55] shows the overall envelope of the 
Teramobile bundle (see Appendix [C]) launched in parallel 
geometry. The initial pulse contains 700 critical powers 
and its transverse profile is scanned at different propa- 
gation distances. Filaments rise from the initial beam 
defects, form a crown of dots growing from the diffrac- 
tion ring and then excite clusters of cells. These clusters, 
some of whose are identified by the labels (1), (2) and 
(3), are faithfully reproduced by the numerics. The fil- 
ament number remains in the order of P m /Pfih where 
Pfii ~ 3 — 5F cr is the power in one filament. At large dis- 
tances, the primary brightest spots decay into secondary 
filaments by exchanging power throu gh the energy reser - 
voir formed by the background field ( Berge et all I2004T) . 

When the focusing geometry is changed, the filamenta- 
tion pattern becomes severely modified. Figure [2"9l shows 
a longitudinal visualization of the filaments nucleated in 
a beam with 760 P cr and focal length / = 40 m. The 
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FIG. 28 Filamentation patterns from the Teramobile beam 
with wo = 2.5 cm, FWHM duration of 100 fs and input power 
equal to 700 P CT at different distances z. (top) Experiments; 
(bottom) Numerical computations. 
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FIG. 30 Open cloud chamber. The cloud spans over 10 m. 
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FIG. 29 Teramobile fluence of a focused beam (/ = 40 m) 
with 760 critical powers, yielding three filamentary strings 
beyond the linear focal point: (top) Experiment, (bottom) 
numerical computation from the 2D model (|71[) . 



beam is directed towards an open cloud chamber scat- 
tering a weak-density cloud. After the focal point, only 
three strings of light emerge and cover about ~ 8 m 
each, while the same beam should produce more than 
one hundred filaments in parallel geometry. Numerical 
computations for this configuration clear up that many 
filaments are created before the linear focus, but they 
fuse into three strings of light acquiring a high directiv- 



ity afterwards (jSkupin et all l2004bl ) . The same experi- 
mental setup (Fig. [30]) put in evidence the robustness of 
femtosecond filaments through a multitude of 1 /im large 
water droplets randomly-distributed at various densities 
in the 10-m long cloud chamber. A filament is still trans- 
mitted through a cloud with an optical thickness as high 
as 3.2. For a cloud optical density of 1.2 or below, corre- 
sponding to cumulus or stratocumulus, the filamentation 
does not seem affected. Filaments remain visible at the 
exit of the fog even for a droplet concentration so high 
(8.6 x 10 4 cm -3 ) that each filament hits on the average 
2000 droplets per propagation meter. The corresponding 
extinction coefficient is 0.2 m -1 . Hence, filamentation 
can be transmitted through a fog over distances compa- 
rable with the visibility length. Energy losses due to ran- 
dom collisions implies an exponential decrease of power, 
which modifies the number and pos ition of the filaments 
in the bundle ( Meiean et all l2005h . Such experiments, 
whose results are recalled in Fig. [3TJ permitted to esti- 
mate the power per filament to about 5 P cr ~ 15 GW in 
air (see bottom panel). 

Pulse propagation in adverse weather becomes an im- 
portant topic for LIDAR applications. In this scope, 
Kandidov et al. addressed the point of the nucleation 
of filaments in a turbulent atmosphere s upporting statis- 
tical fluctuations of the refractive index ( Kandidov et all 
Il999l ). These fluctuations naturally arise as the medium 
becomes, e.g., locally heated. On the basis of a phase 
screen model, the cubic NLS equation including statisti- 
cal variations of the linear optical index (Kolmogorov tur- 
bulence) describes random paths for the nucleation and 
position of the nonlinear focus. This focus appears on the 
average shorter compared with unperturbed air and the 
beam centroid moves along a random path in the (x,y) 
plane. This induces a transverse deflection of the beam 
axis by several mm, whi ch was later observ ed from exper- 
imental averaged data ( Chin et a l, 2002b). Besides this 
beam wandering, femtosecond filaments keep their posi- 
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FIG. 31 Modification of a filamentation pattern with ~ 
120P C r over 50 m: (a) in dry air and (b) after traversing a 
10-m long tube of fog. The bottom panel reports the filament 
number vs. the transmitted power over critical. 



tion stable relatively to the bundle in highly turbulent 
regions, as the index variation caused by the nonlincar- 
ity exceeds turbulent-induce d refractive-index gradients 



k gr 

by two orders of magnitude ( Ackermann et a/ll2006bl) 



Furthermore, the influence of pressure variations on 
TW laser pulses was recently investigated. This issue is 
timely, as filaments propagating vertically over several 
km undergo an exponential decrease of the local pres- 
sure, p(z) = poe~ z / La , where L„ = 8 km is the at - 
tenuation length for air density (jSprangle et all l2002h . 
Pressure variations induce proportional changes in the 
GVD coefficient, the Kerr refractive index and molecule 
density available for ionization, as recalled by Eq. (|82p. 
At reduced pressure (0.7 atm), experiments performed 
at 3230-m altitude above sea level revealed farther fila- 
mentation onset and reduction of the fi lament number 
linked to the decrease of n 2 (p) by 30% ([Mechain et all 
l2005bl) . The filamentation process, especially the fila- 
ment length, was shown to be qualitatively unaffected. 
Numerical simulations in this field specified that a sin- 
gle atmospheric filament does not indeed significantly 
change its self-channeling range and keeps an almost 
constant peak intensity. In contrast, the plasma level, 
the filament width [Eqs. (fTB")) ] and the nonlinear focus 
([50")) evolve with p and y/p, respectively, due to the ef- 
fective increase of P CT (n2). Multifflamentation patterns 
have their onset distance governed by the maximum rate 
for modulation instability 7 max ~ 1/^-2 (p)-^o [ see Eq. 
(I7H)1. This length thus varies linearly with air pressure 



Champeaux and BergeLl2"b"06r) . 



4. White-light generation 

The white-light continuum generated in air by ultra- 
short laser pulses is essential in view of LIDAR applica- 
tions, since it constitutes the light source used in mul- 
tipollutant remote sensing. This white light was char- 
acterized over the last years, with progressively extend- 
ing bandwidths. The supe rcontinuum has firs t been 
characterized in the visible ( Alfano and ShapiroL Il970t 
iNishioka et all Il995h . J. Kasparian et al. next i nvesti- 
gated the infrared region (|Kasparian et a/.l . l2000bT ). Two 
different terawatt CPA laser systems (A: 60-mJ energy, 
35-fs minimal pulse duration, 25-mm FWHM beam di- 
ameter; B: 100-200 mJ, 100-fs minimal pulse duration, 
35-mm FWHM diameters) produced spectra measured 
at a total distance of ~ 30 m from the lens, as the laser 
beam was diffracting after 20 m of filament propagation. 
Illustrated in Fig. I3"27 a). the continuum band developed 
from laser system A is very broad, extending at least 
to 4.5 /im. An almost exponential decay over 4 orders of 
magnitude up to 2.5 /xm is observed, followed by a slower 
decay of one order of magnitude only. As shown from the 
inset plotting results from laser system B, variations in 
the input energy makes spectral intensity change by only 
one decade in some spectral regions. The spectral shape 
of these different pulses remains, nevertheless, quite sim- 
ilar within one decade in spectral intensity. Extension to 
the UV-visible domain down to 230 nm (Fig. [T7|) is not 
represented, as it was detected later. For comparison, 
Fig. 132T b) shows the numerically-computed spectrum of 
a 0.5-mm waisted, 4P cr ~ 10 GW, 100-fs pulse after a 
single femtosecond filament has been generated in air by 
accounting for third-harmonic generation in the limits 
T, T -1 -> 1 (|Berge et a/.l . l2005D and from the full UPPE 
model (|43p . The experimental and numerical spectral 
shapes are similar up to 1.2 /xm. Differences, however, 
occur at larger wavelengths, as space-time focusing and 
sclf-stcepening are included, which shortens the red parts 
of the spectrum to some extent. Recovering the experi- 
mental redshifts from the UPPE model in air is still an 
open issue. 

The similarity of white-light spectra emitted by TW 
laser pulses (experiments) and GW pulses suggests that 
the multifilamentation produces a spectrum analogous to 
that generated by a single filament. This conclusion is 
consistent with Chin et a/.'s experimental observations 
about the coherenc e properties of femtosecond filaments 
( Chin et aLl . [20~0 2a). Observation of the interference pat- 
tern produced by two or more filaments allowed to predict 
that they emerge in phase from the background field and 
they possess the same phase relationship. This means 
that the laser spectrum around the central wavelength 
in the conditions of multiple filamentation is in principle 
identical to that developed by an isolated filament. 

The previous property is important for LIDAR appli- 
cations, because it suggests that the comprehension of 
the spectral dynamics of one filament is sufficient for 
understanding the spectral dynamics of multifilamented 
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FIG. 32 (a) Measured white-light spectrum by 2-TW laser 
pulses (laser system A). Chirp settings correspond to 35-fs 
initial pulse duration without chirp (filled symbols) and 55-fs 
initial pulse duration with negative chirp (open symbols). In- 
set, spectrum measured from laser system B. The two curves 
have the same normalization factor, (b) Spectrum for a sin- 
gle filament numerically computed from Eqs. (|45p . The dash- 
dotted curve includes TH generation in the limits T, T _1 — > 1, 
while the solid curve represents the spectrum calculated from 
the complete equation set. 



beams. In addition to the wide spectral region covered by 
nonlinear femtosecond pulses, the strong enhancement of 
the backscattered photons in filamentation regime makes 
ultrashort laser pulses quite promising tools for the re- 
mote identification of multipollutants in aerosols. Indeed, 
the supercontinuum emitted by a filament is enhanced in 
the backward direction [see Eq. (|14p], i.e., towards the 
laser source, compa red with linear Rayleigh-Mie scatter- 
ing iYu et aZ.I . [200lh . 



B. Remote sensing (LIDAR) applications 

1. Principle of LIDAR: Towards " Femtolidars" 

The LIDAR (Light De tec tion A nd Rangin g) tech- 
nique dFuiii and Fukuchil. 120051: iMeasured . 11984 
iTheopold et all 120051: IWoli l2000f) was demonstrated 
shortly after the advent of the first lasers. In this 
technique, a laser pulse is emitted into the atmosphere. 
The backscattered light is collected on a telescope and 
detected as a function of time, with a typical resolution 
of 1-10 nanoseconds. This temporal window yields 
a high spatial resolution, since the flight time of the 
detected photons is directly proportional to the distance 



where they have been backscattered. Such a spatial 
resolution, combined with the possibility of sweeping 
the laser beam, provides two- and three-dimensional 
maps of measured atmospheric species. This is the 
main advantage of LIDAR over other measurement 
methods for atmospheric trace-gases. One of the most 
popular LIDAR techniques is called DIAL (Differential 
Absorbtion Lidar). It allows to selectively measure the 
concentration of gaseous pollutants by comparing the 
Lidar signals at two wavelengths nearby to one another, 
one being on an absorption line of the pollutant, and the 
other just beneath. However, this method is basically 
limited to pollutants that exhibit a narrow absorption 
line without interference from the absorption spectra of 
other atmospheric compounds. Moreover, the need to 
tune the laser wavelength on the absorption line forbids 
to simultaneously identify more than one pollutant 
within one acquisition. 

Femtosecond Lidars (or so-called " FemtoLidars" ) over- 
come classical DIAL limitations. Exploiting the "Fem- 
tosecon d atmospheri c lamp " discovered in 1997 by Woste 
et al. (jWoste et all [l997). high-power ultrashort laser 
pulses can be shined vertically into the sky, where they 
generate white-light through the filamentation mecha- 
nism. The backscattered light, recorded with a telescope 
linked to a time-gated spectrometer, then provides a fas- 
cinating vector for atmospheric resear c h ove r km ranges 
(|Rairoux et all \200$ , [Rodriguez et oil . l2004h . The spec- 
tral bandwidth developed by fs filaments is very broad, 
since it spans at least from 230 nm in the ultraviolet to 
4.5 /im in the mid-infrared (Figure [3"3"|) . On the spectrum, 
the absorption band of water between 1.8 and 2.5 /xm is 
clearly visible, showing the potential of white-light for 
optical remote sensing in the atmosphere. It also cov- 
ers the absorption band of volatile organic compounds 
(VOCs) between 3 and 3.5 fim. VOCs constitute a family 
of organic compounds with strongly overlapping absorp- 
tion spectra, which prevents any measurement by clas- 
sical DIAL techniques. In addition, the flat continuum 
spanning from the visible down to 230 nm through third- 
harmonic generation (see Fig. [17)) provides a promising 
light source for the measurement of trace gases that ab- 
sorb in the blue or the UV, such as ozone, toluene, ben- 
zene, SO2 or nitrogen oxydes. Recently, femtosecond 
laser-induced filaments were also appl i ed for g reenhouse 
gaz methane (CH4) in air (|Xu et all l2006al) . Intense 
filaments dissociate pollutant molecules into small frag- 
ments, which emit characteristic fluorescence. This can 
be used to remotely measure the pollutant concentration 
at characteristic spectral lines. 

Besides, multiparameter measurements are neces- 
sary to monitor the dynamics of atmospheric physico- 
chemistry or to remotely analyze a cocktail of species 
emitted during a chemical accident. They are also re- 
quested to characterize the nucleation and maturation 
of clouds. Such processes play an important role for at- 
mospheric modeling, both on meteorological and clima- 
tological scales. In particular, the droplet growth and 
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FIG. 33 Spectrum of the white light continuum assembled 
from five spectral regions. Main atmospheric absorption 
bands are marked above the curve. 



related density have a key influence on the forecast of 
both precipitations and earth albedo. Their characteri- 
zation requires continuous measurements of the size dis- 
tribution inside the clouds, with a temporal resolution of 
a few tens of minutes, compatible with the evaporation 
and growth time ranges. While airborne measurements 
are too expensive for routine monitoring, radiosound- 
ing does not provide the adequate repetition of prob- 
ing. Therefore, optical techniques are promising, and 
they can be made optimal with fs laser pulses. As an 
example, a high-resolution absorption spectrum over a 
spectral interval larger than 200 nm is shown in Fig. 
I34[ based on the white-light continuum, providing both 
the water vapor concentration and the tempera ture of 
air ( Bouravou et all l2005t iKasparian et al. , 2003). Once 
combined, these two parameters supply the relative hu- 
midity, which is the relevant factor for atmospheric dy- 
namics. Moreover, angular measurements of the multi- 
ple scattering yield the size distribution and concentra- 
tion within the clouds, through an inversion of the laws 
of multiple Mie scattering based on a genetic algorithm. 
These measurements require the same laser source and 
two independent detectors. With fs Lidar setups, they 
can be implemented simultaneously so as to yield a com- 
plete characterization of the cloud microphysics. 

" FemtoLidars" are not restricted to the white-light Li- 
dar. The high intensities carried in the filaments can 
generate nonlinear effects in situ on a target. Such non- 
linear processes, which include ablation, ionization and 
multiphoton-excited fluorescence (M-PEF), constitute a 
supplementary information channel for remote sensing. 
An example of this technique has been provided by the 
remote detection and i dentification of biolog ical aerosols 
through 2-PEF Lidar (jMeiean et all 120041 ). Here, the 
purpose was to detect and locate rapidly a suspect emis- 
sion, to map the emitted cloud, and to identify poten- 
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FIG. 34 High-resolution atmospheric absorption spectrum 
from an altitude of 4.5 km measured in a Lidar configura- 
tion. The broad spectrum allowed to simultaneously measure 
the air temperature and humidity. 
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FIG. 35 Remote detection and identification of bioaerosols. 
The femtosecond laser illuminates a plume of riboflavin 
(RBF) containing microparticles 45 m away (left). The back- 
ward emitted 2-photon excited fluorescence (2-PEF) , recorded 
as a function of distance and wavelength, exhibits the specific 
RBF fluorescence signature for the bioaerosols (middle) but 
not for haze (right). 



tially pathogen agents among the multiple background 
atmospheric aerosols, including organic compounds such 
as soot or pollens. The 2-PEF allowed, for the first time, 
to remotely identify aerosols simulating biological agents 
in air by nonlinear Lidar (Fig. [35)) . Moreover, collected 
signals based on N-PEF become more efficient than Lidar 
signals based on f-PEF for distances above a few kilome- 
ters. This is due to two effects: (i) the directional emis- 
sion of the N-PEF from the aerosol particles minimizes 
the decrea se of the collection efficiency with increasing 
distances ( Boutou et ali , l2002j ). and (ii) the visible or 
near-infrared wavelengths used for exciting N-PEF, com- 
pared to the UV wavelength required for 1-PEF, experi- 
ences less attenuation in the atmosphere because of the 
1/A 4 dependence of the Rayleigh scattering. 



2. Remote filament-induced breakdown spectroscopy 

The ability of the filaments to remotely deliver in- 
tensities as high as 50 TW/cm 2 also opens the way 
to innovative exploration techniques, such as the "re- 
mote fi lament-induced breakdow n spectroscopy" (R- 
F1BS) (jStermaszczyk et all I2004D . This consists of a 
combination of Lidar and LIB spectroscopy (LIBS). LIBS 
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(ICremers and Kniehtl . l200Ct [C remers and Radziemskil . 
120061) is a versatile tool allowing an e lemental analysis o f 
surfaces of materials such as metals ( Angel et'all . l200lh , 
lastics (Jong-Il et all [20021: iNiessneii 1 1994T) . minerals 



Spectrometer + iCCD 



Knig ht et a/I l2000l : ISharma et aU , 120031: IWiens et al. 



2002), aerosols, biological materials (jKvuseok et 



19971 ) or liquids as well. It relies on the local ioniza- 



tion of the surface by a strongly focused pulsed laser, 
typically a Nd:YAG laser. The emission spectrum of the 
plasma generated at the surface allows a fast analysis, 
either qualitative or quantitative, with detection limits 
down to a few parts per million (ppm) for some elements. 
The use of subpicosecond laser pulses significantly en- 
hances the reproducibility of the measurements, because 
the lower pul s e energy limits th e heat i ng of the sample 
(lAlbert et all 120031: lAngel et all 120011 iDou et all l2003t 
iRohwetter et aZ.1 . 120031 ). The use of broadband detection 
systems makes LIBS a flexible technique which requires 
neither preparation of the sample, nor a priori knowl- 
edge of the elements to be found. Applications such as 
the identification of highly radioactive nuclear waste or 
real-time monitoring of melted alloys in industrial pro- 
cesses require a remote analysis technique. LIBS, which 
is suitable for raw samples, is a good candidate in this 
regard, since it only requires a direct view to the sample. 
However, due to the limited size of the optical compo- 
nents, diffraction intrinsically limits the intensity that 
can be focused on a remote target. On the contrary, self- 
guided filaments can deliver much higher intensities than 
the ablation threshold of many species, at distances of 
hundreds of meters or even kilometers. So, in R-FIBS 
techniques, the laser-generated filaments are launched 
on a remote target, and the light emitted by the ex- 
cited plasma plume is collected through a detection setup 
comparable to a Lidar-based one (Fig. |3"6T). Although is - 
sued from 180 m distances ( Stelmaszczvk et all . 12004) . 
the data suggest that mea surements can be dimen sioned 
up to the kilometer range (|Rohwetter et all [2003) . This 
technique is currently developed for various applications , 
e.g., the mo nitoring of heritage (TTzortzakis et all 20061 ) 
or bacteria (jBaudelet et all 120061 : IXu et all l2006b| ). 



C. Towards a laser lightning rod 

Besides remote sensing applications, ultrashort fila- 
ments may give access to the control of lightning strikes. 
Classical techniques to trigger lightning have employed 
rocket-pulled wires since the 1970's. However, the num- 
ber of rockets available per storm is necessarily lim- 
ited and rockets must be launched synchronized with 
the raise of the ambient electric field. Moreover, the 
wire falling down may pollute the measurement as well 
as the environment. Therefore, the idea emerged to 
apply lasers to control lightning by ionizing the am- 
bient air along the beam and forming a conducting 
plasma "wi r e". First attempts used nanosecond pulses 
(|Balll . 11974 iKoopman and WilkersonI Il97lt iMiki et all 



20 to 90 m 



Sample 



0.85 m 




fs-laser 
(3 TW - 80 fs) 



Typical filamentation pattern 



filaments 




FIG. 36 Principle of Remote filament-induced breakdown 
spectroscopy (R-FIBS). Bottom: beam profile near to the 
sample showing multifilamentation with typ ically 30 filaments 
across the beam (|Stelmaszczvk et all 12004 ). 



1993). They were unsuccessful, because such lasers 
could not produce continuously-ionized plasma chan- 
nels. More recently, this field was renewed by the ad- 
vent of CPA lasers providing higher intensities in shorter 
pulses, therefore avoiding laser absorption by inverse 
bremsstrahlung. Encouraging results have been obtained 
using focused ultrashort laser pulses to trigger and guide 
high- voltage discharges over several meters in the labora- 
tory (IComtois et aZ.Ll2000t iPepin et all \200ll iTing et all 
l2005al ). Others have been obtained on smalle r scales 
using UV ultrashort lasers ( Rambo et al .1 120011) . Since 



the filaments provide a conducting path over several me- 
ters or even longer, they are particularly suitable for 
the atmosphere. Spectacular experiments with the Ter- 
amobile laser installed in a high-voltage facility showed 
that ultras hort filaments can guide discharges over up 
to 4.5 m ( Ackermann et al\ , 1200 6a). Instead of their 
usual erratic path, discharges are guided along the trig- 
gering laser beam (Fig. |3"T|) . Moreove r, the breakdown 
volta ge is typically reduced by 30% ( Rodriguez et all 
I2002IJ . Partly guided discharges also occur in sphere- 
plane gaps, providing valuable information about the 
mechanism of the initiation and propagation of laser- 
trigg ered streamers, wi t h, e.g., plasma lifetimes of about 
1/is ( Ackermann et all , l2006al) . Furthermore, an artifi- 
cial rain does not preve nt the laser filaments fro m trig- 
gering such discharges ( Ackermann et all l2004h . Cur- 
rent research now focuses on the possibility to extend the 
plasma lifetime, in order to increase the guiding length 
and improve scalability to the atmosphere. This ap- 
proach relies on re-heating and photodctaching electrons 
of the plasm a channel by subsequ ent pulses, either in the 
nanosec ond (iRambo et all. 120011) or in the femtosecond 
regime ( Hao et a7] . l2005a ). Although high laser powers 
are usually believed to efficiently detach electrons from 
O^ ions in the plasma, it was recently demonstrated 
that a subsequent Nd:YAG laser pulse of moderate en- 
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FIG. 37 Laser control of high-voltage discharges. (A) Exper- 
imental setup. (B) Free discharge over 3 m, without laser fil- 
aments. Note the err atic path. (C) Straight discharge guided 
along laser filaments (iKasparian el all 120031 ) . 



ergy (sub- Joule) at 532 nm efficiently supports the trig- 
gering of discharges by an infrared, femtosecond laser 
(jMeiean et qfl . l2006al) . This effect was interpreted as re- 
sulting from a positive retroaction loop where Joule heat- 
ing of the plasma channel enhances photodetachment, 
while the resulting higher electron density boosts in turn 
the Joule effect. 



VIII. OUTLOOK 

Femtosecond lasers opened up beam parameter ranges 
from which universal, robust ultrashort filamentary 
structures conveying high intensities can propagate over 
long distances in a rich variety of transparent media. 
These " femtosecond filaments" are characterized by sizes 
of a few tens up to the hundred of microns and can reach 
ultrashort FWHM durations down to the optical cycle 
limit. They design new objects, that can be exploited for 
delivering few-cycle pulses at very high power levels in 
the future. Ionization processes do not break this prop- 



erty. They even amplify it by cutting the back time slices 
of the pulse. An accurate control of plasma generation 
should push the frontiers in nonlinear optics and allow 
for a precise monitoring of high-order harmonics. The 
basic scenario of a dynamical balance between Kerr fo- 
cusing and plasma defocusing explains the existence and 
the long life of femtosecond filaments in most of propa- 
gation configurations. This scenario must, however, be 
revised at visible wavelengths in condensed materials, in 
which chromatic dispersion and nonlinear absorption (at 
high enough power) seem to drive the self-guiding mech- 
anism through an important conical emission managed 
by AT-shaped waveforms in normally dispersive regimes. 
This aspect requires more investigations. Also, propa- 
gation regimes of anomalous dispersion, which promote 
temporal compression and may further improve self- 
compression techniques, should be deepened. 

In the last years, impressive progresses have been made 
in the understanding of the nonlinear propagation of 
high-power laser pulses over long distances in the atmo- 
sphere. For broad beams, multiple filamentation becomes 
more and more controllable by the use of dcformable 
mirrors and genetic algorithms. Besides classical tech- 
niques of linear focusing and phase chirping, this will 
soon enable experimentalists to optimize the propagation 
range of very powerful pulses. Filaments can therefore be 
used to remotely deliver high intensities generating in- 
situ nonlinear effects. The spectrum of the white-light 
continuum is much broader than that expected a few 
years ago, especially in the ultraviolet where the mix- 
ing between the fundamental and third harmonics wave- 
lengths together with the spectral amplification caused 
by pulse steepening give rise to a spectral plateau ex- 
tending down to 230 nm. These progresses open the 
way to applications such as remote sensing of gaseous 
pollutants and aerosols by Lidar, LIBS analysis of solid 
samples, or the control of high- voltage electric discharges 
and lightning strikes. Due to its mobility, the Teramobile 
allowed to demonstrate the feasibility of many of these 
techniques. These applications will probably be pushed 
forward in the future by technological improvements in- 
volving more compact systems, diode pumping, as well 
as spatial and temporal pulse shaping. Novel active me- 
dia, such as Ytterbium doping or OPCPA, also open the 
way to new spectral regions, especially the infrared. This 
spectral region, where eye-safety constraints are easier to 
fulfill, is nearer the absorption region of many pollutants, 
such as volatile organic compounds. These more flexible 
systems will be easier to operate, allowing routine uses 
for future industrial or atmospheric applications. 
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APPENDIX A: Ionization rates for atoms and molecules 

This appendix details different photo-ionization theo- 
ries for gases and dense transparent materials. 



1. Ionization in gases 



over the spatial coordinates result in the expression 
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To model the generation of free electrons by an in- 
tense light field interacting with a gaseous medium, 
photo-ionization theor ies were f i rst de veloped for atomic 
systems by Keldysh ()Keldvshl . 1 19651 ), later by Perelo- 
mov. Popo v and Terent'ev (PPT) (|Perelomov et all . 
Il966l . Il967f ). before being refined by A mmosov, Delone 
and K rainov (ADK) twenty years after ( Ammosov et all . 
198G). More recently, Tong et al. proposed a the- 
ory for molecular species in tunnel ionization regime 
(iTong et qJ.l . l2002r i. 



The function L(p, t) is periodic, with period equal to T = 
2ix j 'to. It can then be decomposed into Fourier series. 
Plugging this series into Eq. (IA4|) leads to 
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with 



a. The Keldysh theory 

Keldysh's theory is limited to hydrogenoid atoms in 
their fundamental electronic state and does not consider 
the Coulomb interaction between the leaving electron and 
the residual ion. The ionization rate W of an atom ir- 
radiated by a laser field E(t) = E p cos(u)t) is evaluated 
by 
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where q e is the electron charge and the ground state of 
the hydrogenlike atom, characterized by the energy E g — 
— Ui (ionization potential), is 
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Here, a = clb/Z includes the Bohr radius clb of hydrogen. 
The continuum electronic states are described with the 
Volkov functions 

where A(t) is the vector potential of the laser field. In- 
sertion of (|A2p and (|A3|) into Eq. (|A1|) and integration 
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Assuming that the electron leaves the atom with a small 

2 

kinetic energy Cff;), the poles in the denominator 
of L{p) are written as 
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where 9 is the angle between the impulsion vector p and 
the electric field E{t). 7 = ujy/2m e Ui/ (\q e \E p ) is the 
adiabaticity Keldysh parameter. It involves the ratio 
of the ionization potential over the ponderomotive en- 
ergy U p = q 2 E 2 /4:m e u! 2 . By means of the Saddle point 
method and the residue theorem, Eq. (|A7[) is next inte- 
grated to yield 
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where Ui = Ui + U p and S(j, x) is defined by 
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To take electron-ion correlation into account, Keldysh 
eventually multiplies the ionization rate by the factor 
Ex pressed in atomic units (a.u . ) 



(C7 i /M7/^/^^ : 7 ^ • 

\q e \ = h = cib = 1 (|Bransden and Joachainl [2003T) ■ 
the resulting ionization rate expresses 
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where E = {2U l f' 2 , 7 = lu^/2~U~/E p , v = U l /[?uj) t 

(3 = 27/VI + 7 2 , a = 2 [sinh- 1 ( 7 )- 7 / v /rTl 

Vq =< v + 1 > and $ m (x) = e~ x2 f£ \x 2 - y 2 )^e y2 dy. 
Equation ()A1 1[) differs from the original Keldysh formu- 
lation by a factor 4, originating from a corrected version 
of the residue theorem. 

This theory was the first one able to describe atom 
ionization by an alternating field in low intensity regime 
(7 3> 1) as well as in high intensity regime (7 <C 1). The 
former regime refers to multiphoton ionization (MPI), 
through which the electron is freed as the atom absorbs 
K =< Ui/ (hio) + \ > photons. The latter one corre- 
sponds to the tunnel regime, for which the electron leaves 
the ion by passing through the Coulomb barrier. In MPI 
regime, the ionization rate is obtained by taking the limit 
7 — > +00 in Eq. (|A11|) . which reduces to 

W = (7 {K) x I K 7 (A12) 

where I is the laser intensity and is the photo- 

ionization cross-section 
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Expressed in s 1 cm 2K /W K , the above "cross- 
section" parameter must be converted as cr( K ) — > 
CT W[a.u.]/[2.42 x 1(T 17 x (3.51 x lO 16 )^]. 



accurate model. First, they included the Coulomb 
interaction between the ion and the electron, when the 
latter leaves the atomic core. Second, they considered 
any atomic bound states as initial. The resulting rate is 
then 
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where n* = Zj \J2Ui is the effective quantum number, Z 
is the residual ion charge, I* = n* — 1 and n, I, m are the 
principal quantum number, the orbital momentum and 
the magnetic quantum number, respectively. The factors 
|C n «,;«| and f(l,m) are 
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Even if Eq. (|A14[) is usually presented as the PPT for- 
mula, the coefficients |Cn*,i*| are m fact extracted from 
the tunne ling theory derived by Ammosov, Dclone and 
Krainov ( Ammosov et all Q.986). Differences between 
PPT and ADK coefficients essentially lie in the fact 
that ADK theory employs electron wavefunctions in a 
Coulomb potential (Volkov states), which are connected 
by continuity with the continuum states at large dis- 
tances (r > l/y/Wl). 



c. The ADK molecular theory 



The PPT rate [Eq. l[AT4]l ] holds to describe photo- 
ionization of atoms. It can lead to some discrepancy 
when it is applied to molecular systems, because the co- 
efficients |C„*,i*|, originally evaluated from atomic wave- 
functions, cannot reproduce molecular peculiarities, such 
as, for example, the suppression of ioniz ation observed 
from the molecule 2 (jDeWitt et aLl . l200ll ). To overcome 
such limitations, we may extend the molecu lar tunneling 
theory by Tong et al. ( Tone et~al\ , l2002h by plugging 
molecular coefficients into the tunnel limit of the PPT 
formula and prolonging the latter to low intensity MPI 
regimes analytically. By doing so, Eq. (|A14p is able to 
describe molecule ionization after the substitution 



b. The PPT theory 

Later, Pere lo mov, Pop ov and Terent'ev 



( Perelomov et al .], 119661 I1967D developed a more 
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where the coefficients Ci have been established for differ- 
ent molecules. For dioxygen, C2 = 0.683 and C 4 = 0.033 
while Ci=C 3 =0, and m = 1 (|Tong et aZ.l . l2002h . 



2. Ionization in dense media 

Plasma generation in dense media is described with 
the ionization rate for crystals developed by Keldysh 
( Keldvshl . 1 19651 ). Its analytical evaluation is identical to 
that applying to atoms, except that the initial states are 
now modeled by Bloch wave functions. Following similar 
procedure, the ionization rate for crystals with energy 
gap Eg irradiated by an electromagnetic field E p cos(wt) 
is then 
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APPENDIX B: Atomic dipole for High-Harmonic 
Generation 

The nonlinear polarization vector PniX't, t) used to de- 
scribe HHG in gases is expressed as 



PNh{r,t) = p nt d(r,t), 



(Bl) 



where d(r, t) is the time-depe ndent atomic dipole ca l- 
culated by Lewenstein et al. ( Lewenstein et all . I1994T) . 
For a single electron interacting with a linearly-polarized 
laser field E(t) = E p xcos(ujt), the atomic dipole is de- 
fined by 



d(f,t) = q e (i;(f,t)\f\i;{f,t)). 



(B2) 



Solving the time-dependent Schrodinger equation, vj)(r, t) 
denotes the electronic wavefunction 
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Here, the functions K(x) = j^ /2 (l - x sin 2 B)- 1 l' 1 dQ 

and £{x) = ^^(1 — xsin 2 9) l / 2 d9 are the com- 
plete elliptic integrals of the first and second kind 
( Abramovitz and Stegunl . Il972t ) and to* is the reduced 
mass for the electron/hole pair. The above equation cor- 
rects a slip of pen occurring in the original Keldysh's 
formula. Note that whereas W for gas [Eq. JA11[) ] is ex- 
pressed per time unit, W for condensed media [Eq. (IA17I) ] 
is expressed per time unit and per cubic meter. This ion- 
ization rate reduces at low intensities to its multiphoton 
limit (7 — > +00) taking the form 
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and no — n(uj — loq) is the linear refractive index. 



where Ui is the ionization potential of the atom, |0 > 
is the ground state and b(v, t) is the amplitude of the 
continuum states \v >, respectively. In Eq. (|B3j) . several 
assumptions have been made, which limits the validity 
domain of this theory: 

• All excited electronic states are ignored, reducing 
the model to the harmonic production of orders 
2K + 1 > &. 

• The depiction of the ground state is neglected such 
that the ground state amplitude is equal to 1. 

• Electrons in the continuum are treated as free par- 
ticles moving in an oscillating laser field, with no 
Coulomb potential. This assumption is valid when 
the ponderomotive energy U p is higher than the 
ionization potential. 



Inserting (|B3p into Eq. (|B2|) . we obtain 
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where d(p + q e A(t)) is the field-free dipole transition ma- 
trix element between the ground state and the continuum 
state characterized by the momentum v = p + q e A(t). 
Here, p is the canonical momentum and A(t) is the vec- 
tor potential of the laser field [E{t) — —dA{t)/dt]. Its 
formulation for transition from state lis) is 
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where m e is the electron mass. S(p, t, if) is the quasi- 
classical action describing the motion of an electron freely 
moving in the laser field with constant momentum p as 
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The harmonic amplitude cfeif+i follows from Fourier 
transforming the time-dependent dipole moment d(t) 
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For isotropic media, only odd harmonics are produced. 
The dominant contributions to cfoif+i comes from the 
stationary points of the Legendre-transformed quasi- 
classical action, for which the derivatives of S(p,t,t') — 
(2K + \)hw x tot with respect to p, t and if vanish. In- 
troducing the returning time r = t — if ', the Saddle point 



equations read 



Vp-S 1 = pstLJT + q e I A(t") d(ut") = (B8a) 

"(t-T) 



dS ^ | [p st +q e A(t-r)} 2 =Q 



d(ijjr) 2m e 

dS _ [p st + q e A(t)} 2 [p st + q e A(t - t)} 2 
d(uji) 2m e 2m e 
= (2K + l)hu 



(B8b) 



Sc) 



Equation (|B8aP reduces to VpS(p,t,t r ) = x(t) - x(t') — 

where x(t) = J Q v(t')dt'. This relation imposes that 
the electron trajectories return at time t to the same 
point they left at the time if of ionization. Using Eqs. 
(|B8bjl and (|B8cj) together with the Fourier expansion of 

[(pst + qeA(t-T)) 2 1 '(m e hw) + 2U l j '(huj)]^ 3 yields the final 
expression for the Fourier component of the dipole, i.e., 



12K+1 



\hu J IF V Hlo ) 

U P C(^) 



E 



r+oo 


7T 


/ d$ 

10 


i ■ * 



6 M ($) e - £ ^z^- M 



B($)J A '-Af+2 



iB($)e l<p J K _ M _ 1 ( 



U P C{§) , 



+i [B($)e 4 * + £>($)] J J 



K-M+l 



U P C(<5>) 

ftuJ 



)+ [B($)+^($)e 1 *] J K -m{ 



(B9) 



where Fk{x) = (Ui + U p — Kftw)x — 2U P [1 — cos(x))/x, 
B(x) = -| + sin(o:)/a; - 2 sin 2 {x/2)/x 2 and D(x) = 
— 2B(x) — 1 + cos(ie). The function 6^,f ($) has a cumber- 
some expression, whose analytical meth ods to compute 
it can be found in (jAntoine et all , [19961) . 



APPENDIX C: The Teramobile laser 

Field experiments are required to characterize the fila- 
mentation over long distances as well as to develop atmo- 
spheric applications in real scale. Such experiments de- 
mand mobility to perform investigations at adequate lo- 
cations. Studies of high-power fs-laser beam propagation 
over km-range distances can only be performed outdoors, 
where relevant aerosol pollutants take place, e.g., in ur- 
ban areas or at industrial sites. Laser-induced lightning 
investigations require spots where the lightning proba- 
bility is high, as well as test experiments at high voltage 
facilities. These considerations clearly define the need 
for a mobile fs-TW laser system, embedded in a standard 
freight container-integrated laboratory equipped with the 
necessary Lidar detection, power and cooling supplies, 
temperature stabilization, vibration control, and an addi- 



tional standard Lidar system to assure eye safety. These 
speci fications were first achieved by the Teramobile sys- 
tem ()Wille et all l2002h . The laser itself is based on a 
TkSapphire CPA oscillator and a Nd:YAG pumped Ti:Sa 
amplification chain. It provides 350 mJ pulses with 70 fs 
duration resulting in a peak power of 5 TW at around 800 
nm and with a repetition rate of 10 Hz. Its integration 
in the reduced space of the mobile laboratory required 
a particularly compact design (Fig. [35]). The classical 
compressor setup has been improved into a chirp gener- 
ator to pre-compensate the group velocity dispersion in 
air. Combined with an adjustable focus, this permits to 
control the location of the onset of filamentation and its 
length. Mechanical and thermal stabilities of the mobile 
laboratory are kept under control, so that the Teramo- 
bile can be transported to any place in the world and 
operated even under adverse weather conditions. 

The Teramobile container also includes a Lidar detec- 
tion chain based on a 40 cm receiving telescope, a high- 
resolution spectrometer equipped with a set of gratings 
and detectors allowing simultaneous temporal and spec- 
tral analysis of the return signal in a wavelength range 
comprised between 190 nm and 2.5 /zm. 
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FIG. 38 Three-dimensional view of the Teramobile. (L): 
Laser system: Ti:Sa oscillator and its Nd:YAG pump laser 
(LI), stretcher (L2), regenerative amplifier, multipass pream- 
plifier (L3) and their Nd:YAG pump laser (L4); Multipass 
main amplifier (L5) pumped by two Nd:YAG units (L6); Com- 
pressor (L7). (S) Beam expanding system; (C) Power sup- 
plies; (D) Lidar detection system. 
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